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Fuel-air  explosions  (FAE)  have  been  investigated  in  the 
context  of  far-field  nuclear  airblast  simulation »  The  objec¬ 
tive  of  the  investigation  is  to  determine  the  feasibility  of  a 
reusable  FAE  simulator  at  the  one  kiioton  level.  Two  issues 
have  been  researched  in  parallel  efforts.  These  are  the  mech¬ 
anisms  by  which  largescale  FAE  clouds  of  controlled  shape  can 
be  reliably  and  repeatedly  formed  and  detonated,  and  the  quality 
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20.  ABSTRACT  (Continued) 

of  nuclear  airblast  simulation  that  is  achieved  when  such  FAE 
clouds  are  detonated. 

The  formation  of  hemispherical  clouds  by  simultaneous,  impulsive 
liquid  fuel  injection  through  a  large  number  of  radially  directed, 
centrally  clustered  nozzles  is  discussed  in  detail.  Specific  fuel 
dispenser  designs  are  considered  and  use  experience  with  two  of 
these  is  described.  Survey  experiments  dealing  with  the  atomization 
and  penetration  characteristics  of  large-diameter,  impulsively 
formed  single  liquid  jets  are  discussed.  Small-scale  hemispherical 
cloud  formation  experiments  are  also  discussed. 

A  new  reusable  facility  for  testing  FAE  simulation  of  nuclear 
airblasts  at  the  1/4- ton  scale  is  described.  Preliminary  surface- 
burst  experiments  at  that  scale  using  propylene  oxide  and  heptane 
as  fuels  are  discussed.  The  results  from  these  experiments  are 
scaled  and  compared  with  1  KT  nuclear  curvefits.  The  agreement  is 
judged  to  be  reasonable  at  this  scale. 

Numerical  calculations  of  the  airblast  that  emerges  from  a 
detonated  heptane-air  cloud  have  also  been  carried  out.  .  The  calcu¬ 
lated  results  are  compared  with  the  nuclear  curves.  The  quality  of 
simulation  at  the  1  KT  level  is  generally  good  and  becomes  excellent 
at  ranges  corresponding  to  overpressures  below  about  70  kPa  (10  psi) 
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SUMMARY 


The  progress  of  an  experimental  and  theoretical  research 
program  is  reported  in  which  simulation  of  fat-field  nuclear 
airblast  by  means  of  fuel-air  explosions  ( FAE )  is  being  inves- 
tigatea.  The  advantages  of  a  full-scale  reusable  FAE  blast 
simulator  include  the  absence  of  cratering,  ejecta,  and  signi¬ 
ficant  grouna  shock,  short  turn-around  times  between  blastweve 
experiments,  and  relatively  lower  costs  per  experiment  when 
compared  with  other  means  of  blast  simulation.  The  existence 
of  such  a  simulator  should  greatly  enhance  the  state-of-the-a r t 
of  blastwave  simulation  and  provide  a  means  for  accelerating 
our  knowledge  of  blastwa ve-structur al  interactions. 

The  present  investiga tipn  has  focused  specifically  on 
explosions  of  hemispherical  fuel-ait  clouds  formed  from  a  point 
source,  using  multiple-nozzle  liquid  fuel  injection.  The 
research  involves  two  areas.  First  is  an  assessment  of  the 
degree  of  correspondence  between  FAE  and  nuclear  axrblasts  in 
terms  of  peak  and  time-resoived  static  and  stagnation  pres¬ 
sures,  static  impulse,  positive-phase  duration,  and  wave-front 
decay  rates.  A  small-scale  instrumented  FAE  facility  (cloud 
diaaieter  ^9  m  (30  ft))  has  been  developed  for  this  purpose. 

The  second  research  area  involves  examining  the  engineering 
requirements  for  scale-up  of  this  facility  to  a  one  kiloton 
nuclear  equivalent.  In  particular,  impulsive  liquid  fuel 
injection  has  been  investigated  with  nozles  ranging  from  2  to 
10  cm  (0.75-4  in.)  in  diameter.  The  ultimate  vertical  reach, 
degree  of  atomization,  fuel-air  distribution,  anu  transverse 
spreading  rate  of  the  jets  ate  considered  to  be  tne  basic 
dependent  variables,  with  initial  jet  velocity,  nozzle  diam¬ 
eter,  nozzle  length,  quantity  of  fuel  dispensed  per  nozzle,  and 
fuel  properties  (uensity,  viscosity,  surface  tension,  and  vapor 
pressure)  taken  as  independent  variables.  Current  estimates 
inaicate  that  a  one  Kiloton  FAE  simulation  using  propylene 
oxide  would  require  a  hemisphe r ica 1  cloud  142  m  (466  ft)  in 
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diameter  if  the  glooal  fuel-air  ratio  were  stoichiometric.  The 
individual  nozzles  that  would  be  required  to  form  such  a  cloud 
must  therefore  be  capable  of  injection  to  71  m  (233  ft). 

The  small-scale  experimental  facility  consists  of  a  con¬ 
crete  test  pad  with  continuous  runways  connecting  inlaid 
instrumentation  modules.  Instrumentation  consists  of  stagnation 
and  static  piezo-electric  pressure  transducers  as  well  as 
high-speed  photography.  Test  sequencing  is  controlled 
electronically.  Two  impulsive  fuel  dispensers  have  been 
built.  In  one  of  these,  high-pressure  gas  (such  as  compressed 
air  or  nitrogen)  is  discnarged.  through  a  sonic  orifice  to  drive 
a  piston  which  forces  fuel  through  the  nozzle (s).  In  the 
other,  the  piston  is  driven  by  a  gas  generator. 

Observations  of  large-diameter,  impulsively  injected 
single  liquid  jets  have  indicated  that  several  breakup 
mechanisms  may  simultaneously  be  operative.  Breakup  at  the 
head  or  the  jet  appears  to  result  from  a  recurring  insta¬ 
bility.  This  might  be  a  Taylor  acceleration  instability. 

Breakup  along  the  sides  of  the  jet  also  occurs  ana  may  result 
from  shear  layer  stripping  or  from  a  Helmholtz  instability. 

Due  to  internal  turbulent  motion  the  jet  ceases  to  possess  a 
contiguous  core  above  some  height.  The  end  or  tail  of  the  jet 
breaks  up  in  what  appears  to  be  a  response  to  periodic  vortex 
shedding.  Material  removed  from  the  jet  as  droplets  or  stranas 
presumably  undergoes  successive  aerodynamic  shattering  ana 
subsequent  vaporization.  By  these  combined  mechanisms  ruel 
vapor  is  added  to  tne  air  along  the  entire  jet  trajectory. 

Tail  breakup  appears  to  limit  jet  reach  if  insufticient  fuel  is 
dispensed  through  the  nozzle. 

Small-scale  FAE  clouas  have  been  formed  by  injection  of 
fuel  through  a 'large  number  (600-1400)  of  radia lly-ai rectea 
nozzles.  Visually  uniform  spray  distributions  nave  been 
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achieved;  however,  interactions  between  the  expanding  spray 
cloud  and  the  air  can  lead  to  shape  distortions  if  the  injec¬ 
tion  velocity  is  too  high. 

The  small-scale  experiments  with  detonated  fuel-air 
clouds  are  believed  to  exhibit  reasonable  nuclear/FAE  airblast 
fidelity.  Experiments  in  this  area  are  continuing  with  the 
recently  upgraded  hardware  and  instrumentation.  The  experimen¬ 
tally-determined  effective  blastwave  energy  of  the  FAE,  with 
respect  to  a  scaled  1  KT  nuclear  explosion,  is  on  the  order  of 
40-75%  of  the  constant-pressure  heat  of  the  combustion 
depending  on  the  fuel  used  and  the  way  in  which  t.he  FAE  cloud 
is  formed.  A  theoretical  calculation  ot  the  airblast  emerging 
from  a  heptane  FAE  has  indicated  that  under  ideal  conditions 
the  higher  values  could  be  anticipated.  Experiments  have  shown 
that  the  FAE  airblast  characteristics  are  repeatable  and 
azimutha.lly  symmetrical.  The  quality  of  the  nuclear  airblast 
simulation  was  found  to  improve  with  range  in  both  the 
experimental  and  the  theoretical  work. 
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SECTION  I 


INTRODUCTION 

There  exists  an  ongoing  defense  requirement  to  assess  che 
response  of  various  objects  to  airblast  loaning  resulting  from 
nuclear  explosions.  Since  tne  1963  ban  on  above-ground  nuclear 
weapons  testing,  it  has  been  necessary  to  conduct  such  experi¬ 
ments  by  simulated  means.  For  far-field  airblast  simulation  of 
surface  bursts,  the  predominant  technique  has  been  the  use  of 
high  explosives.  Grouna-piane  hemispheres  and  surface- 
tangent  spheres  constructed  of  stacked  nigh-explosives  have  been 
detonated  for  this  purpose.  This  technique  has  the  advantage 
of  being  relatively  straightforward,  and  the  quality  of  air- 
blast  simulation  obtained  is  reasonably  satisfactory  at  over¬ 
pressure  levels  below  about  1  MPa  (150  psi) . 

Unfortunately ,  inexpensive  surplus  TNT  has  become  very 
scarce.  As  of  Marcn  1979  an  estimate  of  the  bulk  cost  of  TNT 
was  about  2.75  $/kg  (1.25  $/lom^ ,  Since  approximatley  4.54  x 
lO'3  kg  (10°  lom)  or  TNT  are  required  for  a  1  KT  nuclear- 
simulation,  a  program  of  frequent  simulated  testing  at  this 
energy-release  level  wouia  be  prohibitively  expensive. 

In  addition  to  cost,  high  explosive  (HE)  simulation  has 
otner  drawbacks.  A  substantial  crater  is  formed  at  the  explo¬ 
sion  center  during  a  test.  For  subsequent  testing  this  neces¬ 
sitates  either  relocation  of  all  instrumentation  to  a  new  test 
area,  or  considerable  earthmoving  to  refill  the  crater.  Either 
operation  is  time-consuming  and  costly.  Furthermore,  high 
explosives  generate  a  ground  shock  which  in  the  tar-field  is 
not  representative  of  a  nuclear  ground  snock.  The  ground  shock 
therefore  tends  to  interfere  with  experiments  intended  to  test 
solely  for  airbiast-loading .  Finally,  a  substantial  amount  of 
high  velocity  debris  and  ejecta  accompanies  HE  detonations. 
Impingement  of  these  materials  upon  test  objects  is  clearly 
undesirable . 


For  these  reasons,  concepts  tor  alternate  simulated  far- 
field  airblast  sources  have  been  explorea.  In  particular, 
attention  has  recently  focused  on  the  use  of  fuel-air  explo¬ 
sions  for  this  purpose  (References  1-9).  A  fuel-air  explosion 
results  from  fast  combustion  of  a  cloud  comprising  fuel  that  has 
been  dispersed  into  atmospheric  air.  The  event  includes  both 
the  combustion  process  within  the  cloud  itself  and  the  cumula¬ 
tive  subsequent  interactions  with  the  ambient  air. 

Two  modes  of  FAE  cloud  combustion  are  possible.  Defla¬ 
gration  alone  will  occur  if  the  cloud  is  not  too  large  and  if 
it  is  initiated  with  a  low  energy  blast  source  or  a  thermal 
ignition  source  (such  as  a  flame) .  Detonation  will  occur  if 
the  cloud  is  initiated  by  a  blastwave  or  other  shock  source  of 
sufficient  strength  and  duration.  The  explosions  that  result 
trotn  these  two  modes  possess  similar  characteristics  in  the 
extreme  far-field  (-7000  Pa  (1  psig))  but  they  are  very 
different  at  higher  overpressures.  The  highest  overpressure 
attained  by  a  deflagrative  explosion  is  a  function  of  the 
deflagration  velocity  and  is  in  general  significantly  less  than 
that  generated  by  a  detonative  explosion.  Since  it  is  desir¬ 
able  to  extend  the  blast  simulation  to  the  highest  pressure 
levels  possible,  only  detonative  explosions  are  of  present 
interest. 

The  characteristic  size  of  the  fuel-air  cloud  that  would 
be  required  to  simulate  a  1  KT  nuclear  airblast  is  exemplified 
by  estimates  assuming  propylene  oxide  as  the  fuel.  In  this 
case  at  stoichiometric  proportions  the  cloud  volume  would  be 
7.50  x  105  (2.65  x  10^  ft^)  and  would  contain  85,400 

kg  (188,000  lbm)  of  fuel.  For  a  hemispher ically  shaped  cloud, 
tne  diameter  would  be  142  m  (466  ft). 

Two  fundamental  issues  are  being  addressed  by  the  present 
research  program.  The  first  involves  the  determination  of  a 
satisfactory  means  by  which  a  cloud  having  the  required 


v*  u*  ■:*  ^  V-;  i 

:vSi 


^ 


"7  'V'." 


V  V  V  7  V  V  WW\' 


V  .' 
'  <1 


*•, ,  ♦; 

p  r 

t  i1 


dimensions  can  oe  formed.  The  second  involves  the  determina¬ 
tion  of  the  degree  ot  correspondence,  or  lidelity,  between  FAE 
and  far-fiela  nuciear  aicoiaets. 

For  tests  requiring  peak  overpressure  levels  above  about 
7  kPa  (11  psi),  it  nas  been  considered  preferable  to  torm 
hemispher ically-shapea  clouas.  A  nemispner ical  FAE  is 
tneor etically  the  most  efficient  shape  (Reference  10) .  Its 
symmetry  assures  the  .minimum  number  of  undesirable  secondary 
waves  (due  to  reflections) ,  and  allows  predictable  testing  in 
all  directions  at  the  highest  overpressure  levels  possible. 

Several  techniques  for  forming  hemispherical  FAE  ciouds 
have  been  suggested.  These  include  oalloon  containment, 
explosive  ruel  dissemination  rcom  an  array  or  canisters, 
rockets  propelled  oy  the  fuel  tnac  is  to  be  dispersed,  ana 
large-scale  hydraulic  fuel  injection. 


Of  tnese  techniques,  hydraulic  fuel  injection  appears  at 
present  to  be  the  most  viable.  Exquid  fuel  is  impulsively 
injected  into  the  atmosphere  at  high  velocity  through  a  large 
numoer  of  nozzles.  Each  liquid  jet  is  atomized  hydraulically, 
that  is,  without  induced  swirl  or  coaxial  air  injection.  The 
spray  columns  from  these  many  jets  overlap  to  fill  in  the 
cloud.  The  jets  can  be  injected  either  racially  from  a  central 
(point)  source  or  vertically  from  distriouted  ground  posi¬ 
tions.  The  former  seems  preferaole,  because  the  reach  required 
of  each  jet  would  be  tne  same  ana  oecause  the  amount  of 
plumbing  would  be  minimized  due  to  the  smaller  ground  area 
covered  by  the  dispensing  system. 

With  this  technique,  the  cost  of  materials  expended  per 
test  would  be  minimized,  bo  above-grouna  structure  is 
required,  and  potentially  harmful  debris  is  not  produced.  On 
the  other  hana,  such  a  tecnnique  has  not  been  previously 
developed.  Investigation  has  oeen  required  to  aetermme 
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whether  or  not  it  is  in  fact  possible  to  impulsively  generate 
jets  with  the  necessary  reach,  degree  of  atomisation,  and  fuel 
distribution.  The  technology  needed  to  impulsively  inject  a 
large  number  of  such  jets  simultaneously  has  also  requirea 
development.  The  tecnnique  has,  nowever,  been  shown  to  be 
feasible  in  preliminary  small-scale  experiments  (References 
2-5)  and  is  discussed  in  detail  in  this  report. 

Regardless  of  the  specific  technique  chosen  for  fuel  air 
cloud  formation,  the  quality  of  the  simulated  far-field  nuclear 
airblast  is  of  fundamental  concern.  In  particular,  the  cor¬ 
respondence  between  the  positive-phase  static  ano  dynamic 
pressure  profiles  as  functions  both  of.  radius  and  time  needs  to 
be  evaluated  in  terms  of  peak  pressures,  positive  phase  dura¬ 
tions,  impulse,  and  the  extent  of  secondary  shocks  and  other 
aberrations  in  the  pressure-time  profiles.  Blastwave  symmetry 
and  repeatability  are  also  of  importance.  The  results  of 
investigations  in  these  areas  are  discussed  in  tnis  report. 

The  report  is  structured  as  follows.  First,  experimental 
fuel  dispensing  systems  tnat  have  oeen  considered  and  those 
that  have  been  tested  at  small  scale  are  describe.'’  in  Section 
2.  In  Section  3  the  experimentally  ooserved  atomization 
cha racteristics  of  liquid  jets  are  discussed,  while  Section  4 
deals  with  the  formation  of  hemispherical  clouos  using  a  large 
number  of  such  gets.  In  Section  5  the  results  of  theoretical 
and  experimental  investigations  of  nuclear  airblast  simulation 
by  fuel-air  explosions  are  discussed.  A  concluding  discussion 
is  provided  in  Section  6. 
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SECTION  2 


FUEL  DISPENSING  SYSTEM 

Point-source  fuel  dissemination  requires  a  dispensing 
system  capable  of  impulsively  forming  a  large  number  of 
raaially-dir ectea  jets  simultaneously.  Tne  specific  require¬ 
ments  for  such  a  dispenser  are  closely  coupled  to  the  dynamics 
of  breakup  of  individual  liquid  jets.  (Jet  breakup  is  the 
subject  of  Section  3.)  However,  since  an  impulsive  dispenser 
is  required  to  investigate  jet  breakup,  development  of  the 
dispenser  has  been  iterative.  In  this  section  dispenser 
hardware  concepts  are  described  generally  and  experience  with 
two  specific  designs  is  discussed. 

Several  concepts  have  been  considered  and  scnematic 
sketcnes  of  tnem  are  shown  in  Figure  1.  The  first  of  these, 
whicn  appears  to  be  the  least  complex,  operates  on  a  principle 
similar  to  that  of  a  common  aerosol  dispenser,  A  gas  generator 
raises  the  pressure  above  a  liquid  supply  which  is  then  ejected 
through  a  dip  tube.  The  gas  generator  is  located  at  the  top  of 
the  dispenser  to  facilitate  reloading.  This  design  (Figure 
1(a))  had  been  previously  suggested  as  a  candidate  system 
(Reference  5) . 

There  are  three  problem  areas  with  this  design.  First, 
it  is  questionable  whether  the  liquid  in  the  dispenser  tank  can 
be  fully  ejected.  After  the  liquid  level  drops  bexow  the 
bottom  of  the  dip  tube  the  liquid  column  remaining  in  the  tube 
would  be  unstable.  Secondly,  even  if  this  last  fluid  were  to 
be  dispensed  smoothly,  subsequent  gas  ejection  through  the 
nozzles  would  disturb  the  spray.  This  would  be  detrimental  to 
the  formation  of  a  predictable  spray  pattern.  Finally,  the  hot 
produces  from  the  gas  generator  might  ignite  the  fuel  inside 
the  dispenser,  or  these  gases  could  ignite  tne  spray  cloud  as 
tney  issue  from  the  nozzle. 
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All  of  the  remaining  dispenser  concepts  in  Figure  1 
incorporate  an  interface  between  the  propelling  gas  and  tne 
fuel.  This  avoids  the  potential  problems  ot  the  aerosol 
dispenser  design.  The  first  of  these  concepts  is  illustrated 
in  Figure  1(0).  A  piston  ariven  oy  gas  that  is  admitted 
through  a  solenoid  valve  into  a  free  volume  behind  it  forces 
fuel  through  tne  nozzles(s).  The  valve  acts  as  a  choked 
(sonic)  orifice,  limiting  the  gas  flow  rate  from  a  high  pres¬ 
sure  supply  bottle.  In  this  design  the  valve  opening  time  also 
limits  the  rate  of  initial  piston  acceleration.  This  can  be 
problematic  if  the  acceleration  time  is  comparable  to  the  total 
fuel  dispensing  time.  Also,  rather  high  supply  bottle  pres¬ 
sures  are  requited  to  achieve  useful  working  pressures  behind 
tne  piston.  This  design  might  be  improved  by  incorporating  a 
fast-acting  doubie-oiaphragm  type  of  valve,  such  as  that  used 
in  snock  tubes.  A  sketch  of  this  alternate  arrangement  is 
shown  in  Figure  1(c). 

A  high  pressure  supply  bottle  can  be  avoided  through  the 
use  of  a  gas  generator  in  a  configuration  such  as  that  shown  in 
Figure  1(a).  A  propellant  grain  witn  variable  burn  rate  might 
be  used  so  as  to  tailor  the  dispensing  pressure  to  any. 
time-varying  profile  aesired. 

Another  concept  involving  comoustion  is  shown  in  Figure 
1(e).  In  this  case,  however,  the  free-volume  is  filled  with  a 
combustible  gas  mixture.  Dispensing  is  initiated  by  detonating 
this  mixture.  The  pressure  will  then  decrease  cue  to  gas  expan¬ 
sion  as  the  piston  displaces  the  fuel.  The  pressure  dropoff 
rate  is  governed  by  the  free  volume  size  and  by  gas  properties, 
but  it  is  not  otherwise  controllable. 

Ar  arrangement  that  does  not  involve  suaaen  pressuriza¬ 
tion  by  combustion  or  by  means  of  valving  is  shown  m  Figure 
1(f).  Here  the  tree  volume  pressure  is  raised  in  advance  to 
the  appropriate  initial  level  with  a  suitable  gas  (at  ambient 
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temperature) .  Toe  piston  is  restrained  by  a  cable  that  can  be 
severed  by  an  explosive  device  (such  as  an  exploding  bolt) .  As 
in  the  preceding  device  concept,  this  design  permits  only  a 
limited  amount  of  control  over  the  pressure  dropoff  rate, 
however,  nigh  temperature  gases  are  not  involved.  This  would 
simplify  construction  ana  eliminate  driver-gas  energy  losses 
aue  to  heat  transfer. 

All  of  the  schemes  shown  in  Figure  1  use  gas  pressure  to 
force  the  fuel  out  of  the  dispenser.  However  other,  more  novdl 
techniques  might  also  be  feasible.  As  examples,  tne  piston 
could  be  driven  by  a  rocket  motor  or  perhaps  even  by  a  large 
spring . 


Two  dispenser  concepts  (Figures  1(b)  and  1(d))  have  been 
examined  at  a  size  compatible  both  with  small-scale  fuel-air 
cloud  formation  experiments  ana  with  large  scale  single-ge.t 
reach  and  breakup  experiments. 


2,1  U-TUbE  DISPENSER 

A  rudimentary  dispenser  system  was  usea  exclusively 
during  the  early  pnases  of  the  present  research  program 
(Reference  5) .  It  is  essentially  an  adaptation  of  the 
choked-otif ice  arrangement  of  Figure  1(b).  A  sketch  of  the 
dispenser  is  shown  in  Figure  2.  Because  of  its  shape  it  has 
oeen  referred  to  as  the  U-tube  dispenser.  It  is  2.75  m  (9  ft) 
long  overall,  constructed  ut  15.2  cm  (6  in.)  schedule  40 
wrought  steel  pipe.  The  maximum  amount  of  liquid  that  can  be 

3 

dispensed  is  about  0.03b  m  (9.5  gal,).  The  centerline 
distance  between  tne  two  legs  is  46  cm  (18  in.).  The  piston 
shown  is  i5.2  cm  (6  in.)  long  ana  is  fitted  witn  two  teflon 
seals.  High  pressure  nitrogen  is  discharged  into  the  free 
volume  through  a  solenoid  valve  with  a  1.9  cm  (0.75  in.) 
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Figure  2.  Schematic  diagram  showing  the  U-tube  used  in  the 

experimental  investigation  for  disseminating  fuel  to 
form  hemispherical  fuel-air  explosive  clouds.  One  leg 
of  the  U-tube  was  pressurized  in  order  to  force 
the  fuel  through  a  nozzle  head  attached  to  the 
other  leg  of  the  U-tube. 
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effective  orifice  diameter.  'The  nitrogen  is  supplied  from  two 
manifolded  cylinders  having  a  total  volume  of  0.71  m^  (2.5 
ftJ)  . 

In  use  this  dispenser  was  buried  in  the  ground.  The 
volume  behind  the  piston  was  filled  with  water  which  efrec- 
tively  formed  a  flexible  U-shaped  extension  of  tne  solid 
piston.  The  back  side  of  the  extended  piston  was  then  acces¬ 
sible  .from  the  surface.  The  amount  of  water  behind  the  piston 
could  be  aujusted  so  as  to  permit  some  control  over  the  initial 
free  volume. 

The  pressure  level  requited  in  the  supply  bottles  is  much 
nigher  chan  that  in  the  free  volume  behind  the  piston.  An 
estimate  of  tne  necessary  supply  pressure  was  obtained  by 
assuming  quasi-steady-state  dispensing  (constant  piston 
velocity)  ano  ignoring  the  dropoff  in  supply  bottle  pressure. 

In  that  case  the  pressure  on  both  sides  of  the  piston  is  the 
same  ana  constant,  and  the  rate  of  increase  in  the  free  volume 
equals  the  rate  of  fuel  volume  efflux  through  the  nozzle(s). 
Referring  to  the  notation  in  Figure  2,  this  condition  can  be 
expressed  as 

vFV  -  V*  ,  (1) 


where 


V 


i 


(2) 


is  the  volume  flow  rate  of  the  liquid  fuel  and  m^ ,  ace  the 
fuel  mass  flow  rate  and  the  liquid  density  respectively.  The 
rate  of  increase  in  tree  volume,  V  ,  can  be  written  as 
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Here,  the  perfect  gas  law  has  been  used  to  express  the  instan¬ 
taneous  free  voLume  magnitude  in  terms  of  the  mass  of  gas 
contained  within  it,  m;  the  pressure,  Pp;  molecular  weight, 

W;  temperature,  T  of  the  gas;  and  R,  the  universal  gas 
constant.  Equating  Eqs.  (2)  ana  (3)  according  to  Ej.  (1) 
gives,  wich  Pp  and  T  constant, 


(4) 


where  m  =  dm/dt  represents  the  rate  at  which  gas  flows  from  the 
supply  bottle  into  the  free  volume. 

Now,  if  the  nozzle  area,  ana  the  U-tube  cross- 

sectional  area,  A  are  such  that 

u 

<w2  1 


then  the  liquid  mass  flow  through  N  such  nozzles  can  be  ex¬ 
pressed  as 


-  *V2°*(PI>  -  PA>1 


where  it  has  been  assumed  that  the  nozzle  discharge  coefficient 
is  equal  to  unity.  On  the  other  hand,  the  gas  mass  flow  rate 
through  the  orifice  (solenoid  valve)  is  given  by 
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where  A„  is  the  orifice  area  and 
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The  orifice  Mach  number  in  Eq.  (7)  is 
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and  Y  is  the  ratio  of  specific  heats  of  the  gas.  The  maximum 

value  of  M_  is  unity,  corresponding  to  a  critical  ratio 
£» 

P  /P  that  can  be  calculated  from  Eq.  (8).  This  is  the 

5  XT 

choking  point.  At  larger  ratios  of  P  ,/p  ,  achieved  by 

S  F 

reducing  Pp,  the  mass  flow  m  does  not  increase  further. 
Combining  Eq.  (4) -(6)  gives 


P_  (P_ 
P  P 


(9) 


This  equation  is  not  in  closed  form  since  f(MT?)  is  a  function 
of  Pg.  However,  at  the  choking  point,  Mg  =  1,  and  in  that 
case  Eq.  (9)  becomes 


For  a  fixed  choked  orifice  area  k„,  Equation  (10)  indicates 

hi 

that  the  required  supply  bottle  pressure  increases  linearly 
witn  the  liquid  nozzle  area  and  with  the  number  of  these 
nozzles.  As  an  example,  for  a  single  7.62  cm  (3  in.)  diameter 
nozzle  driven  at  pp  =  1.38  MPa  (200  psi) ,  ana  with  P^ 

»  1  x  105  Pa  (14.7  psi),  W  =  28,  Y  =  1.4,  and  P0  =  1000 

3  3  . 

kg/m  (62.5  lbm/ft  ),  the  supply  bottle  pressure  requirea 

in  conounction  with  a  1.91  cm  (0.75  in.)  diameter  choked 

orifice  at  room  temperature  woula  be  Ps  =  5.47  MPa  (658  psi 

Higher  free  volume  pressures  can  be  achieved  without 
excessive  supply  oottle  pressure  by  increasing  the  soienoia 
valve  orifice  area.  Larger  valves  generally  have  longer 
opening  times,  however,  and  since  fuel  dispensing  times  are 
characteristically  in  the  ran$e  100-500  ms,  valve  opening  times 
must  be  kept  very  short.  In  the  case  of  the  U-tube  system,  the 
solenoid  valve  opening  time  was  estimated  to  be  on  the  order  of 
50  ms.  Together  with  the  large  effective  mass  of  the  dispenser 
piston  (including  the  water)  this  resulted  in  rather  sluggish 
behavior  during  tne  initial  startup  perioa.  Transient 
oscillations  were  also  observed  on  occasion  during  startup. 

Free  volume  pressures  above  1.4  MPa  (200  psi)  could  not  be 
easily  obtained  when  testing  large  nozzles  (or  with  a  large 
number  of  small  nozzles)  because  of  the  unreasonably  high 
supply  pressures  that  would  be  required  according  to  Eq.  (10) . 

The  principal  advantage  of  a  choked-or if ice  driven 
dispenser  is  that  the  free  volume  pressure  does  not  decrease  as 
rapidly  as  does  the  pressure  in  the  supply  bottle  during  the 
dispensing  operation.  This  can  be  seen  in  Eq.  (10).  A  design 
of  this  kind  with  a  fast  valve  (for  example,  a  double 
diaphragm)  and  no  moving  water  mass  remains  open  for  con¬ 
sideration.  However,  it  is  felt  at  present  that  the  use  of  a 
solenoid-actuated  or  similar  valve  is  not  totally  satisfactory 
and  that  the  d-tube  configuration  is  not  optimum. 


2.2  LINEAR  DISPENSER 


t- 


A  dispenser  of  the  type  sketchea  in  Figure  1(d)  was  also 
designed  and  built  and  has  been  employed  in  more  recent 
experimental  work.  A  schematic  to  scale  of  this  device  appears 
in  Figure  3  and  a  photograph  is  shown  in  Figure  4.  This  was 
originally  referred  to  as  tne  linear  dispenser  to  distinguish 
it  from  the  earlier  U-tube.  The  piston  in  this  case  is  driven 
by  pressure  developed  from  a  gas  generator.  Details  of  the  gas 
generator  are  given  in  Appendix  B.  In  addition  to  testing  the 
gas  generator  concept,  this  dispenser  was  designed  to  improve 
performance  and  reliability  over  that  of  the  U-tube.  Water  can 
be  aooed  to  the  tank  below  the  piston  so  as  to  adjust  the 
initial  free  volume;  however  the  water  does  net  move  with  tne 
piston.  This  results  ir.  greatly  improved  inertial  start-up 
response.  The  piston  acceleration  time  was  reduced  further  by 
increasing  its  diameter. 

une  disadvantage  of  the  linear  design  is  that,  unlike  the 
U-tuoe,  the  linea.r  dispenser  must  be  removed  from  the  ground  in 
oraer  to  change  the  water  level  in  the  free  volume  tank  or  to 
reload  the  gas  generator  between  runs.  In  the  present 
installation  this  effectively  constrains  the  overall  dispenser 
length  to  about  3  m  (10  ft)  due  to  the  reach  limit  of  the 
equipment  on  hand  that  is  used  to  r~  ^e  the  dispenser.  For 
this  reason  the  linear  dispenser  capacity  was  limited  to  a 

3 

maximum  of  0.0273  m  (7.2  gal.)  of  liquid.  Tms  is  slightly 
greater  than  the  quantity  of  fuel  needed  for  the  hemispherical 
cloud  detonation  experiments,  but  it  is  somewhat  less  tha.n  the 
capacity  of  the  U-tube. 

For  design  purposes  one  start-up  inertial  response  of  the 
linear  dispenser  was  analyzed  to  first  order.  The  gas 
generator  was  assumed  to  develop  the  initial  pressure  PpQ  in 
the  initial  free  volume,  vfVq'  in  negligible  time.  The  mass  of 
fluid  in  the  nozzle  was  small  compared  with  that  of  the 
remaining  fluid  and  the  piston  mass,  and  was  therefore  ignored. 
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Figure  4,  Photograph  of  linear  fuel  dispenser a  A  10.2  cm 

(4  in.)  I.D.  x  102  cm  (40  in.)  long  nozzle  is  in¬ 
stalled.  The  dispenser  is  standing  on  the  center 
area  of  the  concrete  'test  pad.  In  use  it  is 
lowered  into  the  adjacent  hole. 


The  fluid  between  the  piston  and  nozzle  was  assumed  to  move  as 
block  flow  (uniform  velocity  throughout) .  The  motion  was 
further  assumed  to  be  quasi-steady  so  that  the  Bernoulli 
equation 

Pe  +  Hu3  -  PA  +  *  Ve  (11) 

and  steady-state  conservation  of  mass, 


could  be  used.  In  these  relationships,  points  "3"  and  ”E"  (see 
Figure  3)  designate  the  nozzle  inlet  and  exit  planes  respec¬ 
tively,  Ap  and  are  the  piston  and  r\ozzle  areas,  and 
is  the  ambient  pressure. 

The  equation  of  motion  tor  the  piston  plus  fluid  plug  is 

simply 

du 

(m£  +  mp)  *  (Pp  -  P^)  Ap  ,  (13) 


in  which  m^  and  mp  are  the  masses  of  the  liquid  and  the 
piston,  Pp  is  the  instantaneous  pressure  behind  the  piston, 
and  Up  is  the  piston  velocity.  This  equation  by  itself  shows 
that  the  pressure  drop  across  the  piston  plus  fluid  plug  goes 
to  zero  as  the  piston  attains  constant  velocity.  In  that 
condition  Pp  =  P^  and  the  pressure  drop  to  Pft  takes  place 
entirely  within  the  nozzle. 

Combining  Eqs.  (11)-(13),  with  u^  ~  u  ,  gives 


It  should  be  noted  in  this  equation  that  if  Ap  =  A^  (equal 
piston  and  nozzle  areas)  ,  the  piston  will  continue  to 
accelerate  indefinitely  as  long  as  Pp  >  and  no  steady 
velocity  will  be  reached  ir.  that  case.  It  is  also  clear  that 
increases  in  Ap  will  enhance  the  inertial  reponse  and  that  an 
increase  in  the  ratio  A^/A.,  will  reduce  the  time  to  reach  a 
near-steady  velocity. 

In  reality  the  free  volume  pressure  Pp  drops  off 
continuously  as  the  result  of  gas  expansion  during  dispensing. 
The  rate  of  free  volume  increase  is 


dV 


FV 


dt 


=  Vp 


(15) 


The  volume  increase  is  related  to  the  rate  of  fuel  expulsion  by 


dl% 

dt 


aV 


FV 


dt 


9 


and  assuming  isentropic  expansion,  the  pressure  changes  accord¬ 
ing  to 


Equations  (14) -(16)  were  solved  with  the  intial  condi¬ 
tions  Pp  =  PpQ,  Vpv  =  VFV0'  and  up  =  0  at  t  =  0.  This 
yielded  the  piston  velocity  as  a  function  of  time.  The  nozzle 
inlet  pressure  was  then  calculated  from 


=  PA  +  ^ *UP 


(17) 


32 


0  at  t 


obtained  by  combining  Eq.  (11)  and  (12) .  Since  up  = 

0,  it  is  clear  that  P3  is  initially  equal  to  Pft. 

It  was  originally  felt  that  the  U-tube  could  be  adapted 
for  gas-generator  driven  dispensing.  This  would  be  done  by 
eliminating  the  water  behind  the  piston  and  mounting  the  gas 
generator  at  tne  top  of  the  left-hand  leg  shown  in  Figure  2. 

If  the  U-tube  could  have  been  adapted  in  this  way  a  new 
dispenser  might  not  have  been  needed.  A  calculation  of  the 
type  just  described  was  carried  out  to  estimate  the  response  of 
the  U-tube  if  it  were  used  in  this  manner.  With  a  7.62  cm  (3 
in.)  nozzle,  the  pertinent  values  are  Ap/A^  =  4.0,  Ap 
=  0.0182  m2  (0.196  ft),  and  mp  =  7.51  kg  (16.5  lbm) . 
Characteristic  values  of  the  tollowing  parameters  were  used: 
mp  =  26.5  kg  (7  gal.  of  water),  Pn  =  1.83  topa  (265  psi), 

VFVq  =  0.0708  m  (2.5  ft  ) ,  Y  =  1.3,  PA  =  1  x  10D  Pa 
(14.7  psi) ,  and  pp  =  1000  kg/m3  (62.5  lbm/ft3) . 

Tne  result  is  plotted  in  Figure  5.  As  can  be  seen  the 
pressure  at  the  nozzle  inlet  continues  to  rise  curing  nearly  30 
percent  of  the  dispensing  period.  The  resulting  fluid  efflux 
would  tnerefore  not  be  impulsive  in  character.  The  u§e  of  the 
U-tube  modified  for  gas-generator  operation  was  tnerefore 
rejected.  In  order  to  obtain  a  more  impulsive  response  a 
larger  piston  area  was  clearly  required. 

The  linear  dispenser  was  accordingly  designed  with  a 
piston  diameter  of  29  cm  (11.5  in.),  nearly  twice  that  of  the 
U-tube.  An  even  larger  piston  would  have  been  desirable; 
however,  machining  costs  rise  very  rapidly  with  piston 
diameter,  and  the  availability  of  materials  was  limiting. 

With  this  larger  piston  the  start-up  dynamics  were 
improved.  Tne  actual  area  and  mass  of  the  piston  are  An 

ir 

=  0.0670  m2  (0.721  ft2)  and  mp  =  43.1  kg  (94.9  lbm).  In 
conjunction  with  a  7.62  cm  (3  in.)  nozzle,  the  p'.ston-to- 
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Figure  5. 


4i 

TIME,  t  (ms) 

Nozzla  inlet  pressure  as  a  function  of  time  during  opera¬ 
tion  of  U-tube  fuel  dispenser.  Values  used  in  calc  ula- 
tion  are  =  26.5  kg,  ^Pq 1.33  x  106  Pa,  V0  =  0,0708  m 

Y  =  1.3,  Pa  =  1  x  xo5  paf  Pe  =  1000  kg/m3,  Dp/DN  =  2 ,  Ap 
=  0.0182  m2,  mp  --  7.51  kg. 


nozzle  area  ratio  is  Ap/A^  =  14.7.  With  other  quantities 
the  same  as  in  the  previous  calculation  for  the  U-tube,  the 
inertial  response  that  was  expectea  is  shown  in  Figure  6.  The 
time  to  reach  maximum  nozzle  inlet  pressure  (or  maximum  nozzle 
exit  velocity)  is  in  this  case  onli  six  percent  of  the  uotal 
dispensing  time. 

The  piston  was  constructed  of  6Q61-T6  aluminum  alloy.  To 
minimize  its  weight  it  was  assembled  trom  27.3  cm  (10.75  in.) 

OD  x  0.93  cm  (0.365  in.)  wall  tubing  with  welded  end  caps.  The 
overall  length  of  the  piston  is  56  cm  (22  in.).  A  6.4  cm  (2.5 
in.)  long  nose  with  a  60°  taper  angle  was  included  to  proviae 
cushioning  during  piston  deceleration  (stopping).  The  maximum 
stroke  of  the  piston  is  41  cm  (16  in.).  At  the  upper  limit  of 
its  travel  it  is  stopped  by  an  annular  shoulder.  In  this 
position  the  nose  mates  with  the  nozzle  inlet  wall  with  a  0.5 
mm  (0.020  in.)  gap  between  the  taperea  surfaces. 

Two  piston  seals  are  employee  on  each  of  the  two  piston 
end  caps.  The  three  uppermost  seals  arc  4.3  mm  (3/16  in.) 
diameter  tef lon-encapsuiated  viton  o-rings  manutactured  by 
Chesterton,  Inc.  These  are  more  flexible  and  remain  resilient 
longer  than  solid  teflon.  They  can  enaure  sustained  exposure 
to  solvent  fuels  (in  particular,  propylene  oxide).  A  single 
teflon  spiral-type  backup  o-rmg  is  used  with  each  of  these 
three  seals.  The  lowermost  teal  is  a  6.4  mm  (0.25  in.)  wide 
aiuminum/bronze  step-cut  piston  ring.  This  ring  is  needed  to 
protect  the  other  seals  from  potentially  damaging  exposure  to 
the  hot  products  from  gas-generator  comoustion  and  it  also 
serves  as  a  cylinder  wiper. 

The  cylinder  in  which  the  piston  travels  was  fabricated 
from  32  cm  (12.75  in.)  OD  by  28  cm  (11.1  in.)  ID  A106  grade  B 
seamless  black  pipe.  This  was  bored  to  29  cm  (11.5  in.)  ID, 
honed,  flash  chrome  plated-  and  then  re-honed.  The  chrome 
plating  was  applied  to  provide  corrosion  protection.  The 
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TIME,  t  (ms) 

Figure  6,  Nozzle  inlet  pressure  as  a  function  of  time  during  opera- 

tion  of  linear  fuel  dispenser.  Values  used  in  ca.lcula-  H 

tion  are  mj,  =  26.5  kcr,  -  1.83  x  10^  Pa,  V0  =  0.0708  m3,  pj 

y  =  1.3,  PA  =  1  x  10$  Pa,  p*  =  1000  kg/m3,  DP/Dn  =  3.83,  ^ 

Ap  =  0.0670  m2,  mp  =  4  3.1  kg. 
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piston/cyiinder  assembly  was  attached  to  the  other  aispenser 
parts  using  the  techniques  of  hydrualic  cylinder  construction. 
Eight  1-9  cm  (0.75  in.)  tie  bolts  connect  the  cylinder  to  tue 
free  volume  chamoer .  Tne  free  volume  is  about  0.093  m3  (3.3 

3 

ft  )  with  the  gas  generator  installed  ana  without  any  water 
fill. 

In  use  tne  linear  aispenser  is  installed  oelow  grade  in  a 
steel  cylinder  that  .is  51  cm  (20  in.)  ID  oy  2.6  m  (103  In.) 
long  and  has  a  2.5  cm  (1  m.)  thick  welded  steel  case  plate. 
This  cylindrical  housing  was  constructed  from  the  gas  bot_le  of 
a  surplus  torpedo.  A  fork  lift  truck  is  used  when  it  is 
necessary  to  remove  the  dispenser  from  the  housing.  A  2.7  m 
(9  1 1 )  A-frame  gantry  with  hoist  is  used  when  raising  the 
dispenser  just  far  enough  to  service  the  gas  generator. 

Experience  with  the  linear  aispenser  has  been  for  the 
most  part  positive,  with  the  exception  of  difficulties 
associated  witn  tne  gas  generator  (Appenaix  B) .  An  additional 
problem  was  experienced  with  the  method  used  to  vary  the 
internal  volume  of  the  free  volume  tank.  It  was  found  that 
aaaing  water  to  the  free  volume  tank  would  lead  to  a 
substantial  loss  of  aispensing  pressure.  Tnis  was  apparently 
due  to  energy  absorption  by  the  water  which  was  evidently 
penetrated  by  the  exhaust  jet.  The  effect  was  aggravated  with 
greater  amounts  of  water  present.  No  suitable  alternate  fill 
to  replace  water  was  round.  Any  suitable  alternate  fill  must 
be  unaffected  by  exposure  to  high  temperature,  violently 
turbulent  gases.  It  should  generate  no  dust  or  grit  which 
could  damage  the  piston  or  cylinaer,  ana  it  should  not  be 
propelled  by  tne  high  velocity  get  into  damaging  projectiles. 
Heavy  steel  chain  was  tried  as  a  fill  but  this,  like  the  water, 
apparently  absorbea  much  of  the  gas  energy. 


37 


SECTION  3 


TURBULENT  LIQUID  JETS 

The  projection  and  atomization  characteristics  of  im¬ 
pulsive,  turbulent  liquid  jets  are  of  importance  to  the  point- 
source  FAE  cIouq  rormation  technique  described  in  the 
Introduction.  The  oreakup  of  steady,  laminar  liquid  jets  has 
been  studied  extensively  ana  a  substantial  body  of  literature 
is  available.  Turbulent  liquid  jets  on  the  other  hand  have 
received  comparatively  little  attention.  In  any  case  the 
primary  breakup  characteristic  investigated  by  other 
researchers  nas  been  the  point  of  jet  rupture.  This  is  the 
point  along  the  jet  trajectory  beyond  which  there  is  no 
continuous  path  that  connects  fluid  particles  to  the  nozzle. 
Very  rew  experiments  are  reported  in  which  ultimate  jet 
penetration  distances  (reacii)  have  been  measured. 

Impulsive  hydraulic  fuel  injection  has  also  received 
little  previous  attention  other  chan  in  the  areas  of 
high-pressure  jet  cutting  and  diesel  fuel  injection.  In  both 
of  these  applications  very  small  diameter  jets  are  involvea. 
Except  for  earlier  work  at  Systems,  Science  and  Sortware 
(Reference  5)  and  that  reported  herein,  it  appears  that  no 
prior  systematic  experiments  have  been  conducted  with  large 
diameter  (to  10  cm  (4  in.)),  impulsively  injected,  turbulent 
liquid  jets  m  the  intermediate  pressure  range  (  2.8  MPa  (400 
psi)). 

Impulsive  (single-pulse  or  finite  column  length)  jets  are 
distinguished  from  steady  state  (continuous  or  infinite  column 
length)  jets  on  the  basis  of  their  duration.  If  the  time  for 
an  element  of  jet  fluid  to  traverse  the  distance  from  tne 
nozzle  exit  to  the  point  of  ultimate  jet  reach  is  short 
compared  to  the  total  dispensing  time,  the  jet  can  be  con¬ 
sidered  steady  state.  if  these  two  times  are  comparable,  tne 
jet  is  considered  impulsive. 
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The  dispensing  time  is  approximately 


t 


D 


UDAN 


(18) 


where  Vp  is  the  total  volume  of  fluid  that  is  dispensed  at 

constant  velocity,  u^,  through  a  nozzle  of  area  A^.  The 

time,  tD,  for  a  fluid  element  co  attain  full  reach,  R  ,  is 
K  C 


Therefore  the  condition  for  an  impulsive  jet  is 


or,  in  terms  of  volume  dispensed 


V  ^ 

F 

MAX 


Vc 


(20) 


This  is  just  the  volume  tnat  would  be  contained  in  a  virtual 
liquid  column  of  length  F  and  cross-sectional  area  AN.  If 
the  actual  fuel  volume  dispensed  is  very  much  larger  than  this 
value  the  virtual  liquid  column  is  said  to  be  of  infinite 
length  and  the  jet  is  essentially  steady. 


3.1  STEADY -STATE  JETS 

An  early  classification  of  jet  breakup  characteristics 
was  given  by  Ohnsorge  (Reference  11) .  This  is  reproduced  in 
Figure  7.  The  classification  was  developed  in  terms  of  the 
liquid  Reynolds  number. 
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where  is  the  nozzle  diameter  and  y^  and  a ^  are  the  liquid 
viscosity  and  surface  tension.  As  the  velocity  of  the  jet  is 
increased  at  fixed  diameter,  the  dominating  breakup  mechanism 
changes  from  symmetrical  varicose  surface  wave  growth 
(Reference  12)  to  helicoidal  wave  oevelopement  (Reference  13) . 
At  still  higher  velocities,  Ohnsorge  reported  that  the  jets 
were  fuiiy  disrupted  or  atomizeu  at  the  nozzle  exit. 

For  reference  the  regime  of  interest  to  the  present 
investiation  is  designated  by  the  shaded  area  in  Figure  7.  The 
regime  spans  the  range  0.35  MPa  1  AP  -  2.8  MPa,  1.6  -  1 

100  mm  (50  -  AP  =  400  psi,  1/16  -  DN  -  4  in.).  It  is  noted 
that  nozzle  pressure  drop,  AP,  and  jet  velocity,  uQ  are 
related  by 


AP  -  hp^2  (21) 

assuming  loss-ire^  flow  in  the  nozzle.  As  can  be  seen  on  the 
figure,  ail  of  the  jets  of  interest  lie  in  the  complete 
atomization  regime  according  to  the  Ohnsorge  classification. 

Complete  jet  disintegration  at  the  nozzle  exit  does  not 
however  appear  to  be  an  accurate  descript  ->n  of  turbulent  jet 
breakup  in  many  cases,  kore  recent  investigations  clearly  snow 
that,  although  surface  atomization  begins  in  high-velocity 
turbulent  jets  at  tne  nozzle  exit,  a  self-contiguous  central 
jet  core  persists  for  a  consideraole  distance.  In  terms  of  the 
liguid  Weoer  number 
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the  distance  to  this  rupture  point,  xfi,  has  been  correlated 
empirically  by  Grant  and  Middleman  (Reference  14)  as 


(22) 


for  100  1.  Wg  L  105,  and  by  Phinney  (Reference  15)  as 
x  u 

~  =  55  +  1.085  W J  (23) 

UN  e  A 


<  <  4 

for  the  approximate  range  35  ~  W'e  ~  6  x  10  .  These  two 

& 

relations  ao  not  agree  well  with  one  another.  They  ao  show 
nowever  that  for  example  at  W-  ^  10 ^  the  distance  to  the 
rupture  point  can  be  as  much  as  400  diameters  from  the  nozzle 
exit.  It  should  be  noted  that  in  turbulent  jets,  this  rupture 
does  not  appear  to  occur  abruptly.  Rather,  the  jet  core  is 
dissociated  gradually  ana  the  rupture  point  simply  designates 
the  furthest  position  o.f  self -continuity  of  the  core. 

Many  investigators  have  concluded  that  turbulent  liquid 
jets  are  ultimately  shattered  by  the  unconstrained  radial 
motion  of  internal  eddies  (References  15-20)  .  The  smallest  of 
these  perturb  the  surface  but  appear  to  be  contained  by  it. 
Those  eadies  with  sufficient  radial  momentum  penetrate  the 
surface  and  are  swept  away  by  aerodynamic  forces.  The  amount 
of  material  stripped  from  tne  jet  by  this  mechanism  is 
relatively  small.  However,  the  resulting  spray  partially 
conceals  the  inner  jet  core,  wnich  has  an  irregular,  mottled 
surface.  It  is  the  large-scale,  low-momentum  eddies  wnich 
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eventually  produce  'surface  disturbances  of  sufficient  magnitude 
as  to  induce  utlimate  varicose  breakup  of  the  jet  core. 

The  intensity  and  character  of  turbulence  upstream  of  the 
nozzle  exit  therefore  play  key  roles  in  determining  both  the 
extent  of  early  spray  removal  and  the  ultimate  reach  of  the 
jet.  The  turbulent  intensity  is  increased  oy  surface  rough¬ 
ness,  by  increasing  the  nozzle  lengtn,  and  by  non-smooth 
plumbing  transitions*.  Other  characteristics  of  specific 
nozzle  geometry  seem  to  be  of  little  consequence  (References  16 
and  21) .  Differences  in  the  scale  and  intensity  of  turbulence 
at  the  nozzle  exit  could  perhaps  account  for  tne  differing 
descriptions  giver,  of  turbulent  jet  breakup  by  the  various 
investigators . 

The  distance  to  the  rupture  point  of  a  turbulent  liquid 
jet  and  its  ultimate  penetration  distance  (reach)  are  not 
directly  related.  Rupture  is  associated  witn  tne  formation  of 
liquid  si  rands  or  ligaments  (Reference  22)  whicn  subsequently 

collap'.'  nnder  ^he  action  of  surface  tension  into  droplets. 

* 

Lara^i,  liquid  c  .obules  with  diameters  on  the  order  of  the 
intia*  ;et  diameter  may  also  be  formed  (References  23  and 
24) .  The  droplets  and  globules  are  then  broken  up  into  a  finer 
spray  by  aerodynamic  shattering  if  their  velocity  is  high 
enough.  In  fact,  the  process  of  aerodynamic  shattering  can  be 
expected  to  be  repeated  with  the  formation  of  successively 
smaller  droplets  unti1  the  spray  momentum  has  been  largely 
depleted.  Ti.^  surv:vi».y  droplets  then  presumably  undergo  pnase 
change  oy  convection-assisted  evaporation.  In  jets  with  a 
vertical  velocity  component,  these  may  or  may  not  pass  through 
a  trajectory  extremum  i  fall  to  the  ground  before  being  fully 
vapor ized . 


*It  should  be  noted  that  the  nozzles  used  in  the  development  of 
Eqs.  (22)  and  (23)  were  very  long  (L/D  'v  100),  so  that  the 
turbulence  was  presumably  well  developed. 
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Aerodynamic  shattering  of  liquid  droplets  is  a  complex 

phenomenon.  A  good  summary  discussion  is  given  in  Reference 

25.  Shattering  occurs  when  the  gas-liquid  Weber  number  for  the 

2 

droplet,  w  =  p.u  L/\i  is  greater  than  about  ten.  The 

6  A  36 

mode  of  shattering  changes  as  Wg  is  increased.  At  low  Weber 
numbers  (WQ  5  100)  the  parent  droplet  is  broken  into 
relatively  large  secondary  drops,  while  for  -  1000  the 
parent  droplet  is  atomized  into  a  fine  microspray,  typically 
1/100  of  the  parent  droplet  size.  Breakup  times  and  droplet 
trajectories  have  been  correlated  with  coarse  accuracy  by  a 
number  of  researchers  (Reference  25).  This  information  could 
conceivably  by  used  at  some  future  time  to  predict  jet 
penetration  beyond  the  rupture  point;  however,  neither  the 
droplet  size  distribution  nor  the  droplet  velocity  at  the 
rupture  point  is  presently  known  with  sufficient  accuracy  for 
t'hat  purpose. 

For  these  reasons  the  determination  of  ultimate  jet  reach 
is  best  obtained  experimentally.  Data  of  this  type  are  very 
limited,  even  for  steady  state  jets.  Prior  interest  in  jet 
reach  has  been  confined  largely  to  studies  of  fire  stream 
projection  (Reference  26) .  In  most  cases  the  maximum  distance 
that  the  jet  remains  suitable  for  firefighting  has  been 
measured,  rather  than  the  ultimate  jet  reach.  However,  Box 
(Reference  27)  and  Freeman  (Reference  21)  give  ultimate  jet 
reaches  for  vertical  fire  streams  up  to  5  cm  (2  in.)  diameter 
and  at  pressures  to  0.5  MPa  (70  psi) .  A  correlation  in  this 
range  given  by  Freeman,  rewritten  here  in  dimensionless  format, 
is 

£ 

~  =  A  -  1.12  x  lo'4  X2  ,  (24) 

N 
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where 


AP 


u„ 


p*gDN 


2gD 


N 


Here,  AP  is  the  pressure  drop  across  the  nozzle,  =  1000 
kg/m3  (62.5  lbm/ft3)  ,  g  =  9.80  ni/s2  (32.2  ft/sec2)  is 
the  acceleration  due  to  gravity,  and  is  the  nozzle 
diameter.  It  is  interesting  to  note  that  ultimate 
dimensionless  reach  as  given  by  Eq.  (24)  decreases  witn 
increasing  nozzle  diameter  (at  fixed  velocity) ,  while  the 
dimensionless  rupture  distance,  given  by  Eqs.  (22)  and  (23) 
increases.  However,  the  ranges  of  validity  of  these  three 
empirical  relations  do  not  coincide. 

The  accuracy  of  Eq,  (24)  rapidly  diminishes  for  nozzle 
sizes  larger  than  about  3.8  cm  (1.5  in.),  and  for  pressures 
above  about  0.5  Mpa  (70  psi) .  Increasing  the  diameter  or 
pressure  beyond  these  limits  leads  to  much  less  increase  in 
actual  reach  than  is  predicted  by  Eq.  (24).  At  sufficiently 
high  pressures  a  decrease  in  the  reach  of  these  water  jets  has 
been  observed.  The  physical  constraints  that  limit  the  upper 
bound  on  the  reach  that  could  be  achieved  by  an  arbitrarily 
large  jet  are  not  presently  known.  However  the  greatest 
vertical  jet  reach  that  has  been  documented  is  170  m  (560  ft) 
(Reference  28) .  This  is  a  steady  state  water  jet  that  is 
developed  from  an  annular  nozzle  at  2.8  MPa  (400  psi). 

A  limited  number  of  survey  experiments  were  performed  at 
Systems,  Science  and  Software  with  small  diameter,  quasi-steady 
liquid  jets  at  pressures  to  2  MPa  (300  psi).  These  small 
diameter  jets  are  relevant  to  the  formation  of  small-scale 
experimental  FAE  clouds.  The  individual  jets  were  produced 

3 

using  the  U-tube  dispenser.  A  total  of  0.011  m  (3  gal.)  of 
liquid  was  discharged  through  five  nozzles.  The  nozzles  were 
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drilled  into  a  blind  ena  flange  as  shown  in  Figure  8.  All  of 
the  nozzle  bores  were  13  ram  (0.50  in.)  in  length  so  that  the 
nozzles  naa  xeng cn- to-diameter  ratios  varying  between  3.5  and 

3.5.  The  U-tube  was  pressurized  by  nitrogen  gas  from  a  9.8 

-43.3  , 

x  10  m  (60  m  )  supply  bottle.  During  dispensing  the 

pressure  in  this  bottle  dropped  from  an  initial  value  near  2 

MPa  (300  psi)  to  atmospheric  pressure.  Tnis  took  place  in 

aoout  10  s.  The  reach  of  each  jet  was  measured  as  a  function 

of  time  on  films  taken  during  a  test  and  correlated  with  the 

measured  nozzle  entry  pressure. 

Results  from  tests  with  JP4  kerosene,  NAPTHA,  and  pro¬ 
pylene  oxide  (P.O.)  were  qualitatively  similar.  Characteristic 
oata  from  these  tests  is  given  in  Figures  9  and  10.  below 
about  0.34  MPa  (50  psi)  the  get  reach  w.as  found  to  increase 
rapidly  with  pressure.  However  the  lacx  of  atomization  in  this 
range  made  the  jets  nearly  invisible  on  the  films  and  reach 
measurements  could  not  be  taken.  Tnis  difficulty  accounts  for 
the  considerable  data  scatter  near  0.34  MPa  (50  psi).  It 
snouid  also  be  noted  tnat  in  all  cases,  accurate  definition  of 
the  instantaneous  jet  reach  was  very  difficult  and  could 
account  for  some  of  the  oata  scatter.  Above  0.34  MPa  (50  psi) 
get  atomization  was  substantial  and  increased  subjectively  with 
increasing  pressure.  The  jets  became  visible  on  tne  films 
allowing  reach  to  be  measured.  The  reach  remained  roughly 
constant  and  in  some  cases  decreased  somewnat  in  this  range. 
This  implies  balancing  Competition  between  increasing  jet 
momentum  ano  the  increased  drag  per  unit  mass  tnat  results  from 
the  formation  of  smaller  spray  droplets  at  the  higher  injection 
velocities.  In  essence  the  utlimate  .jet  reach  was  found  to  be 
a  function  only  of  nozzle  diameter  above  some  minimum  jet 
velocity. 

It  should  be  noted  that  during  the  early  experiments  of 
tnis  kind,  only  cursory  attention  was  given  to  prevailing 
atmospneric  winds.  These  were  estimated  at  tne  time  to  be  no 
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more  than  about  1.3  m/s  (3  MPH) .  some  of  the  cests  were  later 
repeated  during  completely  quiescent  weather.  These  later  tests 
were  also  of  greater  duration.  About  0.036  m3  (9.5  gal.)  of 
fuel  was  dispensed  and  the  pressurized  supply  bottle  volv  ae  was 
increased  to  0.033  m3  (1.17  ft3).  These  changes  resulted  in 
a  dispensing  time  of  about  30  seconds.  The  pressure  range 
covered  was  0.7-2  MPa  (100-300  psi) .  The  jet  reaches  measured 
in  these  later  tests  were  substantially  greater  than  those  in 
the  earlier  experiments  (by  as  much  as  a  factor  of  two) ,  but  it 
is  believed  that  the  increases  were  virtual^  It  is  felt  that 
in  a  very  still  atmosphere  the  extremely  fine  spray  droplets 
hovered  in  a  cloud  near  the  jet  peak.  These  fine  droplets  were 
then  displaced  vertically  by  air  currents  induced  by  the  jet. 

As  a  result,  the  apparent  jet  reach  was  increased.  With  any 
slight  breeze  or  with  much  shorter  dispensing  period,  this 
phenomenon  could  not  occur.  It  is  believed  for  this  reason 
tnat  the  earlier  tests  are  the  more  credible. 

Some  of  the  data  scatter  in  Figures  9  and  10  is  no  doubt 
due  to  momentary  wina  gusts  that  occurred  during  the  experi¬ 
ments.  In  consequence  the  trends  indicated  by  the  curves  are 
probably  accurate  only  within,  say,  about  ±10%.  The  general 
sensitivity  of  ultimate  reach  to  wind  has  been  noticed  by  other 
investigators  even  in  iarge  diameter  jets  at  lower  pressures 
(with  comensurately  larger  spray  droplets)  (Reference  10) . 

3.2  IMPULSIVE  JETS 

Impulsively  injected  jets  are  of  greater  interest  to  the 
present  investigation  than  are  steady-state  jets.  A  hemi¬ 
spherical  fuel-air  cloud  having  a  radius  on  the  order  of  70  m 
(230  ft)  is  required  to  simulate  the  far-field  airblast  of  a 
1  KT  ncuiear  explosion.  To  form  this  clouo  a  large  number  of 
radially-directed  jets  are  to  be  discharged  from  a  central 
point  source.  The  cloud  must  be  formed  in  less  than  about  two 


seconds  to  avoid  cloua  distortion  ana  arift  due  to  atmospheric 
winds.  Therefore  steady-state  jets  are  unsuitable. 

An  impulsive  jet  from  a  single  nozzle  would  ideally 
atomize  fully  into  a  uniform,  fine  spray  and  would  attain  the 
required  reach.  Unfortunately,  techniques  to  enhance  jet 
atomization  tend  to  diminish  jet  reach.  For  example,  it  is 
well  known  that  induced  swirl  and  impinging  or  coaxial  air  jets 
promote  early  liquid  jet  atomization  and  the  formation  of  small 
droplets.  These  benefits  however  accrue  at  the  expense  of 
reach.  Maximum  jet  reach  is  instead  obtained  by  sustaining  jet 
coherence  over  the  greatest  possible  trajectory  length  and  by 
encouraging  the  initial  formation  of  large  droplets  upon  jet 
rupture.  For  this  latter  purpose  a  simple,  converging 
hydraulic  nozzle  that  is  terminated  at  the  orifice  nas  been 
found  to  be  tne  most  effective  (Reference  16).  Additives  such 
as  long-chain,  high  molecular  weight  polymers  can  additionally 
be  used  to  dramatically  increase  jet  conerence  (References  19 
and  2y) .  However,  for  the  application  to  spray  cioua 
formation,  excessive  jet  coherence  is  also  no.t  desirable. 

Highly  coherent  jets  would  not  atomize  sufficiently  to  fill  in 
the  voids  between  them  in  clusters  of  jets  and  the  extremely 
large  droplets  would  be  difficult  to  detonate.  Therefore,  a 
hydraulic  nozzle  designed  to  generate  moderate  but  not  minimum 
turbulent  intensity  is  believed  to  be  the  best  practical 
tradeoff.  Tnis  will  produce  a  reasonably  coherent  jet  that 
however  will  ultimately  atomize  into  a  detonable  spray. 

Impulsive  jet  breakup  has  been  investigated  at  Systems, 
Science  and  Software  in  part  with  the  U-tube  facility  and  in 
part  with  the  linear  dispenser.  The  linear  dispenser  driven  by 
the  smokeless-powaer  gas  generator  as  described  in  Appendix  B 
is  cur r ently' limited  to  1.4  MPa  (200  psi) .  In  order  no  conduct 
experiments  at  higher  pressures  an  alternate  technique  was 
used.  In  tests  above  1.4  MPa  (200  psi) ,  a  brass  diaphragm  was 
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placea  over  the  nozzle  inlet  ana  the  dispenser  freevolume  was 
pressurized  with  nitrogen  to  within  about  0.17  MPa  (25  psi)  of 
tne  diaphragm  rupture  point.  To  start  a  test,  nitrogen  gas 
from  a  large-volume  supply  bottle  at  5.5  MPa  (800  psi)  was 
allowed  to  flow  rapidly  into  the  freevolume  tank.  This  quickly 
raised  the  pressure  in  the  freevolume  tank  causing  the 
diaphragm  to  rupture.  This  procedure  was  considered  a 
contingency  measure  since  there  was  some  question  as  to  the 
influence  of  the  diaphragm  on  the  jet  breakup. 

Impulsively  dispensed  jets  from  nozzles  having  diameters 
between  2.5-10  cm  (1-4  in.)  have  been  examined.  In  all  cases 
the  nozzle  design  is  very  simple,  co...  isting  of  a  60°  conical 
inlet  tnat  converges  to  the  required  exit  diameter.  This  is 
followed  by  a  straight  section  that  is  ten  nozzle  diameters  in 
length.  Tne  junction  between  these  two  sections  is  rounaed 
with  a  2.5  cm  (1  in.)  radius  of  curvature.  The  straight 
section  was  adaea  to  allow  a  liraitea  amount  of  turbulence  to 
develop. 


The  characteristics  of  impulsive  jet  breakup  are 
illustrated  stylistically  in  Figure  11,  while  Figures  12 
through  15  are  outline  sketches  tracea  from  high  speed  films 
taken  during  the  experiments.  It  should  be  noted  that  the 
samples  chosen  for  this  group  of  figures  involve  different 
liquids  and  nozzle  sizes. 

Many  breaxup  features  of  impulsive  jets  are  similar  to 
those  of  steady-state  jets.  The  mechanism  of  jet  surface 
penetration  by  turbulent  eddies  is  believed  to  be  primarily 
responsible  for  jet  rupture.  This  is  a  process  that  for  a 
given  moving  fluid  element  is  essentially  independent  of  events 
ahead  of  or  behind  that  element.  It  should  therefore  be  little 
affected  by  the  jet  duration.  Similarly,  those  mechanisms 
leading  to  jet  surface  erosion  as  well  as  those  associated  with 
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Figure  12.  Head  and  surface  breakup.  Prop 
1.91  cm  (0.75  in)  diameter  nozz 
(155  psi)  driving  pressure  drop 
spray  near  the  core  is  shown.) 
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at  x.38  MPa  (200  psi)  driving  pressure  drop.  (Only  the 
spray  near  the  core  is  shown.) 


VWATnW  VWNfV* l >  .  •  A  -  *  W  V.i>  . 
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Figure  14.  Early  jet  spreading.  Water  jet  from  1.91  cm  (0.75  in) 
diameter  nozzle;  nozzle  pressure  1.03  MPa  (150  psi) . 


the  behavior  of  droplets,  ligaments,  and  globules  formed  after 
rupture  appear  to  be  quite  similar  for  both  steady-state  and 
impulsive  jets. 

On  the  other  hand,  impulsive  jets  uifler  signf icantiy 
from  steady-state  jets  because  of  the  existence  of  both  a 
leaning  front  (head  or  nose)  and  a  terminating  tail.  These 
influence  the  quantity  and  distribution  of  residual  spray  along 
the  jet  path. 

As  an  impulsive  jet  initially  emerges  from  a  nozzle,  its 
head  appears  to  undergo  recurring  instabilities  which  result  in 
tne  periodic  shedding  of  spray  sheets.  Qualitatively  these 
have  the  appearance  of  Taylor  instabilities;  that  is,  accelera¬ 
tion  induced  waves  which  grow  within  the  liquid  at  the  jet  head 
until  they  are  removed  by  aerodynamic  drag.  It  appears  that 
the  removed  liquid  forms  a  spray  which  rapidly  decelerates. 

The  suddenly  exposed  liquid  core  at  the  jet  head  passes  through 
this  spray  and  the  instability  cycle  repeats.  This  repeating 
process  contributes  to  the  Christmas-tree  character  of  the 
early-time  spray  immediately  adjacent  to  the  jet  core,  which  is 
apparent  in  Figure  12,  In  other  investigations  similar 
oehavior  has  been  observed  with  impulsively  started  fuel  jets 
of  small  diameter  (200  ym  (0.008  in.))  (Reference  22). 

The  tail  of  an  impulsive  jet  appears  to  have  a  much 
greater  effect  on  jet  breakup  than  that  at  either  the  head  or 
external  surfaces.  An  example  of  the  appearance  of  an 
impulsive  jet  in  the  vicinity  of  its  tail  is  given  in  Figure 
13.  Only  the  dense  region  near  the  core  (not  the  entire  spray 
outline)  has  been  traced  on  this  figure.  The  average  velocity 
of  the  instantaneous  tail  of  the  contiguous  column  exceeds  that 
of  the  head  by  a  factor  of  1.5  to  2.  This  indicates  that 
material  breaks  off  from  the  end  of  the  jet  core  in  a  continu¬ 
ing  process.  The  tail  breakup  could  for  example  be  the  result 
of  unbalanced  lateral  forces  associated  with  periodic  vortex 
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shedding.  In  any  case,  material  removed  from  the  tail  is  very 
rapidly  dispersed  as  a  spray.  The  passage  of  the  tail  is 
always  accompanied  by  a  considerable  widening  of  the  residual 
spray  column  below  the  rupture  point. 

Tail  breakup  is  very  beneficial  in  the  present  context 
because  it  is  primarily  responsible  for  the  majority  of 
residual  spray  in  the  first  several  hunared  jet  diameters  from 
the  nozzle  exit.  The  only  other  contributions  to  the  spray  in 
this  region  are  from  jet  surface  atomization  and  from  remnants 
of  spray  shed  during  passage  of  the  head.  On  the  other  hand, 
the  tail  breakup  process  can  reduce  jet  reach.  As  the  total 
quantity  of  liquid  dispensed  impulsively  through  the  nozzle  is 
reduced,  a  condition  may  be  obtained  in  which  the  advancing  jet 
tail  overtakes  the  head,  if  this  occurs  before  the  head 
reaches  the  normal  rupture  point,  the  jet  will  be  broken  up  at 
a  distance  closer  to  the  nozzle  than  is  usual.  This  would  in 
turn  result  in  diminished  total  jet  travel  or  reach.  The 
minimum  quantity  of  liquid  that  can  be  dispensed  through  a 
given  nozzle  without  reach  reduction  by  this  mechanism  cor¬ 
responds  to  the  condition  in  which  the  advancing  tail  breakup 
point  overtakes  the  head  at  exactly  the  normal  jet  rupture 
point. 


An  estimate  of  this  minimum  quantity  of  liquid  that  must 
be  dispensed  can  be  obtained  from  a  simple  analysis.  Neglect¬ 
ing  jet  surface  and  head  atomization  and  assuming  that  the  jet 
moves  at  a  constant  velocity,  uD,  its  length  at  the  moment  of 
tail  emergence  from  the  nozzle  is 
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Assuming  that  the  tail  bteakoff  point  propagates  along  the  jet 
at  twice  the  jet  velocity  relative  to  the  ground,  the  rate  of 
change  in  the  contiguous  jet  length  L  would  be 
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so  that 


A  plot  of  this  equation  for  water  and  heptane  at  an  injection 
velocity  of  =  60  m/s  (200  ft/s)  is  given  in  Figure  i6. 

Tests  with  heptane  discharged  vertically  through  a  2.5  cm 
(1  in.)  diameter  nozzle  at  uQ  =  60  m/s  (200  ft/s)  were 
carried  out  to  determine  if  this  anticipated  phenomenon  is  in 
fact  operative.  The  minimum  volume  estimated  by  Eq.  (26)  under 
tnese  conditions  is  Vpmin  =  0.012  m3  (3.2  gal).  This  is 
about  half  the  capacity  of  the  linear  dispenser  so  that  testing 
both  at  dispensed  volumes  above  and  below  the  predicted  minimum 
was  possible.  Results  from  these  experiments  are  plotted  in 
Figure  17.  Reach  measurements  were  made  from  films  taken 
during  the  tests.  These  were  checked  against  measurements 
taken  from  a  simple  nanc-helo  sight  stick.  Reduced  reach  was 
experienced  when  less  than  0.0076. m^  (2  gal)  of  heptane  was 
dispensed.  These  results  are  certainly  in  the  predicted 
direction  but  it  is  felt  tnat  more  tests  will  be  required  to 
conclusively  verify  the  tail  breakup  effect. 

before  the  tail  of  an  impulsive  jet  emerges,  tne  shape  or 
the  jet  spray  in  the  immediate  vicinity  of  the  nozzle  is 
conical.  Tnis  is  due  to  surface  atomization  and  subsequent  air 
entrainment  or  droplets  by  the  two-phase  (droplets/air)  mixture 
in  tnat  region.  The  remaining  spray  column  shape  is  roughly 
cylindrical  except  at  the  jet  head.  After  the  tail  emerges  the 
spray  column  also  becomes  cylindrical  near  the  nozzle.  This 
can  be  seen  in  Figures-  14  and  15. 

The  average  spray  column  width  and  height  continue  to 
increase  for  a  considerable  time  after  dispensing  has  ended. 
This  is  evident  in  Figures  18-20.  (The  dispensing  times  for 
all  tests  plotted  in  these  figures  were  less  than  0.6  s.) 

Figure  18  shows  the  increase  in  average  jet  width  as  a  funtion 
of  time  and  Figures  19  and  20  are  plots  of  the  head 
trajectories  of  characteristic  vertical  jets.  In  both  of  the 
latter  two  figures,  a  calculated  trajectory  assuming  simple 
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Figure  16.  Estimated  minimum  dispensed  liquid  volume 
versus  nozzle  diameter. 
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FUEL  QUANTITY  DISPENSED  (10  ‘'iti'3) 

Figure  17.  Ultimate  jet  reach  vs.  fuel  quantity  dispensed 
(2.54  cm  (1  in.)  diameter  nozzle,  heptane, 
nominal  driving  pressure  =  1.23  MPa  (180  psi)) 
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Figure  20 .  Comparison  of  stream  height  as  a  function  of  time 
for  two  experiments  involving  the  projection  of 
propylene  oxide  from  a  3.81  cm  (1.5  in.)  diameter 
nozzle. 
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lateral  aerodynamic  surface  drag  in  conjunction  with 
gravitational  deceleration  has  been  plotted  for  reference. 

Below  some  breakaway  point  this  curve  can  be  maae  to  agree 
quite  well  with  the  measured  trajectories  using  a  drag  coef¬ 
ficient  of  CD  =  0.25.  The  influence  of  gravity  is  small 
compared  to  the  aerodynamic  drag.  The  rapid  reduction  in  jet 
velocity  beyond  the  breakaway  point  is  presumably  associated 
with  aerodynamic  drag  on  spray  droplets  following  jet  rupture. 
The  breakaway  points  in  these  experiments  aid  not  however 
correlate  with  the  rupture  distances  given  by  either  Eq.  (22) 
or  Eq .  ( 23) . 

Ultimately  the  jet  attains  a  trajectory  maximum  after 
which  the  spray  droplets  vaporize  or  else  fail  partially  or 
fully  to  the  ground,  depending  on  the  liquid  volatility  and 
droplet  size.  Jet  widths  in  excess  of  250  nozzle  diameters 
have  been  observed  at  these  late  times.  However  the  jets  are 
nearly  if  not  fully  atomized  at  much  earlier  times.  Since  it 
is  of  importance  to  minimize  spray  formation  times,  the 
earliest  moment  at  which  the  jet  has  developed  into  a  detonable 
cylinaricai  cloud  is  of  interest.  On  the  films  of  impulsive 
jets  it  appears  that  this  moment  occurs  snortly  after  the  tail 
reaches  the  vicinity  of  the  rupture  point  or  the  breakaway 
point.  On  that  basis  a  tentative  definition  for  the  formation 
time  of  an  impulsive  jet,  measured  from  the  beginning  of  the 
dispensing  period,  is 
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where  tD  is  the  dispensing  time.  Maying  use  of  Eq.  (18)  this 
can  be  written 
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This  definition  is  somewhat  arbitrary.  The  second  term 
represents  an  upper  limit  estimate  of  the  time  for  the  con¬ 
tiguous  get  tail  to  reach  the  rupture  point.  In  writing  that 
term  it  was  assumed  that  normal  get  rupture  occurs  at  distances 
less  than  1000  nozzle  diameters  and  that  the  jet  tail  moves  at 
twice  the  velocity  of  the  head.  Since  ultimate  jet  reach  is 
typically  1200-2000  nozzle  diameters  it  represents  a  reasonable 
upper  estimate. 

Some  measurements  of  impulsive  jet  reach,  taken  approxi¬ 
mately  at  the  formation  time,  are  given  in  Figure  21  as  a 
function  of  jet  velocity.  In  this  figure  the  fuel  is  heptane, 

3 

the  nozzle  size  is  1.3  cm  (0.5  in.),  ana  0.0064  m  (1.7  gal) 
of  fuel  were  dispensed.  The  reach  behavior  of  impulsive  jets 
appears  to  be  similar  to  that  of  steady  jets  when  sufficient 
liquid  is  aispensed.  Beyond  a  certain  minimum  jet  velocity  the 
reach  appears  to  be  essentially  constant.  The  knee  of  the 
reach  curve  (the  point  at  which  the  curve  begins  to  level  off) 
also  corresponas  roughly  to  the  minimum  jet  velocity  for  which 
atomization  can  be  visually  observea. 

In  future  work,  it  is  recommenaed  that  jet  reach  ana 
width  at  the  formation  time  be  measurea  systematically  as 
functions  of  jet  velocity  and  diameter.  The  minimum  volumes  of 
fuel  that  can  be  dispensed  in  order  to  achieve  these  reaches, 
and  the  minimum  jet  velocity  for  sufficient  atomization  should 
be  determined. 
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SECTION  4 

HEMISPHERICAL  FUEL-AIR  CLOUD  FORMATION 
4.1  CLOUD  SIZE  REQUIREMENTS 

The  volume  of  the  fuel-air  cloud  and  the  quantity  of  fuel 
required  are  fixed  by  the  nuclear  yield  to  be  simulated.  The 
relationships  are  not  straightforward  but  are  asociated  with 
the  energy  release  available  through  combustion  of  the  fuel-air 
cloud  and  with  the  way  in  which  this  energy  is  disposed  hyaro- 
aynamically . 

The  total  energy  release  resulting  from  the  detonation  of 
a  fuel-lean  or  stoichiometric  fuel-air  cloud  is 


E  =  m„H^  (28) 

o01  ci  F  c 

1 ,  * 

wnere  m.,  is  the  mass  of  fuei  in  the  ciouu  and  H  is  tne 
r  c 

low- tempera ture ,  constant  pressure  heat  of  combustion  per  unit 
mass  of  fuel.  This  expression  assumes  that  the  considerably 
dissociated  detonation  products  achieve  low  temperature 
equilibrium  during  the  expansion  that  follows  airblast  develop¬ 
ment  from  the  fully  reacted  cloua. 


Due  to  the  finite  spacial  ana  temporal  distribution  of 
the  energy  released,  the  airblast  produced  by  a  fuel-air 
explosion  differs  from  a  point  source  olastwave  of  tne  same 
energy,  E0t  However,  at  fixed  range  from  an  FAE,  the 
variables  of  interest  (such  as  overpressure  and  overpressure 
impulse)  may  maiviaually  or  collectively  exhibit  poinc-source 
(or  nuclear  source)  characteristics.  If  so,  it  is  possible  to 
aefine  an  effective  energy,  E0e  p  for  the  FAE  airblast  at 
chat  range.  This  is  the  energy  of  a  point  or  nuclear  source 
which  at  the  same  range  would  yield  the  same  values  for  the 
’ariable  or  variables  of  interest. 


'.V  W\ 


In  general/  the  value  of  Eo^  ,,  so  determinea  will  vary 

ij  f  r 

with  range  and  with  the  particular  blast  variable (s)  being 
considered,  as  well  as  with  the  reference  airblast  itself 
(point  source  or  nuclear) .  Over  limited  range  intervals,  the 
variation  in  Eoc,  „  with  range  may,  however,  be  small  enough 

b  t  r 

that  a  single  average  value  can  be  used  with  acceptable 
accuracy.  The  extent  to  which  a  single  value  or  E0]e  F  is 
successful  in  specifying  FAE/nuclear  airblast  equivalence  nas 
been  examined  in  the  present  investigation.  The  procedure 
followed  to  establish  a  value  for  Eq^,  f  is  discussed  in 
Section  5.  It  has  also  been  found  useful  for  discussion 
purposes  to  normalize  the  effective  blastwave  energy  by 
defining  an  effective  airblast  efficiency;  i.e., 


n 


F 


(29; 


It  should  be  emphasized,  however,  that  the  value  of  nF  does 
not  physically  represent  the  amount  of  energy  that  is 
partitioned  to  hydrodynamic  blast. 


Substituting  Eq.  (28)  into  ^ 2 9 )  and  solving  for  m„  gives 

r 


m  =  F-,F 

F  “  Vc 


(30) 


This  is  the  quantity  of  fuel  required  in  order  to  yield  a 
specified  effective  blast  source  energy  on  the  presumption 
that  n,,  is  known.  The  quantity  of  air  that  is  required  for 

r 

combustion  with  this  amount  of  fuel  is  simply 

“a  “  V*A  '  <31) 

where  <j>A  is  the  fuel-air  ratio  by  mass.  (Stoichiometric 
proportions  are  designated  by  <£*.; 
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The  size  of  the  cloua  that  contains  the  required  amount 
of  air  depends  on  whether  the  fuel  is  dispersed  as  a  vapor  or 
as  a  spray.  If  dispersed  as  a  condensed  spray,  the  volume 
occupied  by  the  spray,  V  ,  is  normally  negligible  compared 
with  the  volume  occupied  by  the  air,  VA-  I.e.,  from  Eq.  (31), 


PAV 


-  3  -4 

This  ratio  is  typically  of  the  order  10  to  10  .  The 

size  of  the  fuel-air  cloud  in  this  case  is  essentially  the 
volume  of  air  containing  mass  m  .  That  is, 
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rn  kT  . 

V  =  v  =  — — — 
A  P  h 
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Using  Eq.  (3(J;  and  (31)  and  solving  for  the  radius  of  a 
Hemisphere  containing  this  volume  leads  to 
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If  on  the  other  hand,  the  fuel  is  dispersed  as  a  vapor, 
or  as  a  spray  which  then  fully  evaporates,  the  volume  of  the 
mixture  is 


where  tne  mixture  mass  and  molecular  weight  are  related  to 
those  of  the  fuel  and  air  by 
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These  relations,  along  with  Eqs.  (30)  and  (31)  give  a  cloua 
radius  o£ 
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(33) 


The  difference  in  cloud  size  calculated  in  these  two  ways 

12 

is  small.  For  example,  for  Eoc  „  =  2.09  x  10  J  (the 

w  t 

effective  far  field  blast  source  energy  of  a  1  KT  nuclear 
explosion),  assuming  a  stoichiometric  mixture  of  propylene 
oxide  and  air,  and  using  nF  s  0.744*,  a  spray  cloud  radius  of 
RCg  =  69.2  m  ana  a  vaporized  cloud  radius  of  RCv  =  70.4  m 
are  predicted.  ^Pnysical  properties  of  selected  fuels  are 
given  in  Table  1.) 


The  temperature  and  pressure  used  in  tne  preceding 
example  were  ordinary  atmospheric  values.  Actually,  it  the 
cloud  is  formed  as  a  spray  ana  then  evaporates,  the  cloua 
temperature  will  be  reduced  somewhat.  To  estimate  the 
magnitude  of  this  effect,  tne  evaporation  can  be  treated  as  a 
simple  heat  removal  process  from  the  ait,  for  which  at  constant 
pressure  the  energy  equation  is 
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Vi 
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Q  =  h2  -  HX 


*Tnis  value  was  determined  experimentally  in  small-scale  FAE 
tests  using  propylene  oxide  and  is  discussed  in  Section  5. 
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The  heat  effectively  added  is  Q  =  -m^L  where  L  is  the 
latent  heat  of  vaporization  per  unit  mass  of  fuel.  Treating 
the  air  as  calorically  perfect,  the  energy  equation  then  gives 


-mFL 


•n  C  AT 
a  PA 


or 


AT  =  -<J>AL/Cp^ 


(34) 


This  effect  can  account  for  an  absolute  temperature  drop  on  the 
order  of  10%  which  may  be  sufficient  to  be  considered  in  cloua 
formation  schemes  involving  fuel  spray  evaporation. 


4.2  FUEL  DISTRIBUTION  UNIFORMITY 

At  trie  present  time  it  is  felt  preferable  to  form  spray 
clouas  that  will  not  largely  vaporize  prior  to  detonation 
initiation.  Tnis  implies  low  vapor-pressure  fuels  which  are 
advantageous  in  that  boiling -induced  jet  breakup  (Reference  i0) 
is  minimized  or  eliminated.  The  more  essential  basis  for  this 
judgement,  however,  derives  from  fundamental  differences 
between  gas-  ana  two-phase  detonations. 

Detonations  which  occur  in  premixed  gas-phase  reactants 
possess  a  hydrodynamic  reaction  zone  that  is  typically  no  more 
tnan  a  few  millimeters  thick.  The  reaction  zone  is  charac¬ 
terized  by  regularly  recurring  localized  ignition  centers. 

Fast  reaction  at  the  ignition  centers  produces  blastwaves  which 
interact  with  the  detonation  front  and  with  adjacent,  similarly 
produced  blastwaves.  Their  collisions  in  the  directions  trans¬ 
verse  to  tne  wavefront  produce  new  ignition  centers.  Combus¬ 
tion  is  also  augmented  behind  these  transverse  shocks.  The 
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propagation  velocity  ana  detonation  limits  of  gas-phase  detona¬ 
tion  depend  on  the  stoichiometry  (at  given  pressure  and 
temperature) .  This  is  because  che  fuel  and  oxidizer  are  mixed 
on  a  molecular  level  in  advance  of  che  detonation  so  that 
reaction  mechanisms  involving  all  of  the  species  initially 
present  can  occur  at  once  in  the  gas-phase  and  are  limited  only 
by  the  chemical  kinetics.  A  fuel-air  explosion  to  be  generated 
from  a  cloud  in  which  the  fuel  is  initially  fully  vaporized 
should  therefore  possess  a  near -stoichiometr ic  global 
equivalence  ratio  to  optimize  the  energy  release  per  unit  mass 
of  mixture.  Moreover,  uniform  mixing  throughout  the  FAE  cloud 
would  represent  an  engineering  objective  for  the  fuel 
uisseminator  in  this  case.  Either  fuel-rich  or  fuel-lean 
pockets  in  the  cloud  would  reduce  the  total  cloud  energy 
release . 

Detonations  in  two-phase  mixtures  in  which  the  fuel  is 
dispersed  as  a  fog  with  very  small  (v2  pm)  droplets  exhibit 
most  of  che  characteristics  of  gas-phase  detonations  (Reference 
31).  However,  when  th*  droplets  are  larger  the  time  for  their 
conversion  to  vapor  is  increased.  The  reaction  zone  of  a 
detonation  in  such  a  spray  is  protracted  and  its  character 
becomes  significantly  altered  (Reference  32-36) .  Mixing  of 
fuel  witn  the  oxidizer  occurs  only  after  it  has  been  removed 
from  the  aroplecs  by  a  process  of  mass  stripping.  This  results 
while  the  droplets  are  accelerating  in  the  convective  oxidizer 
gas  flow  behind  the  deconation  shock  front.  Stripping  appears 
to  involve  liquid  ooundary  layer  removal,  Helmholtz  surface 
waves,  and  shattering  from  a  Taylor  acceleration  instability. 

With  some  fuel/oxidizer  combinations,  fuel  removed  from 
tne  droplets  in  the  early  stages  of  stripping  accumulates  in 
the  recirculation  zones  of  the  individual  droplet  wakes.  After 
a  aelay  period  localized  explosions  occur  in  these  wakes. 

Other  fuel/oxidizer  combinations  do  not  exhibit  explosive 
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ignitioi.s.  In  axl  cases,  however,  the  energy  release  rate  in 
the  reaction  zone  is  limited  by  the  mass  stripping  processes 
rather1  than  by  chemical  kinetics.  The  long  reaction  zones 
which  result  are  responsible  for  a  generally  one-dimensional 
structure  that  is  very  different  from  that  of  gas-phase 
detonations . 

because  the  mass-stripping  process  in  two-phase  detona¬ 
tions  is  rate  limiting,  the  initial  fuel  (droplet)  distribution 
uniformity  is  not  as  critical  as  it  is  in  gas-phase  detona¬ 
tions.  The  fuel  that  is  stripped  from  a  given  droplet  mixes  in 
the  intensely  turbulent  reaction  zone  with  gases  originating  an 
appreciable  distance  away.  The  droplets  can  thus  be  con¬ 
ceptually  viewed  as  moving  sources  of  fuel  vapor  within  the 
well-stirred  environment  of  the  reaction  zone. 

The  propagation  velocity  of  a  fuel-lean  two-phase 
detonation  depends  as  would  be  expected  on  the  global  fuel/air 
stoichiometry.  However,  evidence  exists  in  tne  literature 
indicating  tnat  condensed-phase  fuel  in  excess  of  tne  amount 
required  for  a  globallly  stoichiometric  mixture  does  not 
participate  in  the  detonation.  In  one  study  (Reference  37) 
identical  detonation  velocity  ana  structure  results  were 
obtained  with  stoichiometric  and  with  very  fuel-rich  two-phase 
(spray)  mixtires.  In  a  similar  investigation  (Reference  38) 
the  fuel-rich  limit  could  not  be  determined.  Additional 
studies  (References  31,  39,  and  40)  have  demonstrated 
aetonability  of  fuel  films  at  global  equivalence  ratios  as  high 
as  600.  All  these  results  suggest  that  the  conoensed  fuel 
vaporizes,  mixes  with  the  oxidizer  gas,  and  reacts  as  quickly 
as  it  is  removed  from  the  parent  droplet  (or  liquid  layer,  in 
the  case  of  films) .  Except  in  the  final  stages  of  oxygen 
consumption,  combustion  of  the  fuel  vapor  therefore 
continuously  occurs  in  a  globally  fuel-lean  environment, 
regardless  of  the  quantity  of  fuel  that  is  latently  available 
in  the  conaensed-phase  at  any  given  moment.  Combustion 
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essentially  ceases  when  ail  of  the  oxygen  in  the  oxidizer  has 
been  cc  sumed.  It  appears  that  fuel  vapor  added  to  the 
products  by  further  droplet  stripping  beyond  tnat  point  does 
not  affect  the  detonation,  even  though  the  equilibrium  com¬ 
position  of  the  products  woula  be  sniftea.  Ala.  these  processes 
are  absent  from  a  gas-phase  detonation  in  which  the  fuel  and 
oxidizer  are  premixed  on  a  molecular  level  so  that  all  of  the 
fuel  initially  present  affects  the  instantaneous  combustion 
chemistry,  as  men  .*  i  earlier. 

The  requirements  for  uniform  fuel  distribution  and  local 
stoichiometry  in  the  formation  of  a  two-phase  fuel-air 
cloud  are  therefore  not  exacting,  provided  that  some  moderate 
sacrifices  in  fuel  economy  can  be  tolerated.  Any  quantity  of 
liquid  fuel  can  evidently  be  dispersed  into  a  given  volume  of 
air  provided  that  i 1 )  the  fuel  in  condensed  form  does  not 
occupy  a  significant  fraction  of  the  total  spacial  volume 
either  locally  or  globally,  and  (2)  the  droplets  are  not  spaced 
so  far  from  one  another  in  any  local  region  as  to  preclude 
involvement  of  ail  of  the  adjoining  air  in  the  comoustion 
processes*.  When  these  two  conditions  are  met,  the  total 
energy  release  is  expected  to  be  proportional  to  the  mass  of 
a_i_r  enveloped  by  the  cloud,  and  to  be  relatively  independent  of 
the  amount  of  fuel,  provided  an  excess  above  stoicniometr ic 
exists  everywhere.  The  behavior  expected  of  the  detonation  is 
as  though  the  clouo  were  uniformly  stoichiometric.  Only  the 
expansion  wave  properties  would  be  modified. 


*In  this  context  is  it  noted  that  typical  droplet  spacings  are 
on  the  order  of  ten  droplet  diameters  in  documented  spray 
detonations . 
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4.3  CLOUD  FORMATION  FROM  MULTIPLE  JETS 

A  current  objective  has  been  to  generate  nemispher ical 
clouds  comprising  unevaporated  fuel  sprays  that  are  distributed 
so  as  to  be  at  least  sightly  fuel-rich  throughout.  The  actual 
fuel  distribution  however  is  produced  in  a  very  complex  way 
which  at  present  is  not  fully  understood.  The  fuel  distribu¬ 
tion  is  affected  by  the  breakup  and  spreading  processes  of 
individual  jets  from  which  the  hemispherical  cloud  is  to  be 
formed,  and  by  the  interactions  between  adjacent  jets  including 
the  manner  in  which  they  overlap. 

The  effect  of  adjacent  jet  overlap  on  the  fuel  distri¬ 
bution  in  a  hemispherical  cloud  can  be  assessed  to  first  order 
by  assuming  that  the  spray  cloud  from  each  individual  jet  is 
cylindrical  and  that  adjacent  spray  cylinders  diffuse  into  one 
another  without  interaction  and  geometric  distortion.  These 
assumptions  may  be  reasonable  tor  dilute  sprays  as  is  the  case 
in  FAE  clouds.  It  two  adjacent  spray  cylinders  of  radius  Rj 
originate  from  the  same  point  (the  nozzle  cluster),  and  if  the 
angle  between  their  centerlines  is  0 ,  then  it  is  straight¬ 
forward  to  show  that  they  will  overlap  to  a  distance 


from  the  vertex.  If  R_,  <  R_,  (the  cloud  raaius),  then  at 

£  C. 

radii  between  R_,  and  R  the  spray  cylinders  will  not 
overlap.  In  this  case  the  cloud  will  consist  of  a  core  region 
of  overlapped  spray  cylinders  and  an  outer  region  characterized 
by  non-overlapping  finger-like  projections  (the  appearance  has 
oeen  compared  to  that  of  a  sea  urchin)  .  It  is  clearly  desir¬ 
able  to  design  the  dispensing  system  so  that  K_  =  r 

.c  C 

Since  typically  k^/D^  ^  50  and  R^/D^  n,  1500  at  the 
formation  time-  setting  =  A-,  in  Eq.  (35)  implies  0  ^ 

3.8°. 
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Due  to  jet  spray  overlap,  a  uniform  spray  distribution  in 
the  hemispherical  cloud  would  not  result  from  uniform  distribu¬ 
tion  in  the  individual  spray  cylinders.  The  total  mass  of  fuel 
in  one  isolated  spray  cylinder  is  nip/N  where  N  is  the  total 
number  of  jets  used  to  form  the  cloud.  If  the  special  density 
of  fuel  in  this  cylinder  is  p then  clearly 


/R,, 
C 

o 


p ' A-dR 
F  J 


(36) 


In  this  equation,  AT  is  the  cross-sectional  area  of  the  spray 
cylinder  (net)  and  pp  is  considered  to  be  a  function  of  R. 

Now  the  contribution  by  this  single  jet  to  the  total  mass  of 
fuel  m  a  hemispherical  shell  of  thickness  uR  at  radius 
R  <  Rc  from  the  nozzle  cluster  is 


ci  mi,  =  P'A.,.dR 

E  1  J 

Therefore  the  total  mass  in  the  shell  contributed  by  all  N 
overlapping  jets,  assuming  that  R„  =  R„,  is 

r  C 


am„  =  Nam'  =  NpAAt^R 

E  r  E  J 


The  volume  of  this  differential  shell  is  dV  -  2ttR  dR.  Hence 
the  fuel  spray  density  at  radius  R  resulting  from  tne  over¬ 
lapping  jets  has  an  average  value  of 


_  dm^ 
PF  dv~ 


that  is, 
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It  is  emphasized  that  chis  is  an  average  value.  Since  tne 
degree  of  overlap  varies  with  radius  the  distribution  would 
actually  vary  additionally  with  both  azimuth  and  elevation. 

The  implication  from  Eqs.  (36)  and  (37)  is  that,  as  a 

result  of  jet  spray  overlap,  uniform  average  cloud  density 

obtains  only  when  the  spray  density  of  individual  jet  cylinders 
?  2 

increases  as  R'.  That  is,  if  oc  R  ,  those  two  equations 
lead  to 


Any  otner  variation  in  will  result  in  a  non-constant 

Pp.  This  is  somewhat  fortuitous.  Although  attempts  have  not 

been  made  to  measure  p'  experimentally,  on  the  films  it 

r 

appears  subjectively  that  the  spray  in  fully-formed,  'impulsive 
single  jets  is  in  fact  distributed  non-unitormly .  The  spray 
density  near  the  nozzle  is  low  and  it  increases  gradually  in 
what  appears  to  be  an  accelerating  fashion  with  respect  to 
increasing  radius.  This  has  provided  some  encouragement  that 
hemispner ical  clouds  With  adequate  spray  distribution 
uniformity  can  in-  fact  be  formed. 


4.4  EXPERIMENTAL  INVESTIGATIONS 

Experiments  to  investigate  hemispner ical  cloud  formation 
by  means  of  multiple,  radially  injected  impulsive  jets  have 
been  conducted  on  a  small  scale  (R  =  4.6  m  (15  ft)).  Both 
tne  U-tube  and  linear  dispensers  have  been  used  at  different 
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times.  Nozzle  clusters  have  consisted  of  a  number  of 
countersunk  raaial  holes  drilled  into  machined  hemispherical 
sceel  or  spun  aluminum  domes.  One  of  these  domes  in  shown  m 
Figure  22.  This  particular  dome  nozzle  was  designed  in 
accordance  with  Eq.  (35);  i.e.,  with  Rp  =  Rc.  other  domes 
have  been  tested  having  600  to  1400  holes  (nozzles).  Hole  sizes 
between  1.6-3. 2  mm  (1/16-1/8  in.)  and  dispensing  pressures  from 
0.34-1.4  MPa  (50-200  psig)  have  been  tested.  In  most  of  the 

3 

experiments  the  volume  of  liquid  dispensed  was  0.0265  m  (7 
gal).  Liquids  tested  have  included  water,  propylene  oxide  and 
heptane. 

Examples  of  tests  illustrating  the  features  of  these 
clouds  are  shown  in  Figures  23  and  24.  The  four  clouds  that 
appear  in  Figure  23  were  generated  from  domes  having  600  3 . 2  mm 
(1/8  in.)  diameter  holes.  The  photogrpahs  are  frames  from 
high-speed  films.  The  angle  between  jets  was  about  6°,  whicn 
is  consider aoly  greater  than  the  value  of  3.8°  discussed 
earlier  as  being  required  for  complete  jet  spray  overlap 
throughout  the  cloua.  The  resulting  shape  (which  has  been 
likened  to  that  of  a  sea  urchin)  is  apparent  in  Figures  23(b) 
and  (d)  . 

Since  is  advantageous  to  minimize  the  number  of  nozzles 
(so  as  to  reduce  the  complexity  and  cost  of  the  dispenser),  it 
was  at  first  felt  that  the  problem  of  inadequate  overlap  could 
be  overcome  by  dispensing  at  higher  pressures.  That  would 
produce  individual  jet  sprays  of  greater  diameter,  thereby 
increasing  Rp .  Presumably  R  =  Rc  could  in  that  way  be 
reached  without  changing  9  (see  Eq.  (35)).  Hence  the  number  of 
nozzles  could  be  held  to  600.  Dispensing  at  higher  pressures 
also  would  increase  the  jet  velocity,  reduce  the  dispensing 
time  and  the  cloua  formation  time,  and  produce  a  more  finely 
atomized  spray,  with  little  if  any  effect  on  jet  reach. 


Figure  23.  Examples  of  clouds  formed  with  600-nozzle  dispenser  head.  (a)  Propylene  oxide 
at  1.24  MPa  (180  (psi) ,  (b)  P.O.  at  0.414  MPa  (60  psi) ,  (c)  heptane  at  1.24  MP 

(180  psi),  (d)  heptane  at  0.965  MPa  (140  psi).  Dimensions  given  on  the  photo¬ 
graphs  are  approximate. 


Unfortunately,  it  was  found  that  increases  in  jet 
velocity  lead  to  cloua  shape  distortions.  These  distortions 
are  the  result  of  air  motions  induced  oy  momentum  transfer  from 
the  expanding  spray.  The  character  of  the  distortion  changed 
from  relatively  insignificant  edge  lifting  to  severe,  large- 
scale  dishing  as  the  jet  velocity  was  increased.  The  magnitude 
of  these  effects  can  oe  seen  by  comparing  Figures  23(a)  and  (b) 
or  Figures  23(c)  and  (d) . 

Based  on  these  results  it  was  concluded  that  jet  overlap 
would  have  to  be  accomplished  by  increasing  the  number  of 
nozzles  rather  than  by  increasing  the  jet  velocity.  The  number 
of  nozzles  was  accordingly  increased  from  600  to  1400,  which 
gives  0  'v.  4°.  Clouds  produced  with  the  1400-nozzle  head  were 
then  very  uniform  in  visual  appearance  ano  did  not  have  the 
spiked  sea  urchin  character.  However,  tney  were  now 
surprisingly  dished,  even  at  pressures  as  low  as  0.3  MPa  (45 
psi)  which  corresponds  to  a  heptane  jet  velocity  of  30  m/s  (100 
ft/s).  This  is  illustrated  by  the  example  in  Figure  24. 

The  origin  of  these  clouo-shape  distortions  has  been 
examined  by  observing  the  trajectories  of  spray  droplets  that 
are  visible  on  the  high-speed  films.  In  some  of  the 
experiments,  oblong  (peanut-shaped)  styrofoam  packing  material 
was  stacked  in  a  racial  row  along  the  concrete  test  paa,  in  a 
direction  perpendicular  to  the  camera  line  of  sight.  The 
motion  ot  these  lightweight  chips  was  very  clear  on  the  film. 

On  the  basis  of  ooservations  from  these  experiments,  tne 
following  explanation  is  offered.  A  sketch  or  an  ideal, 
expanding  nemispher ical  cloud  is  shown  in  Figure  25(a).  The 
spray  is  distributee  uniformly  within  the  cloud  and  extends 
fully  to  tne  ground.  As  the  front  ot  the  spray  cloud  advances, 
momentum  is  adaeo  continuously  to  newly  and  previously  engulfed 
air  oy  drag  from  the  droplets.  This  sets  the  air  in  motion 
radially  wnicn  lowers  the  center-cloua  pressure  and  compresses 


the  ait  aheaa  of  the  expanding  front.  A  pressure  profile  is 
therefore  developed  which  rises  from  a  subatmospher ic  value  in 
tiie  center  of  tne  cloud  to  an  above-atmospheric  level  just 
ahead  of  the  front,  and  then  drops  back  to  atmospheric  pressure 
at  a  somewhat  greater  radius.  Air  flow  continues  to  be  outward 
from  the  center  until  the  retarding  internal  pressure  gradient 
becomes  sufficiently  steep  as  to  stop  the  air  movement  locally. 

This  process  cannot  apply  to  the  boundary  layer.  Near 
the  grouna  the  same  caaial  pressure  profile  is  impressed,  but 
the  air  in  this  region  experiences  an  additional,  radially- 
inward  acceleration  due  to  viscous  forces.  Gas  elements  in. the 
boundary  layer  will  therefore  be  stopped  before  the  air  above 
it  stops,  and  under  continued  action  of  the  pressure  gradient 
these  will  begin  to  flow  toward  the  center.  This  movement  is 
directed  upward  upon  collision  with  outward  moving  or  stagnant 
surface  layer  air,  and  a  separation  vortex  results.  In  the 
case  of  an  expanding  hemispherical  cloud  this  voLtex  would  be 
torroidal  in  shape,  as  sketched  m  Figure  25(b).  Inis  vortex 
is  oeiievea  to  be  responsible  for  cioua  eage  lifting. 

At  the  end  of  the  dispensing  period  the  pressure  gradient 
near  tne  center  of  the  cioua  is  no  longer  supported  by  a 
momentum  influx  rrom  the  spray  nozzles.  However  droplets  still 
in  motion  near  the  expanding  front  continue  to  drag  on  the  air 
in  that  region.  Air  begins  to  accelerate  towards  the  cloud 
center.  This  acceleration  is  greatest  in  the  air  layers  on  the 
underside  of  the  cloud  which  do  not  experience  droplet  drag- 
The  result  is  the  development  at  late  times  ot  a  large-scale 
torroidal  vortex,  as  depicted  in  Figure  25(c). 

If  the  dispenser  source  is  raised  above  the  ground, 
boundary  layer  effects  no  longer  influence  the  cloud.  However, 
large-scale  vortices  now  begin  to  form  immediately  after  the 
start  of  aispensing  because  the  reauced  center-cloud  pressure 
in  tnis  case  causes  air  to  be  drawn  under  the  cloud,  toward  the 


center.  Vertical  deflection  of  the  horizontally  directed  jets 
has  been  observed  in  this  case  as  these  jets  emerge  from  the 
nozzle  cluster.  When  the  large-scale  vortex  motion  is  suf¬ 
ficiently  severe,  secondary  vortices  are  visible.  The  result 
is  a  dish-shapea  cloud  as  sketched  in  Figure  25(d). 

The  distortions  described  are  induced  more  readily  in 
clouds  with  greater  jet  overlap,  i,e.,  smaller  0.  Air  paths 
between  the  fingers  of  spiked  clouds  provide  partial  relief  for 
some  of  the  pressure  gradients.  The  distortions  are  also 
aggrevatea  by  increased  spray  momentum  and  by  decreased  oroplet 
size,  bmaller  droplets  tend  to  follow  the  motion  of  vortices 
more  readily,  producing  more  pronounced  cloud  distortions. 

Since  the  size  of  the  droplets  generated  by  an  atomizing  jet 
decreases  with  increasing  jet  velocity,  such  increases  very 
rapidly  proouce  cloud  shape  distortions.  Some  control  might  be 
exercised  over  oroplet  size  by  appropriate  fuel  selection  or 
tnrough  the  use  of  additives,  but  for  a  given  fuel  the  jet 
velocity  must  be  limited.  The  relationship  between  the  maximum 
tolerable  jet  velocity  and  cloud  scale  is  not  presently  known 
since  all  of  the  experiments  to  date  have  been  with  4.6  m  (15 
ft)  radius  clouos.  This  is  an  important  area  which  requires 
further  investigation. 

4.5  DISPENSER  SCALEUP  REQUIREMENTS 

Although  all  questions  regarding  scaleup  of  the  impulsive 
ciouu  formation  concept  have  not  as  yet  been  resolved,  an 
attempt  has  been  made  with  what  is  known  to  anticipate  the 
fundamental  design  requirements  of  large  dispensing  systems.  A 
computer  program  has  been  developed  in  which  all  presently 
recognized  design  constraints  are  considered  and  an  optimum 
design  within  these  constraints  is  sought.  This  program  will 
be  updated  as  the  understanding  of  cloud  formation  improves. 

The  program  does  not  dictate  ail  design  details,  but  at  present 


ic  is  assumed  that  a  piston  is  incorporated  in  the  design  to 
act  as  an  interlace  that  separates  the  fuel  being  dispensed 
from  tne  driver  gases. 

An  outline  flow  diagram  of  the  dispenser  design  code  is 
given  in  Figure  26.  The  mass  of  fuel  and  the  cloud  radius 
required  for  a  specified  far-field,  effective  blast  source 
energy  are  first  calculated  directly  from  Eqs.  (30)  ana  (33;. 
The  individual  nozzle  size  required  is  then  obtained  from 


in  which  5  R^/D^  t^e  dimensionless  reach,  a  func¬ 
tion  generally  or  the  jet  velocity  and  when  sufficient 
fuel  is  dispensed. 

Tne  minimum  number  of  nozzles,  as  aiscussea  earlier, 
approximately  corresponds  to  a  configuration  in  whicn  all  of 
tne  radial  get  cylinders  overlap  and  are  tangent  on  the  cloud 
surface.  The  area  on  the  cloud  surface  that  is  allocated  to 
each  jet  in  this  array  is  that  of  a  spnerical  hexagon.  If 
approximated  as  a  planar  hexagon,  this  area  is  just  A  =  2/3  Rj 
wnere  R,.  is  tne  radius  of  the  spray  cylinder  inscribed  within 
the  nexagonal  space.  Equating  the  total  area  of  ail  such 
hexagons  to  the  surface  area  of  tne  hemispherical  cloud  leads 
to  the  approximation 


N 


N  . 
min 


(39) 


tor  the  minimum  number  of  nozzles,  wnere  R ^  =  rj/°n 
the  aimens ionless  get  radius. 
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SIMULATOR  DESIGN  CODE 
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Figure  26.  Schematic  flow  diagram  of  simulator  design  code 
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The  upper  limit  on  the  number  of  nozzles  is  at  present 
presumed  to  be  fixed  by  the  minimum  volume  that  can  be  impul¬ 
sively  injected  through  a  single  nozzle  without  reach  loss  aue 
to  tail  breakup  interference,  as  describea  in  Section  3;  that  is 
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Using  the  expression  for  Vj?m^n  given  by  Eq.  (26) ,  this  can  be 
written 


m„ 


N 


Nmax  0.428P£UD  Dn7//2 


(40) 


As  a  conservative  measure,  the  maximum  dispensing  velocity, 
UD  “  uDi'  (which  corresponds  to  the  peak  dispensing 
pressure  P3^  in  Figures  5  and  6)  is  presently  being  used  in 
this  equation. 


For  a  given  set  of  design  conditions,  the  maximum 
atmospheric  wind  velocity  that  can  be  tolerated  during  a  test 
is  an  important  consideration.  For  the  present  it  is  assumed 
tnat  the  cloud  will  be  detonated  at  the  formation  time,  tp. 
This  time  is  calculated  from  Eq.  (27)  with  UD  =  Ud2  (corres¬ 
ponding  to  P32  in  Figures  5  and  6)  to  be  conservative.  If 
the  sloud  displacement  due  to  wind  during  this  period  is  to  be 
no  mote  than  a  small  fraction  of  the  cloud  radius,  then 
the  maximum  permissible  wind  velocity  is 
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If  the  value  ot  a* _ calculated  in  this  way  is  not  realistic, 

wmax 

either  tF  or  Rc  must  be  changed. 

During  the  dispensing  perion  the  jet  velocity  muse  be 
high  enough  for  atomisation,  but  it  cannot  be  so  high  as  to 
induce  serious  vortex-induced  cloud  shape  distortions.  Thus, 
there  is  a  minimum  and  a  maximum  jet  velocity.  Since  for  all 
dispensing  systems  presently  being  considered,  UD  >  ,  th;‘-s 

means  that 


max 


and 


Ur 


u. 


min 


(42.1) 


(42.2) 


are  required.  The  lower  jet  velocity  limit  UDmin  is  associ¬ 
ated  with  the  knee  of  the  curve  of  reach  as  a  function  of 
dispensing  velocity  (see  Figure  21) .  The  upper  limit,  Ud  , 

iuaX 

has  not  been  well  established  as  yet  but  it  appears  to  be  net 

far  removed  from  Un  •  . 

“min 


As  presently  conceived,  the  dispensing  system  to  produce 
large  clouds  would  consist  of  a  cluster  of  Np  modules,  close- 
packed  at  the  center  of  the  test  arena  into  a  circular  region 
of  overall  radius  R . •  .  Each  module  would  comprise  a  piston 
of  diameter  Dp  in  a  cylinder  assembly  of  diameter  DpN< 

Behind  eacn  piston  a  time-varying  free  volume  V  is  co  be 
provided.  The  total  instantaneous  free  volume  for  the  entire 
system  is  VQ  -  NpV. 

!jj|  A  particular  i  jaule  would  dispense  an  amount  of  fuel 

mF  =  mF/NP  through  NN  =  N^/Np  nozzles  into  a  designated  .sector 
of  the  cloud.  The  modules  would  be  identical  in  construction; 
chat  is,  all  of  the  pistons  would  have  identical  diameter,  Dp, 
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mass,  nip,  ana  length,  Lp.  The  total  nozzle  cross-sectional 

area  that  is  driven  by  one  piston  is  A^  *  ^nAn'  wnere 

A.,  is  the  area  of  one  nozzle. 

N 


If  all  of  the  nozzles  are  to  be  the  same  size  and  if  they 
are  all  driven  at  the  same  pressure,  they  will  have  the  same 
reach.  If  the  dispenser  cluster  radius  were  comparable  to  the 
cloud  raaius,  the  cloud  from  these  nozzles  could  not  be 
hemispherically  shaped.  To  minimize  such  cloud  shape 
aberrations,  it  is  therefore  necessary  that  R^g  -  C3RC, 
where  is  a  small  number. 

The  dynamics  of  individual  modules  are  predicted  by 
Equations  (11)  through  (17).  (Appropriate  notational 
modifications  are  needed  to  adapt  those  equations  to  the 
present  context.)  The  intial  conditions  are  Pp  =  P-dq,  up 
=  0,  V  =  V0,  and  m^  =  mp  at  t  =  0.  From  an  integration  of 
these  equations  the  time  t-^  to  reach  peak  dispensing  velocity 
UD^  and  also  the  dispensing  time  tD  (which  corresponds  to 
m^  *  0,  UD  =  Ud2)  are  determined.  To  confine  the 
degradation  of  jet  quality  that  is  due  to  piston  acceleration 
to  a  minimum  initial  period  requires 

h  -  c2tD 

where  C2  is  small.  It  is  notea  that  Eq.  (16)  applies  only  to 
systems  in  which  the  freevolume  expansion  is  isentropic.  It 
would  not  apply,  for  example,  to  a  tailored  gas  generator  or  to 
a  choked-orif ice  dispensing  system.  These  designs  would 
require  separate  treatment. 


In  any  dispenser  with  a  driving  pressure  that  decreases 
during  the  dispensing  process,  successive  jet  segments  possess 
progressively  lower  velocities.  This  tends  to  stretch  the 
jet.  That  is,  at  the  nozzle  exit  9p/9x  =  0  along  the  jet  axis, 
ana  so  the  inviscia  equation  of  (one-aimensional)  motion  near 
that  point  is 
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If  SOp/St  <  0  as  the  result  of  aAP/dt  <  0  (see  Eq.  (11)), 
then  9Ud/3x  <  0,  implying  a  stretching.  Since  the  cohesive 
resistance  to  tensile  forces  in  a  liquid  is  very  small,  this 
effect  will  contribute  to  accelerated  jet  disruption  if  it  is 
too  large.  No  satisfactory  criterion  for  the  maximum  tolerable 
velocity  gradient  has  as  yet  been  proposed  ana  no  experiments 
have  as  yet  addressed  this  question  directly.  However,  to 
proviae  some  ceiling  on  the  nozzle-exit  velocity  gradient  (for 
purposes  of  developing  the  dispenser  design  code) ,  the  velocity 
difference  UQ  resulting  from  a  velocity  gradient  that  acts 
over  a  distance  equal  to  half  the  jet  reach  is  arbitrarily 
limited  to  no  more  than  the  average  dispensing  velocity  Ug. 
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This  velocity  gradient  can  be  reduced  by  increasing  the  free 
volume  (which  lowers  Au^)  ,  but  a  very  large  freevolume  tank 
would  be  costly.  An  understanding  of  the  effects  of  a  velocity 
gradient  on  jet  breakup  would  be  helpful  in  establishing  the 
smallest  usable  free  volume. 


To  be  complete,  the  list  of  dispenser  design  considera¬ 
tions  should  include  hardware  and  installation  cost  equations. 
The  design  which  satisfies  all  of  the  engineering  constraints 
at  minimum  cost  would  be  sought.  Although  this  has  not  as  yet 
oeen  done,  some  limitations  based  on  economy  nave  been 
tentatively  incorporated  into  the  code.  To  minimize  the  total, 
numoer  of  nozzles,  Eqs .  (39)  and  (40)  are  set  equal.  This 
effectively  establishes  an  upper  bound  on  which  replaces 

chat  expressed  by  Eq.  (42)  if  it  is  smaller.  In  addition,  the 
diameter  of  the  individual  dispenser  module  pistons  is 
presently  limited  in  the  program  to  29.2  cm  (11.5  in.),  and  the 
length  of  eacn  piston  cylinder  is  limited  to  2.5  m  (100  in.). 
These  limits  were  established  because  it  was  felt  that  the 
costs  of  fabricacing  the  piston  and  cylinder  would  increase 
very  rapidly  for  larger  sizes.  The  freevolume  size  is  also 
minimized  by  setting  Ud2  =  This  replaces  Eq.  (42.2). 


The  program  as  written  iterates  on  Pp0,  VQ,  and 
(which  must  be  an  integer)  until  a  solution  satisfying  all 
constraints  is  round,  if  this  cannot  be  done,  the  constraints 
that  are  violated  are  indicated. 
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Preliminary  calculations  using  the  design  coae  have  been 
carried  out  at  scales  corresponding  to  0.267,  20,  ana  1000  ton 
nuclear  equivalents,  with  propylene  oxide  as  the  assumed  fuel. 
(The  smallest  of  these  corresponds  approximately  to  the  scale 
ot  early  experimental  tests  conoucted  by  Systems,  Science  ana 
Software.)  The  program  output  is  presented  in  Table  2.  In 
the  calculations,  =  1500  and  Rj  =  50  were  assumed. 

Also,  C]_  *  0>2  *  c3  *  0.05  was  used.  It  was  assumed  that  DpN 
»  1.5  Dp,  and  also  that  mp  =  ppApLp,  where  Lp/Dp  =  1.91  and 
«  1150  kg/m3  ( t2  lb/ft3)  are  based  on  the  design  of  the 
linear  dispenser.  A  value  of  np  =  0.74  was  used  as  an 
estimate  of  the  effective  blastwave  efficiency  of  the  propylene 
oxide  FAR  clouds  (see  Section  5.2.3).  The  freevolume  (driving) 
gas  was  assumed  to  have  a  specific  heat  ratio  y  *  1.4. 

The  calculations  that  are  summarized  in  Table  2  were 
performed  after  the  linear  dispenser  had  already  been  fab¬ 
ricated.  For  that  reason  the  calculated  dimensions  do  not 
coinciae  completely  with  the  actual  design.  Values  of  UDm^n 
were  cruaeiy  estimated  for  generating  the  calculation.  The 
results  in  cne  table  should  accordingly  be  regarded  as  very 
rough  approximations. 


One  observation  from  these  approximations  is  that  the 
peak  dispensing  velocities  specified  for  the  two  smaller  scales 
are  quite  high.  Cloua  distortion  effects  due  to  induced  air 
motions  are  therefore  almost  certain  to  occur  at  the  0.267  ton 
scale,  anu  would  be  likely  at  the  20-ton  scale  if  the  values  ot 
Pp  listed  were  to  be  used.  The  dispensing  velocity  required 
at  the  1  KT  scale  is  more  reasonable  and  at  the  same  time  the 
velocity  gradient  is  also  quite  low.  Out  a  very  large  free 
volume  would  be  needed  in  that  case.  Also,  at  the  1  KT  scale 
the  overall  dispenser  cluster  radius  is  nearly  9%  of  the  cloud 
radius.  If  cylinders  longer  than  2.5  m  (100  in.)  were 
permitted  however,  the  dispenser  cluster  radius  could  be 
reduced. 
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Table  2.  Approximace  simulator  system  design  at  three  scales 
using  propylene  oxiae  as  fuel. 

Scale  (tons) 

0.267  20  1000 


CLOUD  CHARACTERISTICS 

Radius,  Rc(m) 

Total  fuel  mass,  m_(kg) 

F  3 

Total  condensed  fuel  volume,  V. ,  (m  ) 

Cloud  formation  time,  tj.(s) 

DISPENSER  MODULE  CHARACTERISTICS 

Piston  diameter,  Dpt1”) 

Piston  mass,  in  (kg) 

? 

Overall  cylinder  length,  L„_(m) 

cp 

Module  overall  diameter,  Dpn(m) 

Fuel  mass  per  module,  m^kg) 

Nozzles  per  moduli  ,  Nn 
Fuel  mass  per  r.o.zle,  uL,(kg) 

Single  nozzle  diameter,  Dn(m) 

Individual  jet  spray  column  diameter,  D. (a) 

~  3  J 

Free  volvrne  per  module,  V3  (m  ) 

Paak  dispensing  velocity,  ^o^tm/s) 

Time  to  reach  0.9^^,  t^  ( 3 ) 

Dispensing  time,  tD(s) 

Dispensing  velocity  at  tD,  ud2(,''/s) 

Ratio  t^/t0 

Initial  driving  pressure,  PpQ(MPa) 

DISPENSER  SYSTEM  CHARACTERISTICS 

Number  of  modules,  N_ 

P 

Dispenser  system  radius, 

Ratio  RdijJ/Rc 

Total  number  of  nozzles,  N_ 

3  n 

Total  free  volume,  V  (m'  ) 


4.43 

18.7 

68.9 

22.7 

1704.0 

85,200.0 

0.0265 

1.99. 

99.7 

0.0601 

0.324 

1.71 

0.292 

0.292 

0.292 

43.0 

43.0 

43.0 

0.952 

2.2 

2.5 

0.439 

0.439 

0.439 

22.7 

94.7 

112.0 

1632.0 

91.0 

2.14 

0.0139 

1.04 

52.2 

0.00296 

0.0125 

0.0459 

0.310 

1.31 

4.82 

0.00413 

0.119 

7.93 

138.0* 

67.7 

35.3 

0.00168 

0.0138 

0.0122 

0.0305 

0.185 

1.05 

50.0 

45.0 

33.0 

0.055* 

0.0743* 

0.0117 

17.3 

2.35 

0.636 

1.0 

18.0 

762.0 

0.219 

0.930 

6.05 

0.0494 

0  „  0  4  9  7 

0.0878* 

1632.0 

1632.0 

1632.0 

0.00413 

2.15 

5080.0 

*Does  not  satisfy  specified  constraint 
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SECTION  5 


FAE  SIMULATION  OF  FAR-FIELD  NUCLEAR  AIRBLAST 

5.1  CHARACTERISTICS  OF  NUCLEAR  AIRBLAST  IN  THE  FAR-FIELD 

The  degree  of  correspondence  or  fidelity  between  FAE  and 
surface  nuclear  airblasts  is  usually  discussed  in  terms  of 
positive-phase  blastwave  characteristics.  These  are  peak 
static  overpressure,  peak  dynamic  pressure,  and  the  associated 
positive-phase  durations  and  impulses,  all  as  functions  of 
radius;  and  the  static  overpressure,  dynamic  pressure,  ana  the 
associated  impulses  as  functions  of  time  at  fixed  radii. 

However,  a  complete  set  of  universally  accepted  nuclear 
data  that  could  be  used  for  a  comparative  basis  is  not 
presently  available.  Sources  or  experimental  ana  theoretical 
nuclear  aata  (as  well  as  some  TNT  airblast  data)  are  listed  in 
References  41-62.  Some  of  the  aata  in  these  sources  are 
conflicting.  Presumably  a  certain  amount  of  the  conflict  is 
due  to  differences  in  instrumentation,  data  reduction,  and 
cur  efitting  procedures. 

For  example,  mucn  of  t  ,.e  nuclear  static  overpressure  data 
were  obtained  with  variable  reluctance  or  strain  gage 
transducers  in  conjunction  with  magnetic  tape  recording 
systems.  The  overall  risetime  capability  of  these  systems  was 
on  the  order  of  one  millisecond.  Self-recording  mechanical 
gages  (developed  by  Ballistic  Research  Laboratory)  with 
risetimes  of  3-5  ms  were  also  used.  Characteristic  samples  of 
semi-smootheo  nuclear  data  are  given  in  Figures  27  and  28. 

High  frequency  noise  has  been  filtered  out  of  these  samples, 
but  lower  frequency  pressure  excursions  (which  carry  real  or 
apparent  impulse)  have  been  retained. 


In  reporting  peak  overpressures,  the  first  crest  of  the 
initial  pressure  rise,  as  recorded  by  the  instrumentation,  was 
used  by  Reference  44.  On  the  other  hand,  in  Reference  41  an 
extrapolation  procedure  to  compensate  for  the  equipment  rise- 
times  was  followed.  On  a  semi-logarithmic  plot  with  overpres¬ 
sure  on  the  logarithmic  scale  ana  time  on  the  linear  scale,  a 
straight  line  was  drawn  through  the  data  points  near  t  =  0. 
This  was  extrapolated  to  time  zero  in  order  to  infer  the  peak 
overpressure . 


Experimental  r.uclear  airblast  data  for  measurements  ether 
than  static  overpressure  (ana  impulse)  are  more  scarce  and 
contain  more  scatter.  This  is  particularly  true  of  the  time 
variation,  positive-phase  duration,  and  impulse  associated  witn 
dynamic  pressure. 

In  the  current  investigation,  the  nuclear-data  curvefits 
given  by  brode  (Reference  52)  were  used  as  the  reference 
against  which  FAE  airblasts  ate  compared.  Examples  of  overpres¬ 
sure-time  cucvetits  from  this  reference  ate  incluaea  in  Figures 
27  ana  28.  These  are  fully  smoothed  curvefits,  in  the  sense 
that  all  excursions  from  tne  average  profiles  have  been 
excludea.  The  Brode  nuclear  curves  for  peak  static  over¬ 
pressure,  peak  dynamic  pressure,  static  overpressure  impulse 
and  positive  pnase  duration  appear  in  Figures  29-32,  The  peak 

static  overpressure  AP__„  and  peak  dynamic  pressure,  q  , , 

luox  max 

have  been  made  dimensionless  in  these  plots  with  respect  to 
P^  *»  1.01  x  IQ5  Pa  (14.7  psi)  .  The  positive  phaso  duration 
has  been  normalized  with  respect  to  the  characteristic  time 
R/c ^ ,  where  R  is  tne  range  from  the  explosion  center  and 
=  341  m/s  (1116  ft/s)  is  the  sea  level  speed  of  sound  in  air. 
Impulse  is  made  dimensionless  with  respect  to  both  the  time 
R/ajL  and  pressure  p^.  Finally,  the  range  is  expressed  as 
the  dimensionless  variable 


R  = 


104 


(44) 


kT 


0 
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R/R 

o 

Figure  30.  Nuclear  Surface  Burst  Reference  Curve  for 
Positive-Phase  Static  Overpressure  Impulse 
Versus  Range  (Reference  52) . 
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Here  the  character istic  radius  RQ  is  defined  by 


(45) 


where 

'  2,  hemispherical  geometry 

Vs*  - 

„  4,  spherical  geometry 

The  energy  EQ  in  this  expression  is  the  effective  or 
apparent  blast  yield,  which  for  a  nuclear  explosion  is 
designated  by  EoE  Effective  yield  is  a  concept  that  was 

discussed  in  Section  4.1.  Values  of  Eqe  are  usually  given 
in  terms  of  equivalence  with  respect  to  TNT  airblasts.  The 
effective  blast  yield  varies  with  range  and  with  the  variable 
chosen  as  a  basis  for  comparsion.  However,  over  a  limited 
range  interval  the  variation  m  Eq^,  ^  may  be  small  enough 
that  a  single  average  value  can  be  used. 

It  has  been  found  useful  to  define  the  effective  airblast 
efficiency  (or  effectiveness)  of  a  nuclear  explosion  by 


(46) 


This  definition  parallels  Eq.  (29)  which  was  written  for  a 
fuel-air  explosion.  Based  on  overpressure  decay  in  the  tar 
field,  a  suitable  average  value  of  is  approximately  0.5. 

This  value  has  been  used  for  the  purposes  of  this  reporrt.  For 
a  1  KT  nculear  explosion,  the  total  energy  release  is  E$T  N 
=  4.18  x  lO1^  J  (3.96  x  10^  Btu)  (Reference  45)  and  so  for 


v  *  2,  R*  ,,  =  149  m  (489  ft)  ,  assuming  P,  *  1.01  x  105 

hi  f  N  i- 

Pa  (14.7  psi) .  (The  asterisk  is  used  to  designate  the  1  KT 
level.)  This  value  of  R  ,  whicn  corresponds  to  the  effective 

Q 

blastwave  energy  of  a  1  KT  nuclear  explosion,  was  usea  in 
preparing  Figures  29-32.  To  the  extent  that  cuoe-root  scaling 
as  implied  by  Eq.  (45)  can  be  considered  valid,  these 
dimensionless  figures  then  apply  to  any  scale  of  nuclear 
explosion. 


5.2  EXPERIMENTAL  INVESTIGATIONS 
5.2.1  Test  Setup  and  Procedure 

The  experimental  facility  presently  being  usaa  to  inves- 
vestigate  small-scale  FAE  airblasts  (^ 1/4  ton  nuclear  equiv¬ 
alent)  is  described  in  Appendix  A.  The  facility  comprises  two 
perpendicular  concrete  runways  extending  from  the  explosion 
center.  Along  these  runways  there  are  instrumentation  troughs 
located  as  listed  in  Table  A.l.  The  principle  instrumentation 
consists  of  stagnation  and  static  pressure  measurements  at  each 
trough.  In  addition,  high-speed  films  are  taken.  FAE  clouds 
are  impulsively  generated  at  the  explosion  center  using  the 
linear  dispenser  (described  in  Section  2)  with  a  multi-nozzle 
hemispherical  head  of  the  type  shown  in  Figure  22.  Tests  are 
normally  conducted  in  the  early  morning  hours  to  minimize  wind 
effects.  Following  a  suitaole  aelay  arter  impulsive  cloud 
formation,  the  cloud  is  detonated.  Detonation  is  initiated 
with  a  small  (50-500  gm)  H.E.  charge  suspended  directly  above 
the  dispenser  head  at  a  minimum  height.  The  minimum  height  is 
determined  by  considering  the  number  of  jets  that  are  blocked 
by  the  H.E.  cnarge,  or  oy  the  potential  for  damage  to  the 
dispenser . 


5.2.2  Data  Reduction 


Pressure  data  gathered  during  a  test  are  digitized  ana 
the  bulk  of  the  aata  reduction  is  then  machine  processed.  A 
schematic  flow  chart  outlining  the  data  reduction  procedure  is 
given  in  Figure  33. 


At  each  instrumentation  position,  the  time  or  olastwave 

arrival,  the  peak  static  overpressure,  and  the  peak  reflected 

stagnation  pressure  are  measured.  The  positive-phase  impulse 

is  obtained  from  the  static  pressure  record  at  that  location  by 

integration.  From  these  four  measurements  an  average  value  of 

Kq  is  determined.  This  is  done  as  follows.  A  curve  is 

smoothed  through  the  time-of -arr ival  versus  range  data.  From 

this  curve  the  wave  front  velocity  li  at  the  instrumentation 

s 

location  in  question  is  determined.  The  shock  Mach  number  M 

s 

is  computed  from  =  Us/a^.  Separately,  the  Mach  number 
is  calculated  rrom  tne  peak  static  overpressure  AP  using 


Ms  -  1  +  fr) 


(47) 


and  from  the  peak  reflected  head-on  overpressure  afk  through 


=  2  ( P  -  1)  + 


(Y  +  1)  (P  -  ir 
2Y  +  (Y  -  1)  (P  -  1) 


where 


t 


and  Y  =  1.4  is  the  specific  heats  ratio  for  air.  The  value  of 

P  is  determined  from  Eq.  (48;  ana  then  M  is  calculated  from 

s 
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DATA  REDUCTION  PROCEDURE 


Figure  33.  Schematic  flow  diagram  of  data 
reduction  program. 


Eq.  (47) .  The  three  values  of  ate  averaged.  Eq.  (47)  is 
next  used  to  infer  a  mean  of  value  of  AP_=V.  The  corre- 
sponding  characteristic  radius  KQ  can  then  be  obtained  from 
Figure  29  since  the  measurement  range  R  is  of  course  known.  In 
the  same  way,  Rq  is  obtained  from  the  measured  positive  phase 
impulse,  by  reference  to  Figure  30.  These  two  values  of  Rq 
ate,  then  averaged.  The  same  proceaute  is  followed  at  all  other 
instrumentation  locations  and  a  final,  overall  average  value 
for  Rq  is  generated.  This  is  the  effective  characteristic 
radius  for  the  FAE  airblast  with  respect  to  the  Erode  nuclear 
curves  and  i3  desiqnateo  Ror,  _.  The  corresponding  effective 

f  £ 

olastwave  energy  is  EOT,  calculated  from  Eq.  (45)  .  The 
measured  peak  overpressure,  positive  phase  duration  ana 
impulse,  and  the  dynamic  pressure  can  then  be  plotted  versus 
R/ROE  and  compared  directly  with  Figures  29-32.  The 
overpressure  as  a  function  of  time  is  also  plotted  and  compared 
with  tne  Brode  curvef it  (Reference  52) . 

Dynamic  pressure  is  derived  from  the  stagnation  and 
static  pressure  measurements.  The  dynamic  pressure 

q  5  1/2  Pu2 

can  readily  be  shown  to  be 

q  ■  1/2  YPM2  (49) 

for  a  perfect  gas,  where  P  is  the  absolute  static  pressure  (not 
overpressure)  and  the  local  flow  Mach  number  is  M  -  u/a.  The 
peak  value  of  q,  that  is,  q2  or  gmax»  can  be  calculated 
directly.*  This  is  because  =  P-,  =  APm=„  +  P.  is  measured, 

*Point  ”2"  is  used  to  aesignate  conditions  just  being  the  blast- 
wave  shock  front  in  theoretical  work  described  later  in  this 
report. 


and  the  flow  Mach  number  just  behind  the  blastwave  shock,  M2 , 
is  related  to  the  shock  Mach  number  by 

2(M?  +  1) 

M2  =  — -  - - -  - - -  .  (50) 

/(2YM?  -  (Y  -  1)  1  [2  +  (Y  -  D  nl  ] 

The  shock  Mach  number  is  determined  from  Eq.  (47),  so  the 

stagnation  pressure  measurement  is  not  needed* to  obtain 

q  „  .  However,  at  later  times  the  local  flow  Mach  number  at 
wax 

a  fixea  radial  position  is  not  related  to  the  shock  front 
strength.  The  stagnation  probe  pressure  measurement  is, 
therefore,  needed  in  order  to  obtain  q(t) .  It  the  flow  into 
the  stagnation  pressure  probe  is  subsonic,  the  pressure  as 
measured  by  that  probe,  Pq, is  equal  to  the  f reestream  stag¬ 
nation  pressure,  PQ.  The  flow  Mach  number  (needed  in  Eq. 

(49))  is  in  this  case  found  from 


where  P  and  P  are  absolute  pressures, 
o 

On  the  other  hand,  if  the  flow  into  the  stagnation 
pressure  probe  is  supersonic,  then  the  measured  pressure  Pq  is 
not  the  actual  stagnation  pressure  PQ.  At  the  first  moment 
following  the  passage  of  the  blastwave  leading  shock,  the 
pressure  at  the  stagnation  probe  reaches  its  peak  reflected 
value;  i.e.,  Pq  -  PR  at  that  moment.  Within  a  few  micro¬ 
seconds  (depending  on  the  probe  dimensions) ,  the  reflected 
shock  moves  upstream  of  the  probe  face  and  a  quasi-stationary 
bow  shock  is  formed.  Thereafter,  the  absolute  static  and 
stagnation  prooe  pressures  are  related  by 


> 


1 

Y-l 


(52) 
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It  can  be  demonstrated  chat  if  at  any  time  the  measured 

ratio  P/Pq  exceeds  0.526  (for  y  *  1.4),  then  the  flow  at  that 

moment  is  subsonic.  Conversely,  if  P/Pd  <  0.528,  then  the  flow 

is  supersonic.  This  criterion  is  invoked  to  determine  which  of 

Eqs.  (51)  or  (52)  should  be  used  to  compute  M  at  a  given 

instant,  once  M  is  known,  Eq.  (49)  is  used  to  calculate  q. 

The  moment  in  time  at  which  q  -  0  (M  ■  0,  u  »  0)  defines  the 

positive  phase  dynamic  pressure  duration,  t  . 

M 

The  effective  airblast  efficiency  of  the  FAE  is 
calculated  by  combining 


(53) 


with  Eqs.  (28)  and  (29),  which  gives 
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It  should  be  notea  that  other  investigators  have  defined 
blastwave  efficiencies  aifferently.  For  example  (Reference 
63),  a  blastwave  efficiency  can  be  definea  as  the  ratio  of  the 
actual  mechanical  positive-phase  work  done  by  the  expanding, 
burned  cloud/air  interface  to  the  total  conu-ustion  energy  based 
on  a  heat  of  combustion.  However,  by  itself  this  particular 
definition  of  efficiency  cannot  be  used  to  determine  the  amount 
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of  fuel  needed  to  simulate  the  airblast  characteristics  of  a 
point  blast  or  a  nuclear  source  of  specified  strength.  In 
conjunction  with  this  efficiency,  an  "equivalent  blast  energy 
ratio"  must  be  used  in  order  to  calculate  the  necessary 
quantity  of  fuel.  When  point  source  ana  FAE  decay  curves  are 
matched  over  the  same  range  of  biastwave  radii,  the  product  of 
the  interface-work  efficiency  and  the  "equivalent  blast  energy 
ratio"  should  yield  a  value  that  is  approximately  equal  to  the 
effective  blast  efficiency  defined  in  this  report. 

An  explosion  scale  factor  with  respect  to  a  1  KT  nuclear 
explosion  is  also  calculated  as  part  of  the  data  reduction. 
Cube-root  seeling  implies,  for  fixed  atmospheric  conditions, 
that 


where  again  the  suoscript  "L\i" 
designates  "FAE".  The  scale 
oy 


S  5 


designates  "nuclear"  and 
factor,  s,  is  accordingly 


II  p  ;i 

oef ined 


(55) 


That  is,  the  scale  factor  is  written  in  terms  of  the  effective 
biastwave  energies  of  the  fuel-air  and  1  KT  nuclear  explosions. 
This  can  also  te  expressed  as 


s 


by  using  Eq.  (45) . 


(56) 
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5.2.3  Experimental  Results  and  Discussion 

As  noted  in  Appendix  A,  preliminary  experiments  were 
carried  out  witn  somewnat  less  sophistication  than  the  more 
recent  work  described  in  the  present  section.  The  early 
results,  however,  were  sufficient  to  demonstrate  the  funda¬ 
mental  feasibility  of  nuclear  airblast  simulation  by  means  of 
two-phase  fuel-air  explosions,  when  the  FAE  cloud  is  generated 
from  a  central  source  by  means  of  multiple  nozzle,  impulsive 
hydraulic  fuel  injection.  Those  experiments  are  discussed  in 
detail  in  Reference  5  and  only  a  brief  summary  is  given  here. 
Approximately  stoichiometric,  4.6  m  (15  ft)  radius  hemisphe¬ 
rical  propylene-oxide/air  clouds  were  detonated  in  the  experi¬ 
ments.  Relatively  long  delay  times  (up  to  2500  ms)  between  the 
start  of  the  fuel  dispensing  and  cloud  detonation  were  used. 
These  delays  were  in  considerable  excess  of  the  minimum  cloud 
formation  times  (Eq.  (27))  so  that  due  to  settling  of  the  fuel 
sprays  the  clouds  were  usually  flatter  than  the  optimum 
hemispherical  shape. 

Characteristic  records  of  static  overpressure  versus  time 
from  one  of  these  tests  are  given  in  Figures  34-37.  The  fully 
smoothed  nuclear  curvefits  of  Brode  are  superimposed  on  the 
figures  for  reference.  The  corresponding  experimental  static- 
overpressure  impulse  plots  and  smoothed  nuclear  curves  are 
given  in  Figures  38-41,  The  degree  of  repeatability  between 
tests  can  observed  in  Figures  42  and  43,  The  symmetry  of  the 
blastwaves  produced  by  the  FAE  source  is  represented  by  the 
samples  in  Figures  44  and  45.  In  Figures  46  and  47,  the  peak 
overpressures  and  total  positive  phase  impulses  versus  range 
are  plotted  and  compared  against  tne  nuclear  curves.  Also 
shown  on  these  two  figures  are  tne  theoretical  point-source 
blastwave  curves  (Reference  62) ,  which  coincide  closely  with 
the  ttrode  nuclear  curves  in  this  range  of  pressures. 
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Figure  35. 
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Comparison  of  Measured  Overpressure  as  a 
Function  of  Time  with  a  Nuclear  Blast  Wave¬ 
form.  The  experimental  data  were  taken  from 
Test  2536  (t„  is  arrival  time) . 
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Figure  36. 
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Comparison  of  Measured  Overpressure  as  a 
Function  of  Time  with  a  Nuclear  Blast 
Waveform.  The  experimental  data  were  taxen 
from  Test  2536  (ta  is  arrival  time) . 
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Function  of  Time  with  a  Nuclear  Blast  Wave¬ 
form.  The  experimental  data  were  taicen  from 
Test  2536  (ta  is  arrival  time; . 
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Figure  43.  Plot  of  measured  static  overpressure  as  a  func¬ 
tion'  of  time  at  the  24.4  m  (80  ft)  station  from 
three  separate  experiments  (t  is  ‘arrival  time) 


Figure  44 . 


Plot  of  measured  static  overpressure  as  a  func 
tion  of  time  from  a  single  experiment.  The 
gauges  were  both  located  at  a  range  of  12.2  m 
(40  ft)  but  were  separated  by  90  degrees.  The 
experimental  data  were  taken  from  Test  2571  (t. 
is  arrival  time) .  ‘ 
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Figure  45. 
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Plot  of  measured  static  overpressure  as  a  function 
of  time  taken  from  a  single  experiment.  Both 
gauges  were  located  at  range  of  12.2  m  (40  ft) 
but  separated  by  90  degrees.  The  comparison  is 
similar  to  that  shown  in  Figure  15  except  that 
the  data  plotted  are  from  Test  2574  (t^  is 
arrival  time) . 
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The  general  agreement  between  the  FAE  and  nuclear  daca  as 

well  as  the  reasonaole  degree  of  symmetry  ana  repeatability 

were  encouraging  m  these  early  tests.  However,  the  strain 

gauge  pressure  transducers  were  too  slow  to  yield  accurate  peak 

overpressure  values  at  tne  small  experimental  scale.  Those 

plotted  in  Figure  46  were  accordingly  considered  to  be  low  by 

some  unaer termined  amount.  Rather  than  extrapolate  to  infer 

peak  values,  the  peaks  from  the  raw  data  were  reported  as 

AP  .  dcale  factors  and  effective  blast  coupling  e£- 
max 

ficiencies  were,  therefore,  determined  from  the  positive  phase 

impulse  measurements  alone,  since  impulse  is  less  sensitive  to 

inaccuracies  in  AP  .  The  averages  of  these  two  parameters 

max  3 

among  all  tests  were  s  =  16  and  n„  =  0.74. 

r 

In  addition  to  rounded  peaks,  the  static  pressure  data 
also  exhibit  numerous  excursions.  These  precluded  credible 
dynamic  pressure  and  positive-phase  duration  determinations 
which  are  quite  sensitive  to  small  measurements  errors.  It  is 
true,  oi  course,  that  raw  nuclear  airblast  data  are  also 
irregular  in  time  (see  Figures  27  ana  28).  However,  it  was  not 
known  if  the  pressure  excursions  in  the  FAE  tests  were  inherent 
to  the  explosion  process.  It  was  considered  likely  that  the 
excursions  were  at  least  in  part  the  result  of  reflections  from 
surface  irregularities  in  the  vicinity  of  the  transducers  or 
reflections  from  ground  slopes  surrounding  the  test  area  that 
was  in  use  at  the  time.  Questions  of  this  nature  were  con¬ 
sidered  sufficient  justification  to  warrant  the  improvement  and 
relocation  of  the  test  facility. 

A  limited  numoer  of  FAE  tests  utilizing  the  new  facili¬ 
ties,  hardware,  instrumentation,  and  data  reduction  software 
wer  ’  conducted.  These  tests  were  intended  to  exercise  the  new 
equipment  and  procedures,  and  to  attempt  detonation  of  clouds 
after  minimal  rormation  time.  A  problem  was  experienced  in 
this  test  series  with  wnat  appeared  to  be  oscillatory  noise 
caused  by  vibration  of  the  transducer  mounting  plates  (see 


Figures  a-3  and  A-4) .  Unfortunately  this  problem  was  not  cor¬ 
rected  by  simple  modifications  of  the  plates  (such  as 
installing  the  transducers  in  aelrin  or  teflon  inserts). 

A  sample  of  the  raw  data  taken  from  a  test  (CD-6)  in  this 
ser.es  appears  in  Figure  48.  The  semi-smoothed  profile  used 
for  data  reduction  is  super imposea  on  the  raw  data  in  this 
figure.  Experiment  CD-6  was  a  test  in  which  a  nominally 
stoichiometric,  propylene-oxide/air  cloud  was  detonated.  The 
linear  dispenser  was  used  with  a  600-hole  nozzle  head.  The 
peak  fuel  dispensing  pressure  was  about  1.3  MPa  (190  psi) . 
Detonation  of  the  cloud  was  initiated  by  a  0.5  kg  (1.1  lb)  ball 
of  H.E.  (C-4)  mounted  0.38  m  (15  in.)  above  the  nozzle  head. 

The  initiator  was  fired  0.6  s  after  the  start  of  fuel 
dispensing,  at  which  time  the  appearance  of  the  fuel/air  cloud 
was  similar  to  that  in  Figure  23(a). 

The  data  from  test  CD-6  was  fully  reduced.  Character¬ 
istic  explosion  lengths  as  inferred  from  peak  static  over¬ 
pressure  and  from  overpressure  impulse  measurements  are  plotted 
in  Figure  49.  Tne.  ar ithmetic  average  of  these  is  RQ  =  7.54  m 
(24.7  ft),  which  leads  to  n c  =  0.392  and  s  =  19.8.  The  iow 

r 

efficiency  and  high  scale  factor  are  attributed  to  the  short 
time  allowed  for  cloud  formation.  That  is,  the  fuel  spray  jets 
had  not  spread  sufficiently  to  engulf  the  air  between  them  at 
the  cloud  boundaries.  The  short  cloud  formation  times  were, 
however,  consistent  with  the  requirement  that  the  FAE  simulator 
be  capable  of  operation  in  moderate  winds.  Improvements  in 
technique  should  permit  fast  cloud  formation  with  increased 
effective  biastwave  efficiency. 

The  reduced  data  from  test  CD-6  are  presented  in  Figures 
50-61.  The  peak  static  overpressure,  overpressure  impulse, 
overpressure  positive  phase  duration,  and  dynamic  pressure  are 
plotted  as  functions  of  range  in  Figures  50-53.  Semi-smoothed 
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Figure  52.  Static  overpressure  positive  phase  duration  versus 
range  for  test  CD-6 . 
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Figure  55.  Static  overpressure  versus  time  during  test  CD-6  at  R/R 
=  1.31,  transducer  position  no.  3  (t  is  arrival  time) . 


0.1 


0.2 


0.3 


(t-ta)ai/R 


Figure  56.  Static  overpressure  versus  time  during  test  CD-6  at  R/R 
*  1.84,  transducer  position  no.  4  (t  is  arrival  time).' 
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Figure  61.  Dynamic  pressure  versus  time  during  test  CD-6  at  R/R 
-  2.51,  transducer  location  no.  5  (t  is  arrival  tim 
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static  overpressure  versus  time  curves  are  given  in  Figure 
54-59  at  the  various  transducer  locations.  Fully  smoothed 
nuclear  form  curves  are  superimposed  on  these  rigures  for 
compar ison . 

Semi-smoothed  dynamic  pressure  versus  time  curves  and  one 
corresponding  nuclear  reference  curves  appear  for  two  instru¬ 
mentation  locations  m  Figures  60  anu  61.  The  ragged  character 
of  the  FAE  curves  is  aue  to  the  sensitivity  of  calculated 
dynamic  pressure  to  small  errors  in  measured  static  ana 
stagnation  pressures,  when  the  two  pressures  are  nearly  equal. 
Although  three  stagnation  pressure  probes  were  used  during  test 
CD-6  the  oscillations  in  dynamic  pressure  at  the  most  distant 
were  so  great  that  dynamic  pressure  at  this  location  has  been 
excluded  from  final  data  presentation,  except  for  the  peak 
value  (Figure  53) .  A  fourth  stagnation  pressure  probe  was  not 
functional  at  the  time  CD-6  was  conducted. 

5.3  THEORETICAL  INVESTIGATIONS 

A  centrally  initiated  fuel-air  cloud  can  be  exploded 
either  as  deflagration  or  as  a  detonation.  In  the  case  of  a 
defiagrative  explosion,  compression  waves  induced  by  the 
accelerating  flame  front  ao  not  fully  coalesce  into  a  shock  for 
an  appreciable  distance  from  tne  cloud  boundary.  The  explosion 
does  not  assume  a  blastwave  character  until  the  shock  over¬ 
pressure  is  quite  low  (typically,  'v  7  KPa  (1  psi)).  On  the 
otner  hand,  if  the  cloud  is  detonated,  transition  to  blastwave 
benavior  occurs  at  much  higner  pressures.  Detonative  explosion 
is,  therefore,  the  preferred  mode  if  the  range  of  the  simulated 
nuclear  far-field  airolast  is  to  be  as  wide  as  possiole. 

When  tne  fuel  is  dispersed  uniformly  throughout  the  cloud 
and  initiated  with  an  adequately  strong  shock  source,  a 
spherical  (or  hemispherical)  detonation  having  constant  radial 
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propagation  velocity  aR  /dt  =  U  will  quickly  develop.  The 
front  velocity  corresponds  closely  to  the  well-known  Chapman- 
Jouget  condition.  The  structure  of  the  flowfield  between  the 
explosion  center  and  the  detonation  front  is  approximately 
self-similar  as  the  front  expands. 

Upon  crossing  the  interface  at  the  cloud  boundary,  the 
shock  front  strength  drops  somewhat  as  a  result  of  the  change 
in  acoustic  impedance.  A  U-shaped  expansion-wave  pressure 
profile  forms  near  the  shock  front,  centered  behind  the 
cloud/air  interface.  A  part  of  this  expansion  wave  moves 
inward  toward  the  explosion  center.  The  other  part  interacts 
with  the  leading  shock.  The  shock  decays  because  it  is  no 
longer  protected  from  this  interior  rarefaction  by  a  sonic 
barrier  at  the  Chaproan-Jouget  point.  The  inward  moving  wave 
eventually  reflects  at  the  explosion  center  causing  a  dramatic, 
momentary  drop  in  pressure  there.  Combustion  products  then 
flow  inward  toward  the  center.  The  pressure  again  rises  at  the 
center  until  a  weak  compression  wave,  which  later  steepens  to  a 
weak  shock,  emerges.  Somewhat  after  this  occurs,  the  positive- 
phase  pressure  profile  behind  the  front  and  the  front  decay 
rate  begin  to  develop  point  blastwave  characteristics.  Further 
wave  interactions  within  the  cloud  remnant  continue  beyond  this 
time  but  with  diminishing  magnitude.  A  considerable  residual 
amount  of  the  energy  released  by  combustion  remains  at  later 
times  in  the  expanding  burned  cloud,  primarily  in  the  form  of 
internal  energy.  In  an  equivalent  central  region,  tnis  is  also 
true  of  ideal  point  blastwaves,  but  rough  calculations  have 
indicated  that  in  fuel-air  explosions  the  fractional  amount  of 
residual  energy  may  be  relatively  greater. 

A  calculation  has  been  carried  out  of  the  airblast  that 
develops  from  a  detonated  hemispherical  fuel/air  cloud.  This 
was  done  in  order  to  assess  the  magnitude  of  secondary  shocks 
from  these  explosions,  to  investigate  the  initial  curvature  of 


the  pressure-time  aecay  at  fixed  raaii,  to  provide  a  means  to 
test  the  experimental  data  reduction  progr  and  to  determine 

theoretical  values  for  n  and  s. 

F 

The  calculation  was  executed  in  three  st.  s.  First,  the 
detonation  front  characteristics  were  determined.  These  were 
then  used  as  boundary  conditions  to  compute  the  self-similar 
detonation  wave  interior  profiles.  Finally,  the  self-similar 
detonation  profiles  at  the  time  the  front  reaches  tne  cloud 
boundary  were  used  as  initial  conditions  in  a  calculation  of 
the  subsequent  time-varying  biastwave  flowfield. 

5.3.1  Detonation  Front  Characteristics 

The  properties  of  Chapman -Jouget  detonation  waves  were 
calculated  using  the  Systems,  Science  and  Software's  ORAKL 
code.  This  program  is  a  modification  of  the  TIGER  code  written 
for  SRL  by  SRI  (Reference  64).  Cnemical  equilibrium  is  assumed 
at  cne  C-J  point.  The  program  calculates  the  detonation  front 
velocity,  the  nroouct  composition,  the  thermal  and  thermo¬ 
dynamic  product  mixture  properties,  ana  the  product  velocity. 
Some  of  the  calculated  characteristics  of  propylene  oxide  and 
heptane  detonations  in  air  are  plotted  in  Figures  62-67  as 
functions  of  equivalence  ratio  <j> .  This  is  the  ratio  of  the 
actual  fuel/air  mass  ratio,  4>A,  to  the  stoicniometr ic  fuel/air 
mass  ratio,  <j>A;  i.e., 

*  (57) 

The  calculations  were  carried  out  for  the  two-phase  case?  that 
is,  initially  the  fuel  was  disposed  as  an  unevaporated  spray  of 
liquid  droplets.  Chemical  species  considered  in  the  product 
gases  include  H,  H2,  o,  02,  H20,  OH,  H02,  CO,  C02,  n2,  NO, 


Figure  64.  Absolute  static  pressure  behind  unconfined  two-phase  heptane  and 
Dropylene-oxide  detonations  in  air  . 


Figure  66 .  Product  mixture  specific-heats  ratio  behind  unconfined  two-phase 
heptane  and  propylene  oxide  detonations  in  air. 


5.3.2  Initial  Conditions 


For  most  purposes  it  is  not  necessary  to  numerically 
evaluate  the  fiowfield  behind  the  aetonation  front  throughout 
the  period  in  which  it  propagates  from  the  explosion  center  to 
the  cloud  boundary.  During  this  period,  the  internal  flowfield 
(expansion  wave)  is  usually  self-similar  to  a  reasonable  aegree 
of  approximation.  Therefore,  the  self-similar  profiles  need  to 
be  calculated  only  once.  This  calculation  was  carried  out 
following  the  development  in  Reference  65  which  is  summarized 
here . 


The  physical  pressure,  density,  and  gas  velocity  are 
written  in  terms  of  dimensionless  quantities  II,  F,  and  ft  by  the 
definitions 


?  =  p1X2Ugn(X)  (58) 

P  =  Pj_r  (X)  (59) 

u  =  xusn(X)  (60) 

where 

X  =  R/Rs  .  (61) 


The  detonation  front  velocity  is  Ug  and  its  instantaneous 
position  is  Rg.  A  dimensionless  sound  speed  is  also  defined 
by 

yn 

2  =  Y*  .  (62) 

The  Euler  equations  of  mass  and  momentum  conservation  can  then 
be  reduced  to 


dz  „ 
d« 


z  [2  (ft-1)  2  +  (w-1)  (Y2-l)  n(n-l) 
nt(sw)2  -  uz] 


-  2z] 


(63) 
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dUnA)  _  2  -  (fl-1) 


dH 


(£1-1)  ^  -  0)2] 


(64) 


where  w  =  1,2,3  for  planar,  cylinorical,  and  spherical  (or  hemi¬ 
spherical)  geometries,  respectively.  Also,  since  the  ilowfield 
in  the  rarefaction  behind  the  detonation  front  is  assumed  to  be 
isentropic,  the  entropy  equation,  derived  from  energy 
conservation,  integrates  to 


J2 


*  constant 


or  in  terms  of  tne  dimensionless  variables, 


A2z 


r(r2"1)  r2(Y2-i> 


(65) 


In  these  equations,  the  subscript  *'2"  refers  to  conditions  at  R 
55  R_  (A  =  1) .  This  is  the  C-J  point  and  so 

b 


y2+i 


f  2 

1  +  — 

(  Y-! 


-1 


(66) 


and 


Y2  /  z,> 

Z2 - “  2  +  y* 

Z  (Y2+1)Z\ 


(67) 


where  the  subscript  "1"  refers  to  conditions  in  the  unburned 
fuel-air  mixture,  in  particular, 


-2 

becomes  z.  =  M  '  since  R  "  R  at  the  front  and  R  -  U  t. 

x  s  S  S  w 

To  generate  the  self-similar  profiles,  the  detonation 
properties  computed  previously  were  used  as  input.  First,  z^ 
was  calculated  from  Eq.  (68),  from  Eq.  (67),  and  f 2  from 
Eq.  (66) .  Next,  the  value  of  ft2  was  f°und  from 


(69) 


which  expresses  the  C-J  condition  at  point  2;  i.e.,  the  gas 
velocity  relative  to  the  front  is  sonic  at  the  end  of  the 
reaction  zone.  Now  Eqs.  (63)  and  (64)  are  integrated  starting 
from  the  point  (z2,  ft  2)  and  going  through  decreasing  values 
of  ft  (with  increasing  values  of  z)  to  the  point  (z  =  1,  ft  »  0) . 
The  value  of  y  used  in  tms  integration  is  y2*  Note  that 


at  (z 2,ft2)  £°r  w  =  2,3,  and  that  this  derivative  is  indeter¬ 
minate  at  that  point  w  =  1.  This  implies  ou/dR  -*■  «  at 

(z2,ft2)  for  m  =  2,3  3,  note  that  the  point  (z  =  l,ft  =  0) 

does  not  occur  at  the  ^Aplosion  center  (center  of  symmetry) , 
but  rather  at  X  ^  0.5.  For  values  of  X  less  than  this,  all 
properties  of  the  flow  are  constant  and  u  =  0.  There  is  a 
discontinuity  in  the  slope  of  all  variables  at  (z  =  l,ft  =  0). 

Integration  of  Eqs.  (63)  and  (64)  effectively  yielas 
z  =  z(X)  through  quadrature.  Also,  from  Eq.  (64),  ft  =  ft(X)  is 
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obtained.  Then  F  =  r (X )  can  be  calculated  from  Eq.  (65)  and 
II  *  II  (X )  can  be  obtained  from  Eq.  (62)  witn  y  -  Y  2'  ^  f 

r 1  and  ft  known  as  functions  of  X,  the  physical  variables  P,  p, 
and  u  can  be  determined  as  functions  of  X  by  means  of  Eqs. 

(58) -(60).  Finally,  when  a  value  for  R  is  specified,  these 
profiles  can  be  expressed  in  terms  of  physical  radius  R  in  view 
of  Eq.  (61)  . 

The  self-similar  profiles  of  p,  p,  and  u  behind  a  stoich¬ 
iometric  heptane-air  detonation  at  R  =  4.42  m  (14.5  ft)  are 
plotted  in  Figures  68-70.  This  fuel  and  radius  were  chosen 
because  it  is  expected  in  forthcoming  experimental  FAE  work  to 
test  neptane  clouds  of  that  size.  The  time  for  the  detonation 
to  reach  R  =  4.42  m  (14.5  ft)  is  2.47  ms.  The  self-similar 
profiles  at  this  time  were  used  as  initial  conditions  for  the 
solution  of  the  subsequent  airblast. 

Two  observations  regarding  these  initial  conditions 
should  be  noted.  First,  tne  peak  values  for  the  variables  as 
shown  in  Figures  68-70  are  actually  those  which  obtain  at  the 
end  of  the  detonation  reaction  zone,  that  is,  at  the  Chapman- 
Jouget  point.  In  reality,  however,  the  peak  pressure,  density, 
and  particle  velocity  in  a  detonation  occur  just  behind  the 
leading  shock,  before  any  reaction  is  initiated.  The  actual 
peak  pressure,  tor  example,  is  about  twice  that  shown  in  Figure 
68.  This  is  usually  referred  to  as  the  von  Neumann  spike.  The 
flowtieid  between  the  von  Neumann  spike  ana  the  C-J  point  can 
generally  be  ignored  when  the  reaction  zone  width  is  small 
compared  to  tne  overall  wave  radius,  including  the  expansion 
wave  (i.e,,  R  ) .  The  impulse  imparted  by  this  very  small 

w 

segment  of  the  overall  wave  is  negligible  in  that  case.  For 
that  reason,  tne  von  Neumann  spike  was  not  included  in  the 
detonation  profiles.  However,  in  terms  of  initiating  the  air- 
blast  calculation,  it  should  be  noted  that  the  changes  in 
(P,p,u)  across  the  detonation  front  as  shown  in  Figures  68-70 
are  not  consistent  with  those  of  any  simple  shock,  wave.  Thus, 
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in  the  first  moment  following  t  =  2.47  ms  a  perturbation  was 
expected  in  the  numerical  airblast  calculation.  This  did  in 
fact  occur,  but  the  perturbation  was  rapidly  damped. 

The  second  comment  concerns  a  related  observation.  The 
flowfiela  comprising  a  detonation  and  its  trailing  expansion 
wave  cannot  actually  be  exactly  self-similar  except  in  the 
limiting  case  of  vanishing  reaction  zone  width.  For  non-z  o 
reaction  zone  widths,  a  certain  quantity  of  mass  is  carried  at 
all  times  within  the  reaction  zone.  Tnis  mass  is  effectively 
removed  from  the  expansion  wave  flow  which  must,  therefore, 
have  a  density  profile  that  differs  from  the  self-similar 
profile.  That  is,  the  expansion  wave  must  be  stronger  than  the 
self-similar  wave  (more  rapid  density  drop-off)  as  the  reaction 
zone  width  increases.  Normally,  in  gas-phase  detonations  the 
amount  of  mass  in  the  reaction  zone  can  be  ignored  except  at 
very  small  radii.  However,  two-phase  detonations  in  sprays  of 
2  mm  (0.08  in.)  droplets  (for  example)  have  reaction  zones  on 
the  oroer  of  0.5  m  (20  in.)  in  length.  In  that  case,  tne 
amount  of  mass  in  tne  reaction  zone  is  not  negligible  wnen  the 
wave  radius  is  4.42  m  (14.5  ft). 

This  observation  may  be  of  some  significance  to  the 
detonation  initiation  question.  Presumably,  there  must  exist  a 
minimum  wave  radius  below  which  it  wouio  not  be  possible  to 
establish  a  C-J  detonation.  This  is  Decause  tnere  would  not  be 
enough  mass  at  smaller  radii  to  fill  the  reaction  zone  to  C-J 
levels  even  if  a  vacuum  existed  throughout  the  region  from  the 
explosion  center  to  the  reaction  zone.  To  illustrate  this 
point  it  was  assumed,  for  simplicity,  that  the  density  varia¬ 
tion  between  the  von  Neumann  spike  ano  the  C-J  point  is  linear. 
Values  of  tne  densities  at  these  points  in  a  stoichiometric 
vapor-pnase  heptane-air  detonation  are  p  =  6.25  kg/m"^ 

(0.39i  lbm/tt^)  and  p9  =  2.23  kg/ni"^  (0.139  ibm/ft^).  The 

^  3 

upstream  mixture  density  is  =  1.23  kg/m  (0.0769 
lbm/ft^) .  it  is  readily  shown  by  mass  conservation  that  the 


detonation  wave  radios  for  which  all  of  the  mass  initially 

between  0  -  R  -  R  would  be  contained  within  the  C-J  reaction 

s  * 

zone  in  this  example  is  R,  =  9.46  <5,  where  5  is  the  reac- 

tion  zone  thickness.  Thus,  if  the  heptane  were  dispersed  as  a 

vapor  so  that  5  ^  5  mm  (0.2  in.),  tnen  a  Chapman-Jouget  det- 

* 

onation  could  not  exist  at  radii  smaller  than  R  ^  47  mm 

s 

(1.9  in.) . 

At  the  other  extreme,  if  the  heptane  were  dispersed  as  a 

fairly  coarse  spray  such  that  6  ^  0.5  m  (20  in.),  then  p-^ 

=  1.25  kg/m3  (0.0781  lbm/ft3),  P2s  ■  6.01  kg/m3  (0.376  lbm/ft3), 

and  P_  =  2.25  kg/m3  (0.14  lbm/ft3)  .  a  Chapman-Jouget.  deton- 
^  * 
ation  could  not,  in  tnis  case,  be  achieved  within  R  ^4.5  m 

s 

(i5  ft) .  It  is  noted  that  in  this  last  case,  h*  would  exceed 

s 

the  experimental  cloud  size.  The  structure  within  the  expan¬ 
sion  wave  of  a  vapor-pnase  detonation  would  be  well  represented 
by  the  self-similar  profiles  of  Figures  68-70  at  a  radius  of 
4.42  m  (14.5  ft),  on  the  other  hano,  very  significant  depar¬ 
tures  from  seif -similar lty  would  be  expected  in  the  two-phase 
case.  The  impact  of  these  departures  on  the  airblast  generated 
is  not  presently  Known.  However,  the  reduction  in  the  rare¬ 
faction  wave  pressures  cind  densities  would  seem  certain  to 
influence  cne  airblast  development.  It  seems  clear  that  to 
produce  truly  self-similar  detonation  waves  in  two-phase  media, 
very  fine  droplets  are  required  if  the  clouds  are  of  small 
scale.  In  clouas  at  are  very  large  (R,  >>  R*)  ,  the  effect 

of  drop  size  would  b  1,  but  the  airblast  would  differ  from 

* 

that  generated  by  a  s..  cloud  (R  ^  R  ). 

C  s 


5.3.3  Airblast  Calculation 

The  airblast  that  develops  from  the  detonation  was  calcu¬ 
lated  using  the  S3  SKIPPER  code  (Reference  66) .  This  is  a 
one-dimensional  Layrangian  code  operable  in  any  of  the  three 
geometries.  Tne  code  was  executed  for  t.  present  calculation 
u^ing  quadratic  artificial  viscosity.  No  linear  artificial 
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viscosity  was  required.  Shock  waves  thereby  retained  fairly 
steep  slopes  in  the  solution.  The  calculation  was  begun  at  R 

S 

*  4.42  m,  t  =  2.57  ms  using  the  self-similar  detonation 
profiles  as  initial  conaitions.  It  was  stopped  when  the 
positive  phase  ended  at  the  position  ft  =  30  m  (98.4  ft).  This 
is  the  position  of  the  most  distant  instrumentation  trough  on 
the  experimental  test  pad.  The  peak  overpressure  at  that 
location  is  about  35  kPa  (5  psi) . 

The  calculated  pressure  profiles  within  the  airblast  at 
various  stages  of  its  decay  are  presented  in  Figures  71  and  72. 
The  profiles  were  plotted  at  times  that  correspond  to  the 
arrival  of  the  front  at  radii  equal  to  the  experimental  trans¬ 
ducer  locations.  The  times  shown  in  the  figures  are  times  of 
arrival.  Note  the  change  in  scale  between  figures.  Beyond 
about  three  cloud  radii  the  positive-phase  pressure  profile 
assumes  the  general  character  of  a  point. source  blastwave. 
Secondary  and  tertiary  shocks  and  expansions  are  confined  to 
tne  negative  phase. 

The  static  overpressures  experienced  at  each  transducer 
position  are  plotted  as  functions  of  time  in  Figures  73-78. 
Again,  it  can  be  seen  that  secondary  shocks  are  confined  to  the 
negative  phase.  Also,  the  curvature  of  the  initial  pressure 
drop-off  is  positive  even  close  to  the  cloud.  This  is 
important  in  simulating  nuclear  airblasts  because  the  way  in 
which  impulse  develops  in  the  early  period  of  pressure  decay 
may  affect  damage  mechanisms  in  some  cases. 

The  path  of  the  interface  that  separates  combustion 
products  from  the  surrounding  air  is  plotted  in  Figure  79.  The 
detonated  cloud  first  expands  to  roughly  twice  its  initial 
radius  and  then  contracts  slightly.  The  outward  travel  of  the 
interface  extends  to  a  location  at  which  the  peak  overpressure 
is  about  0.283  MPa  (41  psi). 


Static  Overpressure 


Figure  73.  Calculated  static  overpressure  versus  time  at 
fixed  range  from  hemispherical  FAE  airblast. 
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Calculated  static  overpressure  versus  time  at 
fixed  range  from  hemispherical  FAE  airblast. 
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Figure  77.  Calculated  static  overpressure  versus  time  at 
fixed  range  from  hemispherical  FAE  airblast. 
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The  calculated  results  were  evaluated  further  by  treating 
them  as  though  they  had  been  acquired  as  experimental  data. 
Stagnation  probe  pressures  were  determined  from  the  computed 
flow  Mach  numbers  and  static  pressures  using  Eq.  (52).  The 
theoretical  data  were  then  input  to  the  data  reduction  program 
that  was  described  earlier.  The  results  appear  on  the  next 
several  figures. 


It  was  found  that  Ro„  -  values  that  were  determined 

from  matching  overpressure  positive-phase  duration,  t"T, 

fc' 

differed  widely  from  those  inferred  from  peak  overpressures  and 
from  overpressure  impulses.  It  was  decided  as  an  interim 


measure  to  use  only  impulse  and  peak  overpressure  to  determine 

Rog  F  until  somewhat  greater  confidence  in  the  nuclear 

reference  curve  for  t*  is  established.  The  values  of  R0- 

p  b  ,  r 

are  plotted  at  each  transducer  position  in  Figure  80.  The 


average  of  all  these  values  is  R0£  F  =  9.83  m  (32.2  ft,. 
This  value  was  used  in  the  remaining  data  reduction.  The 


corresponding  scale  factor  and  effective  blast  coupling 

efficiencies  were  found  to  be  s  -  15.2,  h_  =  0.867. 

F 


Peak  static  overpressure  and  total  static  impulse  are 

presented  in  Figures  81  and  82.  The  peak  pressures  generally 

fall  below  and  the  impulses  generally  above  the  nuclear 

reference  curves.  This  is  aue  to  the  choice  of  average  R  . 

o 

Agreement  with  the  nuclear  curves  improves  with  radius  up  to 
the  point  at  which  the  calculation  was  stopped.  This  was  also 
the  case  with  peak  dynamic  pressure  and  with  overpressure 
positive  phase  duration,  which  appear  in  Figures  83  and  84. 


in  the  next  series  of  figures,  Figures  85-96,  the 
dimensionleos  static  overpressures  and  dynamic  pressures  are 
plotted  as  functions  of  dimensionless  time  and  compared 
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Peak  Static  Overpressure,  &P/P 


Figure  81.  Calculated  peak  static  overpressure  from  FAE 

versus  range,  compared  against  nuclear  reference 
curve . 
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Figure  87.  Calculated  neptane-air  FAE  static  overpressure 
versus  time  after  shock  arrival  compared  with 
nuclear  reference  curve,  at  R/R  =  1.42. 


184 


Static  Overpressure  (AP/P 


NUCLEAR 


Time  (t-ta)a^R 

Figure  90.  Calculated  heptane-air  static  overpressure  versu 
time  after  shock  arrival  compared  with  nuclear 
reference  curve,  at  R/R  =  3.04. 
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against  the  fully  smootnea  nuclear  curves.  The  theoretical  FAE 
aata  have  also  been  fully  smoothed.  The  agreement  between  the 
FAE  and  nuclear  profiles  Decoines  quite  reasonable  at  ranges 
corresponding  to  overpressures  below  about  70  kPa  (10  psi) . 
Closer  to  the  explosion  center  the  FAE  overpressure  curves  are 
flatter  than  the  respective  nuclear  profiles.  That  is,  the  FAE 
peak  overpressures  are  lower  and  the  overpressure  positive- 
phase  durations  are  longer.  The  greater  durations  result  in 
impulses  that  are  above  the  nuclear  values  in  spite  of  the 
lower  peak  overpressures. 

Finally,  it  should  be  noted  tnat  two  simplifications  in 
the  FAE  calculation  may  have  a  significant  impact  on  these 
results.  In  the  calculation,  the  chemical  composition  of  the 
cloud  was  assumed  to  be  invariant  with  time,  and  the  specific 
heats  withir  the  cloud  ana  in  the  air,  while  not  equal,  were 
assumed  to  be  constants.  The  chemical  composiion  of  the  cloua 
was  taken  to  be  tnat  calculated  oy  the  GRAKL  code  for  the 
aetonation  products.  The  fixed  specific  heats  (ana  y )  in  the 
cloua  were  taken  to  correspona  to  these  proaucts. 

In  reality,  the  substantial  pressure  and  temperature 
excursions  in  the  exploued  FAE  cloua  would  result  in 
continuously  shifting  local  chemical  equilibrium.  This,  in 
itself,  would  alter  local  specific  heats  and  would  result  in 
additional  local  heat  release  as  dissociated  species  recom¬ 
bine.  Tne  specific  heats,  both  within  tne  cloud  and  in  the 
surrounding  air,  would  also  change  as  the  local  temperatures 
change . 

The  accuracy  of  gas-phase  detonation  calculations  has 
long  Deen  known  to  be  strongly  dependent  on  proauct  composition 
assumptions  (Reference  67).  This  is  due  to  tne  considerable 
dissociation  that  occurs  at  the  very  high  temperatures  attained 
in  detonations.  The  initial  shock  compression  of  air  surround¬ 
ing  an  FAE  cloud  also  produces  large  temperature  increases 
during  the  early  stages  of  blastwave  aevelopment.  Subsequent 


expansion  then  reduces  these  temperatures  markedly.  Therefore, 
a  dissociation/reassociation  cycle  can  be  expected  in  the  air 
outside  the  cloud  as  well  as  in  the  aetonation  produces.  A 
more  accurate  FAE  calculation  would  include  local  chemical 
equilibrium  and  variable  specific  heats,  both  within  the  cloud 
and  in  the  surrounding  air. 
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SECTION  6 


CONCLUDING  REMARKS 


A  reusable  1/4-ton  nuclear  equivalent  FAE  airblast 
simulator  has  been  built.  Symmetrical  fuel-air  clouds  of 
propylene  oxide  and  heptane  have  been  formed  impul  sively  from 
clusters  of  600-1400  nozzles  located  at  the  cloud  center. 

These  clouds  have  been  routinely  detonated  and  the  scaled 
blast  waveforms  are  in  reasonably  good  agreement  with  nuclear 
waveforms  at  the  scaled  ranges.  Theoretical  calculations  have 
yielded  similar  results.  Results  fron.  survey  experiments  with 
large-diameter  single  nozzles  are  preliminary  but  indicate  that 
it  should  be  possible  to  form  clouds  large  enough  to  simulate 
the  far-field  airblast  of  a  1  KT  surface  nuclear  explosion. 

Scaleup  to  the  1  KT  level  involves  many  issues.  Theo¬ 
retically,  the  scale  of  an  FAE  cloud  should  not  affect  its 
airblast  characteristics  unless  the  detonation  reaction  zone 
contains  a  significant  proportion  of  the  total  cloud  moss 
(when  the  det  ^nation  reaches  the  cloud  boundary) .  This  is  a 
possible  problem  at  the  1/4-ton  scale  with  coarsely  atomized 
liquid  jets.  The  two-phase  detonation  in  these  small  clouds 
may  possess  expansion-wave  flowfields  that  are  very  different 
from  the  self-similar  flowfields  in  larger  detonated  clouds. 

As  a  result,  the  effective  airblast  coupling  efficiency  may 
differ  with  scale. 

Beyond  this  effect,  however,  fuel  droplet  size  should 
have  little  impact  on  the  fuel-air  detonation  or  on  the  air- 
blast  that  is  produced  by  it,  provided  that  initiation  is 
possible.  Initiation  is  a  relatively  greater  problem  in 
small-scale  than  in  large-scale  clouds.  The  droplets  of  the 
atomized  jets  do  not  scale  linearly  with  cloud  size,  although 
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somewhat  larger  mean  initial  droplet  diameters  are  expected 
from  the  breakup  of  .larger  jets.  Therefore  the  initiation 
requirements  of  large  clouds  should  not  be  far  removed  from 
those  for  much  smaller  clouds.  If  the  fuel  were  allowed  to 
vaporize  the  initiation  requirements  would  be  identical.  The 
total  energy  release  of  the  initiator  and  its  physical  size 
can  be  made  negligible  in  both  cases  but  this  is  especially 
so  in  the  case  of  large  clouds. 

The  size  Df  the  fuel-air  cloud  and  the  minimum  quantity 
of  fuel  needed  to  simulate  a  nuclear  explosion  of  specified 
yield  is  determined  by  the  effective  blast  efficiency.  This 
efficiency  was  found  experimentally  to  be  about  74%  for 
propylene  oxide  FAE  clouds .  Theoretically  it  was  found  to  be 
87%  for  a  heptane  cloud.  Experiments  with  small  scale  heptane- 
air  clouds  have  yielded  efficiencies  of  about  40%.  These  values 
are  defined  with  respect  to  the  constant-pressure  heat  of  combus¬ 
tion  which  is  used  as  a  reference.  It  should  be  emphasized  that 
the  effective  airblast  efficiency  is  in  essence  determined  by  a 
process  of  best  fit  matching  of  theoretical  or  experimental  FAE 
data  with  the  respective  nuclear  curvefits.  The  values  that  are 
obtained  in  this  way  are  very  sensitive  to  the  nuclear  curvefit 
choice.  For  example,  the  nuclear  curvefits  used  presently  are 
those  of  Reference  52  (Brode,  1964)  .  If  instead  the  curvefits  of 
Reference  51  (Brode,  1970)  are  used,  the  effective  blast  effi¬ 
ciency  of  the  theoretical  heptane  FAE  is  found  to  be  68%  rather 
than  87%. 

Nuclear  airblast  simulation  by  fuel-air  explosions  is 
limited  to  ranges  at  which  the  overpressures  are  below  the 
Chapman- Jouget  detonation  pressure.  As  examples  the  C-J  over- 

/ 

pressures  of  stoichiometric  heptane/air  and  propylene  oxide/air 
detonations  under  standard  conditions  are  1.78  MPa  (258  psi) 
and  1.86  MPa  (270  psi)  respectively.  The  peak  detonation  pres¬ 
sure  increases  with  decreasing  initial  air  temperature  so  that 
simulation  could  be  extended  to  somewhat  higher  overpressures 
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if  tests  are  conducted  during  cold  weather.  It  might  also  be 
feasible  to  increase  the  center-cloud  pressure  by  dispersal 
of  a  liquid  high  explosive  with  the  fuel  spray. 

FAE./nuclear  airblast  simulation  may  also  be  limited  by 
the  extent  of  expansion  of  the  detonated  fuel-air  cloud.  The 
degree  of  this  expansion  was  not:  investigated  experimentally, 
but  in  the  theoretical  heptane  FAE  calculation  the  cloud 
expanded  to  a  radius  at  which  the  peak  overpressure  was  about 
0.283  MPa  (41  psi) .  Testing  would  have  to  be  confined  to 
ranges  having  peak  overpressures  below  this  level,  if  the 
effect  of  hot  combustion  products  on  the  target  being  tested 
were  a  matter  of  concern. 

Ultimate  limitations  on  the  maximum  scale  of  hemisphe¬ 
rical  FAE  clouds  are  not  presently  known.  It  seems  likely 
that  cost  would  in  practice  limit  the  scale.  While  it  appears 
feasible  to  achieve  a  1  KT  simulation,  it  is  not  known  what 
additional  engineering  difficulties  would  have  to  be  overcome 
to  reach,  say,  the  10  KT  level.  One  possible  technique  to 
attain  higher  level  simulation  would  be  to  locate  the  fuel-air 
cloud  in  the  vertex  of  a  large  natural  earth  formation  in  the 
shape  of  a  spherical  wedge.  This  would  reduce  the  size  of  the 
required  dispenser  hardware  but  all  testing  would  have  to  be 
confined  to  regions  inside  the  wedge. 

A  second  alternative  to  reach  larger  scales  is  to 
abandon  the  hemispherical  FAE  cloud  geometry.  The  hemispherical 
shape  is  the  most  efficient  in  terms  of  ground-level  airblast 
coupling.  It  also  minimizes  unwanted  secondary  shocks  and 
permits  point  symmetrical  test  hardware  layout,  maximizing  the 
data  yield  per  test.  If  these  advantages  can  be  compromised, 

FAE  cloud  shaping  could  be  used  to  produce  large-scale  (long 
positive-phase  duration)  simulation  in  localized  spacial 
volumes  of  restricted  size. 
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It  should  also  be  noted  that  in  the  extreme  far-field 
( i 14  KPa  (2  psi)  overpressure)  the  influences  of  cloud  shape 
are  almost  entirely  lost.  Large-diameter,  pancake-shaped 
FAE  clouds  of  virtually  unlimited  yield  could  be  formed  by  a 
dispensing  system  comprising  for  example  a  number  of  point 
source  dispensers  such  as  that  described  in  this  report.  Very 
large-yield  testing  of  aeronautical  systems  in  tne  extreme 
far-field  might  be  accomplished  in  this  way. 

In  all  cases,  FAE  olcuds  formed  from  the  impulsive  dis¬ 
pensing  of  a  large  number  of  liquid  fuel  jets  are  subject  to 
many  constraints,  some  of  which  are  conflicting.  In  most 
respects  dispensing  at  the  highest  practical  jet  velocity  is 
advantageous.  This  improves  atomization  which  enhances  fuel 
distribution  uniformity.  It  also  reduces  the  cloud  formation 
time  which  is  important  in  minimizing  cloud  shape  distortions 
due  to  winds,  and  it  increases  the  jet  spray  column  width  which 
decreases  the  required  number  of  jets  needed  to  fill  in  the 
cloud.  Unfortunately,  high  velocity  injection  leads  to  inter¬ 
actions  between  the  .spray  and  the  air  which  can  severely  dis¬ 
tort  the  cloud.  Whether  or  not  accurate  cloud  shape  is  con¬ 
sidered  a  critical  requirement,  the  volume  of  air  engulfed 
by  these  distorted  clouds  differs  from  that  intended  and  this 
influences  the  effective  FAE  yield.  One  possible  remedy  to  the 
problem  of  distortions  due  to  high  velocity  jets  may  lie  in 
tailored  aiming  of  the  injector  nozzles,  with  some  nozzles 
directed  downward.  The  distortion  can  also  be  reduced  by 
decreasing  the  number  of  nozzles,  but  if  too  few  nozzles  are 
used  the  jet  sprays  will  not  overlap  at  the  outermost  portions 
of  the  cloud.  Overlap  in  this  region  is  necessary  since  a 
large  fraction  of  the  total  cloud  mass  is  contained  in  it. 

It  appears  that  1400-lb00  nozzles  are  needed  for  complete 
overlap  in  a  cloud  produced  by  point  source  injection. 
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While  it  is  probably  not  necessary  to  fully  understand 
all  details  of  jet  breakup  and  cloud  formation  in  order  to 
build  a  successful  FAE/nuclear  simulator,  basic  engineering 
data  relating  to  certain  rudimentary  features  are  essential. 

The  reach  and  width  of  individual  spray  columns  formed  under 
conditions  that  produce  adequate  atomization  are  needed  as  is 
the  formation  time.  Reach,  width,  and  formation  time  are  in 
general  functions  of  injection  velocity,  nozzle  diameter,  fuel 
properties,  and,  it  appear's,  the  quantity  of  fuel  dispensed. 

For  each  fuel  and  nozzle  size  there  is  a  minimum  injection 
velocity  below  which  atomization  is  poor.  The  minimum  appears 
to  increase  with  nozzle  diameter.  Above  this  velocity  atomiza¬ 
tion  continuously  improves  while  the  jet  reach  changes  only 
slightly.  Increasing  the  nozzle  diameter  increases  bov.h  reach 
and  spray  column  width.  Jet  reaches  are  typically  1000-2000 
nozzle  diameters,  while  jet  spray  widths  are  typically  100 
nozzle  diameters.  Ic  may  be  possible  to  increase  jet  reach 
through  the  use  of  additives  such  as  long-chain  polymers  but  it 
is  anticipated  that  this  would  adversely  affect  atomization, 
particularly  near  the  nozzle.  Fuel  distribution  in  individual 
jet  sprays  is  unknown  but  the  spray  density  appears  visually  to 
increase  with  distance  along  the  jet.  This  should  permit 
rather  uniform  clouds  to  be  formed  when  a  large  number  of  such 
jets  overlap  after  point  source  injection. 
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One  feature  of  jet  breakup  that  has  emerged  during  the 
survey  experiments  is  the  influence  on  jet  reach  of  jet  tail 
breakup.  When  the  end  of  the  jet  emerges  from  a  nozzle, 
material  from  the  jet  tail  appears  to  continuously  be  torn 
off.  It  has  been  postulated  that  this  phenomenon,  which  /aav  be 
due  to  vortex  shedding  at  the  tail,  will  reduce  the  jet  reach 
if  an  insufficient  quantity  of  fuel  is  dispensed.  Some 
experimental  evidence  has  been  generated  in  support  of  this 
postulate  which,  however,  is  presently  considered  unproven. 

If  the  effect  is  real  it  could  have  a  major  impact  on  dispen¬ 
ser  design.  The  implication  would  e  that  there  is  a  minimum 
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quantity  of  fuel  that  must  be  dispensed  through  each  nozzle 
in  order  for  the  jet  to  attain  full  reach.  Thus  the  minimum 
quantity  of  fuel  needed  to  form  a  fuel-air  cloud  of  specified 
radius  would  equal  the  number  of  nozzles  needed  multiplied  by 
the  minimum  quantity  of  fuel  per  nozzle.  This  could  well 
exceed  the  quantity  of  fuel  required  calculated  on  the  basis 
of  effective  airblast  efficiency. 

It  is  recommended  that  systematic  experiments  with 
single  jets  be  conducted  in  which  jet  velocity,  nozzle  diam¬ 
eter,  and  fuel  quantity  are  varied  over  a  reasonably  wide 
range.  This  will  require  impulsive  dispensers  that  are  larger 
and  smaller  than  the  existing  U-tube  and  linear  dispensers. 

It  has  been  found  that  each  dispenser  is  limited  to  a  narrow 
range  of  jet  sizes  that  can  be  tested.  The  largest  jet  that 
can  be  produced  is  limited  by  the  fuel  capacity  of  the  dis¬ 
penser.  The  smallest  jet  is  limited  by  flow  non-uniformities 
when  too  little  fuel  is  dispensed.  In  the  latter  case  multiple 
nozzles  can  be  attached  to  the  dispenser  in  a  single  test  to 
increase  the  quantity  of  fuel  dispensed,  provided  the  layout  is 
such  that  the  jets  do  not  interfere  with  one  another. 

It  is  also  recommended  that  FAE  tests  using  the  1/4- 
ton  reusable  facility  be  continued  to  generate  a  better  data 
base  for  comparison  with  nuclear  explosions.  Additional  areas 
that  must  be  addressed  before  a  full-scale  FAE  blast  simulator 
can  be  constructed  are  more  engineering-oriented  and  include 
hardware  configuration,  fuel  handling  techniques,  repeatability, 
safety,  and  construction  and  life-cycle  costs. 
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APPENDIX  A 


EXPERIMENTAL  FACILITY 


Preliminary  experimental  work  performed  at  Systems, 
Science  and  Software  (Reference  2)  examined  the  basic  tech¬ 
nological  feasibility  of  a  point-source  dissemination,  FAE 
airblast-simulation  concept.  For  that  purpose  relatively 
coarse  measurements  were  acceptable  and  an  undeveloped  test 
area  adjacent  to  an  existing  bunker  was  used  for  the  experi¬ 
ments.  Although  reasonably  level,  this  test  area  had  not  been 
fully  graded  or  surfaced.  Strain  gage  transducers,  buried 
directly  in  the  soil  were  used  to  measure  static  pressures. 

This  simple  arrangement  proved  adequate  to  establish 
concept  feasibility.  However,  it  was  felt  that  additional 
refinements  would  be  necessary  in  order  to  generate  more  con¬ 
clusive  data.  A  number  of  facility  improvements  were  accord¬ 
ingly  undertaken.  The  present  facility  is  described  in  this 
section.  The  facility  is  located  at  the  Green  Farm  Test  Site 
which  is  managed  by  Systems,  Science  and  Software  for  the 
Defense  Nuclear  Agency.  Green  Farm  is  a  part  of  the  Camp  Elliot 
Naval  Reservation  neat  San  Diego,  California. 


A. 1  TEST  PAD 

A  concrete  test  pad  was  built  to  provide  a  smooth 
blastwave  path  from  the  explosion  center.  This  was  intended 
to  minimize  measurement  aberations  that  might  otherwise  result 
from  wave  reflections  in  the  vicinity  of  pressure  transducers. 
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A  schematic  of  the  pad  layout  is  shown  in  Figure  A-l  and 
a  photograph  appears  in  Figure  A-2.  The  test  area,  which  has 
the  approximate  overall  dimensions  of  37  m  x  46  m  (120  ft  x 
150  ft) ,  is  graded  level,  smoothed,  and  covered  with  a  2  1/2  cm 
(1  in.)  layer  of  pea  gravel.  Two  perpendicular  steel-reinforced 
concrete  legs  extend  continuously  from  an  octagonal  center  area. 
The  shorter  of  these  two  legs  was  included  to  permit  blastwave 
symmetry  testing. 

Abrupt  slopes  surround  the  test  area.  However,  it  was 
possible  to  find  an  orientation  for  the  concrete  pad  such  that 
the  arrivc.l  of  reflected  waves  from  these  slopes  would 
not  interfere  with  measurements  during  a  20  ms  minimum  period 
following  the  arrival  of  the  blast  front  at  each  location  on 
the  pad .  This  exceeds  the  longest  positive  phase  duration  of 
blastwaves  at  the  experimental  scale  (1/4-ton  nuclear  equiv¬ 
alent)  . 

The  average  thickness  of  the  concrete  pad  is  15  cm 
(6  inches) .  A  rectangular  opening  in  the  center  area  is 
provided  for  emplacement  of  the  fuel  dispenser  in  the  ground. 
Since  rain  water  runoff  from  adjacent  hills  could  flood  the 
test  area,  ditching  is  provided  and  a  sump  pump  system  has 
been  installed. 

Instrumentation  troughs  were  installed  at  several  loca¬ 
tions  along  each  runway,  with  approximate  positional  corres¬ 
pondence  between  the  two  legs.  A  cross-section  of  a  trough 
is  shown  in  Figure  A-3  and  a  photograph  is  given  in  Figure  A-4 . 
In  each  trough,  a  15.2  cm  x  61  cm  x  1.9.  cm  (6  in.  x  24  in. 
x  3/4  in.)  steel  cover  plate  is  securely  fastened  with  six 
recessed  machine  screws.  The  cover  plates  were  individually 
milled  for  each  trough  so  that  their  upper  surfaces  are  flush 
with  the  concrete.  Side-on  transducers  and  stagnation  pressure 
probes  (stings)  are  installed  on  these  plates.  A  transducer 
location  schedule  is  given  in  Table  A.l. 
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m  Figure  A—4.  Photograph  of  stagnation  pressure  sting  mounted 

|  in  instrumentation  trough.  A  static  pressure 

A  transducer  is  flush  mounted  directlv  below  sting 

v  face . 


Table  A.l.  Location  of  transducers  on  FAE 
test  pad  (see  Figure  A.l) . 


DISTANCE  FROM  PAD 


INSTRUMENTATION 

TROUGH  NUMBER 

CENTER 

METERS 

(DISPENSER) 

FEET 

1 

2.91 

9.54 

2 

5.92 

19.40 

3 

9.93 

32.60 

4 

13.90 

45.60 

5 

18.90 

62.00 

6 

23.90 

78.30 

7 

29.90 

98.00 

8 

3.05 

10.00 

9 

6.04 

19.80 

10 

10.10 

33.00 

11 

14.10 

46.10 

The  troughs  are  joined  with  5.1  cm  (2  in.)  ID  conduit 
that  is  inlaid  in  the  concrete  runways.  Coaxial  instrumenta¬ 
tion  cables  are  permanently  installed  in  the  conduit,  with 
two  terminations  (BNC-type)  available  at  each  trough.  The 
troughs  also  have  drain  openings  into  the  rain  ditches. 

Two  additional  conduit  lines  service  the  center  pad 
area.  These  are  buried  in  the  ground  adjacent  to  the  longer 
runway.  One  of  these  carries  the  high  voltage  line  used  with 
the  high  explosive  detonators  (to  initiate  FAE  clouds) .  The 
other  is  used  to  route  auxiliary  110  VAC  power,  intercom,  and 
control  cables.  A  buried  0.7  MPa  (100  psi)  compressed  air 
line  also  services  the  center  pad  area. 

A. 2  CONTROL  AND  INSTRUMENTATION 

The  event  sequence  during  a  particular  test  is  con¬ 
trolled  by  a  set  of  adjustable,  precision  time-delay  relays. 
These  are  part  of  a  control  console  that  also  includes  previ¬ 
sions  for  low  and  high  voltage  detonator  ignition,  a  low- 
voltage  power  supply,  and  a  fiducial  pulse  generator.  A 
photograph  of  this  console  is  given  in  Figure  A. 5. 

Instrumentation  during  an  FAE  test  currently  consists 
of  static  and  stagnation  pressure  history  measurements  at 
various  ranges  from  the  cloud  center.  The  static  transducers 
are  flush  mounted  in  the  trough  cover  plates .  The  stagnation 
pressure  transducers  are  mounted  in  stings  fabricated  from 
1.27  cm  (1/2  in.)  schedule  4C  steel  pipe.  The  stings  are 
76  cm  (30  in.)  high  and  have  a  43  cm  (14  in.)  horizontal 
extension  to  the  transducer  face.  The  stagnation  pressure 
transducers  are  located  directly  above  companion  static 
pressure  transducers  on  each  plate. 


Conditioner 
for  PCB 
Transducers 
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The  static  pressure  transducers  presently  in  use  are 
PCB  Model  102A12.  These  are  high  frequency  piezoelectric 
transducers  with  a  built-in  voltage-follower  amplifier.  The 
rise  time  of  these  transducers  is  ^1  us.  Resonant  ringing 
resulting  from  step  pressure  inputs  is  suppressed  by  built-in 
electronics.  The  head-on  transducers  have  equivalent  per¬ 
formance  specifications.  They  were  supplied  by  PCB  as  Model 
113M49,  including  aerodynamic  sting  termination  mounting.  The 
transducers  are  powered  by  a  modified  twelve-channel,  20  MA 
constant  current  supply,  PCB  Model  483A.  For  low  frequency 
pressure  measurements  a  Validyne  Model  DP-15  variable  reluc¬ 
tance  transducer  is  used  with  a  model  CD12  indicating  signal 
conditioner.  The  output  from  all  transducers  is  recorded  on 
an  EMI  Model  SE7000A  14-channel  tape  deck,  using  Scotch  2.54 
cm  (1  in.)  instrumentation  tape.  At  a  recording  velocity  of 
305  cm/s  (120  in./s) ,  this  instrument  has  a  frequency  response 
to  80  KHz  and  is  the  limiting  component  in  the  data  acquisition 
chain.  When  faster  response  is  needed  signals  are  recorded 
individually  on  oscilloscopes  or  fast-transient  recorders. 

All  transducers  are  periodically  tested  in  a  shock  tube  to 
verify  that  they  are  operating  properly. 

A. 3  PHOTOGRAPHY 

In  addition  to  pressure  measurements,  high  speed  films 
are  taken  during  many  of  the  tests.  For  framing  rates  to 
400  fr/s  a  Milliken  Model  DBM-4A  is  used.  Framing  rates  up 
to  7500  fr/s  are  obtainable  with  a  Fastax  Model  WF45.  However, 
this  highest  speed  requires  a  120  m  (400  ft)  length  of  film. 

For  a  38  m  (125  ft)  roll  of  film,  the  framing  rate  reaches 
3500  fr/s.  Negative  films  are  processed  at  the  test  site 
using  a  Cramer  Mark  I  Cine  Processer. 


The  filming  of  high-reach  vertical  jets  has  presented 
some  difficulty  when  portions  of  the  jets  stand  against  a 
background  of  sky.  Little  success  has  been  realized  in  ob¬ 
taining  hiah  resolution  films  on  cloudy  or  hazy  days,  in  spite 
of  many  attempts  with  various  film,  filter,  and  polarizer  com¬ 
binations.  However,  when  the  sky  is  blue  it  has  been  found 
that  satisfactory  results  can  be  achieved  by  using  a  film 
with  extended  red  sensitivity  in  conjunction  with  a  deep  red 
filter.  Polarizing  filters  were  found  not  to  be  effective  in 
this  application.  At  present,  Kodak  PAR  2479  or  2475  film  and 
a  Tiffen  #29  red  filter  are  used.  The  film  is  processed  for 
five  minutes  in  Kodak  D-19  developer  at  25 °C  (75°F) .  Under 
these  conditions  an  assumed  effective  f : lm  speed  of  ASA  200 
gives  the  necessary  contrast  of  the  jets  against  the  sky. 
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APPENDIX  B 


LINEAR  DISFENSER  GAS  GENERATOR 

The  gas  generator  presently  in  service  with  the  linear 
dispenser  burns  Hercules  Red  Dot  smokeless  powder  in  a  15  cm 
(6  in.)  OD  by  33  cm  (13  in.)  long  breech.  The  size  of  the 
combustion  chamber  is  about  5  cm  (2  in.)  ID  by  18  cm  (7  in.) 
long.  An  end  cap  or  nut  having  a  right-angle,  2.5  cm  (1  in;) 
diameter  nozzle  passage  directs  the  product  gases  downward  into 
the  dispenser  free  volume  (the  gas  generator  is  mounted  hori¬ 
zontally)  .  The  gas  generator  was  fabricated  from  4340  steel, 
heat  treated  to  1.17  GPa  (170,000  psi)  and  Rockwell-C  hardness 
of  38.  An  outline  sketch  is  given  in  Figure  B-l,  and  a  photo¬ 
graph  is  provided  in  Figure  B-2 . 

Smokeless  powder  was  chosen  over  more  exotic  fast¬ 
burning,  stable  propellants  because  of  the  high  cost  of  the 
latter.  However,  the  low  burning  velocity  of  smokeless  powder 
presented  considerable  difficulty.  The  burning  velocity  of 
Red  Dot  is  given  in  Figure  B-3.  Even  at  138  MPa-  (20,000  psi) 
is  is  only  about  10  cm/s  (4  in./s) .  This  means  that  if  the 
breech  combustion  chamber  were  filled  with  this  powder  and 
ignited  at  one  end,  combustion  would  be  complete  in  about  1.75 
seconds.  Since  the  required  fuel  dispensing  times  a..e  on  the 
order  of  a  few  hundred  milliseconds  at  most,  this  lengthy  burn¬ 
ing  time  could  not  be  tolerated.  Therefore,  to  enable  the  use 
of  smokeless  powder  as  a  propellant,  the  breech  is  mounted 
horizontally  and  the  combustion  chamber  is  only  partially 
filled  with  powder.  The  powder  depth  is  then  about  2.5  cm 
(1  in.).  The  igniter,  which  is  mounted  at  one  end,  produces 
a  flame  that  traverses  the  entire  upper  surface  of  the  powder, 
initiating  it  in  many  locations  simultaneously.  A  sketch  of 
the  igniter  appears  in  Figure  B-4.  The  depth  of  the  powder 
then  determines  the  burning  time,  reducing  it  by  about  an  order 
of  magnitude. 
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In  addition,  the  Red  Dot  products  are  prevented  from 
immediately  exhausting  from  the  breach  by  a  1.8  mm  (0.070  in.) 
thick  copper  diaphragm  that  separates  the  combustion  chamber 
from  the  right-angle  exhaust  nozzle.  The  diaphragm  is  designed 
to  burst  after  combustion  is  completed  thus  delaying  the  start 
of  blowdown.  The  higher  combustion  pressures  prior  to  blowdown 
also  improve  the  combustion  efficiency  and  produce  a  cleaner 
product  exhaust. 

The  gas  generator  size  was  determined  from  a  simple 
thermodynamic  analysis  for  the  quantity  of  smokeless  powder 
required.  All  of  the  powder,  mass  n^,  is  assumed  to  first  burn 
to  completion  in  the  breech.  At  that  time  the  temperature  and 
pressure  in  the  breech  are  Tg  and  Pg.  The  diaphragm  then  rup¬ 
tures  and  the  breech  gases  blow  down  into  and  mix  uniformly  with 
the  air  that  is  initially  in  the  free  volume  tank  (mass  mA) . 

This  process  is  adiabatic.  The  final  mixture  contained  in  both 
the  breech  and  free  volume  has  uniform  composition,  and  the 
final  pressure  and  temperature ,  P2  and  T2 ,  are  uniform  through¬ 
out.  For  the  system  comprising  the  breech  plus  free  volume, 
the  process  is  work-free  so  that  the  First  Law  is  simply 


or 
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mAaA2  +  mBeBj  =  Va-l  +  mBeB1 

where  e.  is  the  mass-specific  internal  energy  of  specie  i  and 
states  1  and  2  correspond  to  conditions  just  prior  to  and  after 
breech  blowdown .  Thus 


mA(eA2  -  eA11  +  ”B(eE2  -  eBl>  =  ° 
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or  assuming  that  both  the  air  and  the  breech  gases  behave  as 
calorically  perfect  gases, 


mACV.  (T2  “  Tl)  +  mB,  CV  (T2  TB, J  0 
A  1  B  1 


where  is  the  initial  air  temperature  and  is  the  constant- 
volume  specific  heat  of  specie  i.  Solving  for  T2  gives 


mBlCvBTBi  +  mACVATl 


(B.l) 


The  final  state  equilibrium  pressure  is  simply 


m2RT2 


p2  =  w,v: 


(B .  2) 
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where  m2  =  mA  +  mB,  and  W2  is  the  final  mixture  molecular 
weight.  However,  for  the  mixture, 
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Combining  Eqs.  (B.1)-(B.3)  gives 
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To  calculate  P2,  the  breech  temperature  reached  after  combus 
tion  but  before  blowdown,  Tb^,  must  be  known.  This  can  be 
estimated  by  treating  the  combustion  process  as  an  external 
heat  addition.  The  energy  equation  for  that  process  is 


*B, 


-  E 


B, 


That  is,  the  increase  in  energy  inside  the  breech  above  its 
initial  value  equals  the  heat  added.  On  a  unit  mass  basis, 
this  is  approximated  by 


B. 
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where  Hc  is  the  heat  of  combustion  per  unit  mass  of  smokeless 
powder.  Treating  the  combustion  products  as  calorically 
perfect,  this  becomes 
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or  with  T  =  T,  , 
B0  1 
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(B.  5) 


The  corresponding  pressure  in  the  breech  at  this  time  is 


mB1RTB1 


(B .  6 ) 
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in  which  is  the  volume  of  the  breech  combustion  chamber. 
Finally,  the  mass  of  air  in  the  free  volume  tank  behind  the 
piston,  also  required  in  Eq.  (B.4)  ,  is  simply 


(B .  7) 


where  is  the  initial  free  volume  pressure. 

During  preliminary  design  of  the  dispenser,  a  maximum 

3  3 

free  volume  of  V  ~  0.085  m  (3  ft  )  was  assumed  for  these 

r  v 

estimates.  It  was  hoped  to  achieve  a  maximum  initial  pressure 

of  P~  =  2.8  MPa  (400  psig)  in  this  volume.  Using  data  for  Red 

^  6 
Dot  powder  supplied  by  the  manufacturer  (H  =  5.14  x  10  J/K3 

(2210  btu/lbm) ,  CVg  -  1471  J/(KG-°K)  (0.352  btu/ (lbm-°R) ) ,  and 

Mg  =  24.9),  the  quantity  pf  powder  required  according  to  Eq. 

(B.4)  is  mBl  =  0.16  kg  (0.35  lbm) .  This  calculation  assumes 

Tx  -  300 °K  (540°R)  ,  so  that  TBl  =  3794°K  (6830°R)  from  Eq.  (B.5). 

The  combustion  chamber  volume  required  for  the  breech  was 

-4  3  3 

V3  =  mg-^/Pg  =  2.72  x  10  m  (16.6  in  ),  since  the  bulk  density 

of  Red  Dot  is  p  =  5.88  kg/m^  (0.368  Ibm/ft^).  The  actual 

“  -4  3  3 

breech  volume  as  built  is  approximately  3.6  x  10  m  (22  in  ). 

The  pressure  developed  in  the  breech  with  mB^  =  0.16  kg  (0.35 

lbm)  of  powder  is  PBl  =  56.1  MPa  (81,400  psi)  ,  according  to 

Eq.  (B.6) .  Theoretical  free  volume  pressures  corresponding  to 

other  quantities  of  smokeless  powder  are  plotted  in  Figure  B-5 . 


The  gas  generator  blowdown  time  after  rupture  of  the 
copper  diaphragm  is  a  function  of  the  exhaust  nozzle  orifice 
size.  Since  a  maximum  blowdown  time  of  10-20  ms  was  considered 
requisite  for  impulsive  dispenser  operation  (relative  to  the 
shortest  liquid  dispensing  times  of  about  100  ms) ,  this  would 
suggest  as  large  an  orifice  as  possible.  However  there  was 
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concern  that  with  too  large  an  orifice  the  blowdown  times 
might  be  so  short  as  to  produce  blastwaves  of  damaging 
strength  inside  the  free  volume  tank.  It  was  therefore 
decided  to  limit  the  blowdown  rate . 

To  calculate  the  required  orifice  size,  the  gas- 
generator  blowdown  was  assumed  to  be  isentropic  and  quasi¬ 
steady.  For  a  control  volume  that  coincides  with  the  gas 
generator  surface,  conservation  of  mass  is 


*”b 

dt 


=  -m 


(B .  8) 


where  the  mass  flow  through  the  orifice  of  area  A*  is 


m  = 


P0A* 
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(B .  9 ) 


and  y  is  the  ratio  of  specific  heats  of  the  combustion 

13 

product  gases.  Eq.  (B.9)  assumes  that  the  orifice  is  choked. 
The  breech  pressure  at  any  time  during  blowdown  is  related  to 
the  instantaneous  temperature  and  mass  of  product  gases  by 


P  _  mBRTB 

B  ‘  WBVB 


(B. 10) 


In  addition,  conservation  of  energy  can  be  expressed  by 


where  h*  and  u*  are  the  enthalpy  and  velocity  of  the  products 
when  they  reach  the  orifice.  Since  Eg  =  mgeg,  and  since 

h*  +  u*2/2  =  h£  , 

where  hg  is  the  enthalpy  inside  the  combustion  chamber,  the 
energy  equation  can  be  written  as 

de  ditig 

■b  dT  +  eB  5T  ’  '^B 

Making  use  of  Eq.  (B.8)  ,  and  with  hg  =  eQ  +  P0vB  (where  vR  is 
the  specific  volume  of  the  breech  gases) ,  this  last  equation 
can  be  written  as 


m 


deB  _  • 

B  dt  _mPBVB 


Finally,  assuming  a  calorically  perfect  gas  so  that 


deB  _ 

air  -  cvB 


dTB 

dt 


r 


and  using  PgVg  =  RgTg,  the  energy  equation  can  be  written  as 


dt 


mg  B 


-  1 )  T 


B 


(B.ll) 


Equations  (B.8MB.11)  can  be  integrated  analytically. 
The  solution  for  the  breech  pressure  as  a  function  of  time  was 
found  to  be 


where 


Pg^,  Tg^  are  values  at  t  =  0  (following  combustion  but  just 
prior  to  blowdown) . 

Estimates  of  the  nozzle  orifice  size  required  for  a 
specified  blowdown  time  were  obtained  by  solving  Eq.  (B.12) 
for  A*  in  terras  of  PB  and  t.  The  blowdown  time  was  defined 
as  the  time  at  which  Pg  =  P2.  In  fact  this  equality  occurs 
only  as  t  approaches  infinity.  Eq.  (B.9)  and  therefore  Eq. 

(B.12)  become  invalid  before  Pg  =  P2  because  the  orifice  becomes 
unchoked.  However,  extrapolation  of  Eq.  (B.12)  to  P_  =  P~  does 
provide  a  characteristic  time.  The  values  of  Pg1  and  are 

related  to  P2  and  VFV,  according  to  Eq.  (B.4)  and  (B.6).  With 
TBl  computed  from  Eq.  (B.5)  and  with  yfi  =  1.23,  the  required 
nozzle  orifice  diameter  for  20  ms  blowdown  time  was  calculated 
as  a  function  of  P2  and  Vpy.  The  result  appears  in  Figure  B-6 . 
On  the  basis  of  these  calculations,  a  gas  generator  exit 
orifice  diameter  of  13  mm  (0.5  in.)  was  selected. 

No  attempt  was  made  to  measure  pressures  developed  in¬ 
side  the  breech.  However,  an  approximate  measure  of  the  blow¬ 
down  time  was  obtained  by  observing  the  gas  generator  exhaust 
with  a  photoelectric  cell.  The  exhaust  gases  remained  luminous 
for  about  18  ms  in  most  tests,  in  reasonable  agreement  with  the 
predicted  value  of  blowdown  time. 


1 


2 


FINAL  FREE  VOLUME  PRESSURE,  P2  (MPa) 

Figure  B-6 .  Required  gas  generator  orifice  as  a  function  of 
free  volume  and  pressure  in  the  free  volume,  for 
20  ms  blowdown  time. 
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While  the  gas  generator  performance  was  encouraging  in 
this  respect,  it  was  found  to  be  marginal  in  others.  One  minor 
difficulty  resulted  from  the  absence  of  high  pressure  seals 
between  the  breech  body  and  the  orifice  retainer  nut.  When 
stainless  steel  and  aluminum  diaphragms  were  tested  initially, 
their  vapor  condensed  while  leaking  between  the  threads  of  the 
nut.  The  gas  generator  was  then  very  difficult  to  disassemble. 
This  partd  cular  problem  was  eliminated  simply  by  changing  to 
copper  d  .aphragms . 

A  more  serious  problem  was  experienced  during  operation 
of  the  gas  generator  with  a  full  load  of  Red  Dot  powder  (mg-^ 

=  0.16  kg  (0.352  lbm) ) .  On  separate  occasions  both  the  nut 
and  the  breech  body  failed  and  were  damaged  beyond  repair. 

There  has  been  some  speculation  as  to  whether  in  these  cases 
the  smokeless  powder  may  have  detonated,  although  the  deflagra¬ 
tion  to  detonation  transition  distance  for  Red  Dot  powder  is 
much  greater  than  the  18  cm  (7  in.)  chamber  length,  according 
to  the  powder  manufacturer.  In  any  case,  the  gas  generator 
was  redesigned  to  eliminate  all  sharp  corners  before  being 
rebuilt.  (The  current  version  is  shown  in  Figure  B-2.) 
Subsequent  operation  has  also  been  limited  to  mg  =  0.10  kg 
(0.22  lbm)  of  smokeless  powder. 

The  most  troublesome  problem  with  this  type  of  gas 
generator  is  a  consequence  of  the  relatively  slow  burning  rate 
of  smokeless  powder.  The  total  burning  time  in  the  breech  is 
estimated  as  250  ms.  However  luminous  exhaust  is  detected 
beginning  at  only  80  ms  after  ignition  (which  indicates  dia¬ 
phragm  rupture),  and  the  blowdown  time  is  then  about  20  ms. 
Presumably  a  substantial  amount  of  powder  is  expelled  in  an 
unburned  or  partially  burned  condition.  The  result  is  that 
pressure  levels  actually  reached  in  the  free  volume  tank  are 
about  20-40%  below  the  predicted  values. 
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